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PREFACE 

This  report  is  the  sixth  issue  of  the  Air  Force  Academy  Aeronautics  Digest.*  Our 
policy  is  to  print  articles  which  represent  recent  scholarly  work  by  students  and  faculty 
of  the  Department  of  Aeronautics,  members  of  other  departments  of  the  Academy  and  the 
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USAFA-sponsored  projects,  and  authors  in  fields  of  interest  to  the  USAFA. 

In  addition  to  complete  papers,  the  Digest  also  includes,  when  appropriate,  abstract 
of  lengthier  reports  and  articles  published  in  other  formats.  The  editors  will  consider 
for  publication  contributions  in  the  general  field  of  Aeronautics,  including: 

•Aeronautical  Engineering 
Flight  Mechanics 
Propulsion 
Structures 
Instrumentation 

•  Fluid  Mechanics 

•Thermodynamics  and  Heat  Transfer 

•  Biomechanics 

•Engineering  Education 

•Aeronautical  History 

Papers  on  other  topics  will  be  considered  on  an  individual  basis.  Contributions 
should  be  sent  to: 

Editor,  Aeronautics  Digest 

D  FAN 

US  Air  Force  Academy,  CO  80840 

The  Aeronautics  Digest  is  presently  edited  by  Maj  A.M.  Higgins,  PhD;  Maj  E.J.  Jumper 
PhD;  and  Capt  J.M.  Kempf,  Department  of  English,  who  provided  the  final  editorial  review. 
Our  thanks  also  to  our  Associate  Editor,  Barbara  J.  Gregory,  of  Contract  Technical 
Services,  Inc. 


*  The  first  five  issues  of  the  Digest  can  be  ordered  from  the  Defense  Documentation  Cen¬ 
ter  (DDC) ,  Cameron  Station,  Alexandria,  VA  22324.  Use  the  following  AD  numbers: 
Aeronautics  Digest  -  Spring  1978,  ADA060207;  Aeronautics  Digest  -  Fall  1978,  ADA069044; 
Aeronautics  Digest  -  Spring  1979.  ADA075419;  Aeronautics  Digest  -  Fall  1979,  ADA085770; 
and  Aeronautics  Digest  -  Spring/Summer  1980,  ADA096678. 
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CANARD  WAKE  MEASUREMENT  AND  DESCRIPTION 
G.  Sisson*  and  R.  Crandall** 

Abstract 

This  paper  describes  a  method  for  measuring,  describing,  and  visualizing  the  com¬ 
plex  flow  fields  that  occur  when  canard-configured  aircraft  models  are  tested  in  the 
wind  tunnel.  This  simple,  rapid,  and  inexpensive  technique  is  a  source  of  complex  em¬ 
pirical  aerodynamic  data  suitable  for  the  preliminary  design  of  future  fighter  aircraft. 
The  technique  makes  use  of  a  miniature  seven-hole  probe,  a  computer-driven  traverse  me¬ 
chanism,  computerized  data  acquisition  equipment,  and  computer  graphic  displays.  Numer¬ 
ical  and  graphic  descriptions  of  wakes  from  three  canard  models  are  presented. 

I .  Introduction 

Since  the  summer  of  1978  the  U.S.  Air  Force  Academy  and  NASA's  Ames  Research  Center 
have  participated  in  a  joint  research  effort  to  understand  the  aerodynamics  of  canard- 
equipped  aircraft.  These  configurations  are  similar  to  proposed  fighter  aircraft  designs 
for  the  1990's,  which  have  pitch  control  surfaces  located  in  front  of  the  wings  in  lieu 
of  a  conventional  aft-located  tail.  The  program  began  as  summer  research  projects,  with 
Academy  cadets  and  officers  attempting  to  analytically  predict  the  aerodynamics  of  these 
configurations  using  potential  flow  computer  analysis.  Their  efforts  were  largely  un¬ 
successful  because  the  location  of  the  canard  wake  was  not  known  in  advance.  To  solve 
this  problem,  NASA  Ames  contracted  the  USAF  Academy  Department  of  Aeronautics  to  develop 
a  rapid,  inexpensive  technique  to  precisely  measure  the  location  of  canard  wakes. 

The  approach  taken  used  specially-designed  seven-hole  pressure  probes  to  collect 
pressure  measurements  in  the  model  wake.  One  characteristic  of  canard  wakes  is  a  region 
of  lower  fluid-mechanical  energy  seen  as  a  drop  in  total  pressure.  This  wake  exists 
at  each  point  in  the  airstream  where  the  measured  total  pressure  differs  significantly 
from  the  free-stream  value.  This  region  is  located  by  placing  the  pressure-measuring 
probe  at  many  points  downstream  of  the  model  using  a  three-dimensional  positioning  mech¬ 
anism  known  as  a  traverse.  The  measured  pressures  were  plotted  as  contour  graphs,  where 
the  steep  pressure  gradients  (like  those  occurring  at  the  edge  of  the  canard  wake)  were 
indicated  by  a  collection  of  closely-spaced  contour  lines.  A  series  of  these  graphs 
compose  a  mapping  of  the  flow  field. 

This  entire  procedure  was  reduced  to  four  computer  programs  run  in  sequence.  Each 
program  prompts  the  computer  operator  for  all  necessary  information.  The  time  and  expense 
required  to  perform  the  procedure  is  primarily  a  function  of  the  three-dimensional  tra¬ 
verse  mechanism  and  how  quickly  it  can  move  the  pressure-measuring  probe.  The  reduction 
and  display  of  the  data  takes  less  than  two  additional  minutes  after  data  collection 
is  complete.  The  procedure  has  been  successfully  performed  by  many  Academy  cadets  and 

*Captain,  USAF,  Department  of  Aeronautics,  DFAN 

**2nd  Lieutenant,  USAF,  Department  of  Aeronautics,  DFAN 
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officers  on  a  variety  of  aerodynamic  models.  On  the  basis  of  these  results  the  research 
can  be  considered  successful  in  developing  a  new  technique  which  can  simply,  rapidly, 
and  inexpensively  measure  and  describe  the  wakes  of  wind  tunnel  aerodynamic  models. 

This  paper  consists  of  four  parts:  (1)  a  theory  section  describing  the  expected 
results  of  mapping  a  theoretical  vortex  in  the  free-stream  direction,  (2)  an  apparatus 
description,  (3)  a  procedure  description,  and  (4)  a  discussion  of  the  results  of  mapping 
three  canard  wakes.  An  appendix  consisting  of  the  user  manual  for  the  four  computer 
programs  is  included. 


II.  Theory 

One  of  the  simplest  pieces  of  instrumentation  for  studying  steady  flow  field  proper¬ 
ties  is  the  pressure  probe.  It  is  easy  to  understand  and  manufacture,  its  use  is  wide¬ 
spread,  and  it  directly  measures  one  of  the  most  important  aerodynamic  properties.  The 
seven-hole  pressure  probe  extends  these  characteristics  to  include  direct  measurement 
of  local  flow  angles  and  static,  total  and  dynamic  pressure,  and  allows  calculation  of 
local  velocities  at  angles  up  to  80  degrees  off  axis  without  iteration  (Ref.  1).  These 
properties  make  the  seven-hole  probe  the  best  instrumentation  for  this  research. 

Total  pressure  is  a  measure  of  total  fluid-mechanical  energy  and  results  from  slow¬ 
ing  the  fluid  velocity  to  zero  without  loss.  Static  pressure  is  the  pressure  exerted 
on  an  aerodynamic  surface  parallel  to  the  free  stream.  Dynamic  pressure  is  total  pres¬ 
sure  minus  static  pressure.  All  three  pressures  can  be  used  to  locate  the  canard  wake 
by  observing  where  large  changes  or  gradients  exist  in  the  flow.  Since  these  pressures 
are  affected  by  slight  changes  in  wind  tunnel  velocity  or  room  temperature,  pressure 
coefficients  will  be  used  instead,  where  the  differences  in  local  pressures  and  reference 
pressures  are  divided  by  tunnel  dynamic  pressure,  resulting  in: 


where 


Static  Pressure 
Dynamic  Pressure 


Static 


Total  Pressure 
Dynamic  Pressure 


P  x  -  P 
oL  o 


'TOTAL 


PotL  =  Local  Static  Pressure  (measured  by  probe  flow) 

P^  =  Free  Stream  Static  Pressure 

P  ^  =  Local  Total  Pressure  (measured  by  probe  in  flow) 

PQ  =  Free  Stream  Total  Pressure 

To  anticipate  what  our  coefficients  will  resemble  in  the  canard  wake,  we  will  start 
with  a  simpler  example  of  an  ideal  two-dimensional  vortex  (Ref.  2).  In  an  ideal  two- 
dimensional  vortex,  angular  velocity  increases  exponentially  as  the  distance  to  the  vor¬ 
tex  filament  (radius)  decreases  (Figure  1).  In  reality  the  viscosity  of  the  fluid  reduces 
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Figure  1.  Fluid  Velocity  Versus  Distance  from 
Vortex  Center  for  an  Ideal  Vortex 

the  angular  velocity  long  before  it  reaches  infinity  (Ref.  3).  In  this  region  where 
viscous  effects  are  substantial,  the  angular  velocity  drops  until  it  is  zero  at  the  vor¬ 
tex  center.  Figure  2  shows  the  velocity  distribution  of  a  vortex  with  a  viscous  core 


Figure  2.  Fluid  Velocity  Versus  Distance  from 
Vortex  Center  for  a  Real  Fluid 


(Ref.  3).  At  distances  greater  than  rt,  the  flow  is  approximately  inviscid. 
The  conventional  aerodynamic  pressure  coefficient  is  defined  as 

PA  -  P 

„  _  A  -  . 


pa  %  P„V„2 

where  PA  is  local  static  pressure  at  any  Point  A.  In  a  steady,  incompressible  flow  we 
have 

h  p  V  2  =  P„  -  P 


and  then 


“STATIC 


P  T  -  P 
«L  “ 


P  -  P 
o  « 


P~L  ~ 

*sp  v  2 


which  yields 


cstatic  cp  . 


Cp  is  also  defined  by 


Cp  =  1 


Since  the  conditions  at  infinity  approach  free-stream  conditions,  Cp  and  q  approach 

zero.  As  the  radius  decreases  from  infinity  to  rt,  the  angular  velocity  increases.  This 
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increase  in  velocity  causes  the  static  pressure  to  drop.  Hence,  as  we  go  from  rB  to 

rt>  ^STATIC  ®oes  from  zero  to  a  negative  value. 

Out  to  rt,  viscous  effects  preclude  the  use  of  the  previous  equations.  This  vis¬ 
cous  region  has  recently  been  examined  experimentally,  and  Figure  3  shows  the  typical 
'"STATIC  variat:i-on  in  this  region  as  well  as  the  CS^AT1C  variation  in  th^  inviscid  region 
which  was  previously  discussed.  Thus  the  effect  of  vorticity  shed  from  the  canard  of 


r 

Figure  3.  Typical  Behavior  with  Distance 

from  the  Vortex  Center 


our  wind  tunnel  model  induces  higher  velocities  and  therefore  causes  a  decrease  in  local 
static  pressure  at  points  in  the  flow  near  that  vorticity.  Unfortunately,  the  CgTATjC 
coefficient  does  not  allow  us  to  clearly  define  the  edge  of  the  wake.  This  is  evident 
in  Figure  3  where  rt  has  been  arbitrarily  located  near  the  inflection  point  of  the  curve. 

The  CT0TAI  coefficient,  however,  can  be  used  to  locate  the  edge  of  the  wake.  Ber¬ 
noulli's  equation  shows  that  total  pressure  is  constant  in  isentropic,  incompressible 
flow.  Using  our  definition. 


o 


it  follows  that  the  coefficient  of  total  pressure  in  the  inviscid  region  is  zero.  Since 
viscosity  causes  a  loss  in  total  mechanical  energy,  PQ  decreases  at  radii  less  than  rt. 
As  P  drops,  C  will  decrease  to  a  negative  value.  Figure  4  approximately  shows 

OLi  1U 1AL 

this  behavior.  We  thus  expect  the  canard  wake  to  appear  as  a  mound  with  a  depression  in 


Figure  4.  cjoTAL  (Theoretical)  Versus 
Distance  from  the  Vortex  Center 


it  on  the  contour  maps.  If  we  define  the  edge  of  the  wake  as  points  where  the 
total  pressure  begins  to  deviate  from  free  stream,  i.e.,  rt,  we  have  the  indicator  of 
the  wake  location  we  desire. 
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An  aerodynamic  problem  associated  with  the  canard  wake  is  the  lack  of  dynamic  pres¬ 
sure  resulting  from  viscous  losses.  If  there  is  a  connection  between  the  local  dynamic 
pressure  and  our  definition  of  cSfATIQ»  might  also  serve  to  locate  the  edge  of  the 
wake.  Our  definition  is 

CDYN  =  CTOTAL  ~  CSTATIC 


Manipulating  these  equations  gives 


CDYN  CTOTAL  ~  CSTATIC 


P  -  P  is  the  local  dynamic  pressure,  while  P  -  P^  is  the  free-stream  dynamic  pres- 

O  L  O 

sure.  Thus,  is  proportional  to  the  ratio  of  local  dynamic  pressure  to  free-stream 

dynamic  pressure. 

Since  the  conditions  at  infinity  approach  those  of  free  stream,  0^^  will  be  approx¬ 
imately  zero.  CDYN  *  ICstaticI  fr°m  infinit>r  1:0  rt  since  SoTAL  =  °'  At  the  center 
of  the  vortex,  local  velocity  and  thus  local  dynamic  pressure  is  zero,  causing  CDYN  to 
equal  -1.  From  rt  to  the  center,  goes  from  lcSTATICl  to  as  in  Figure  5.  Thus 


Figure  5.  ^dynamic:  CTheoretical)  Versus 
Distance  from  the  Vortex  Center 

0^^  also  provides  an  indication  of  wake  location  when  it  reaches  the  point  of  its  local 
maximum. 

To  this  point  we  have  considered  only  a  single,  isolated  two-dimensional  vortex. 

The  actual  canard  wake  is  much  more  complex  and  includes  many  other  factors  such  as  flow 
separation,  multiple  vortices,  and  fuselage  interference.  We  therefore  must  recognize 
that  this  theory  only  provides  us  an  approximate  guide  to  explain  this  complex  wake  be¬ 
havior. 
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To  measure  this  complex  flow  we  used  a  seven-hole  probe.  Then,  using  the  procedure 
outlined  in  Ref.  1,  we  calculated  the  angle  of  attack,  a,  sideslip  angle,  6,  and  two 
pressure  coefficients  CQ  and  directly  from  the  seven  measured  pressures.  C'i'OTAL  and 
CSTAT1C  are  calculated  using  Co  and  ,  and  CDyN  is  calculated  from  and  CgTATIC. 

After  the  data  is  collected  and  all  coefficients  are  calculated,  we  graphically 
display  a  three-dimensional  surface  consisting  of  the  Y  and  Z  location  of  each  point 
and  the  desired  coefficient.  The  graphs  are  displayed  with  either  an  axonometric  or 
contour  projection.  (The  1980  Spring/Summer  Edition  of  the  Aeronautics  Digest  displays 
an  axonometric  projection  on  its  cover.) 

III.  Apparatus 

The  apparatus  we  used  to  perform  these  tests  consisted  of  a  model,  seven-hole  probe, 
traverse  mechanism,  wind  tunnel,  and  data  acquisition  system. 

A.  Model 

The  model  used  in  the  wind  tunnel  was  a  variable-height  canard  with  swept-back 
wings.  The  canard  was  positioned  at  either  a  high,  middle,  or  low  position.  Because 
of  the  size  of  the  subsonic  wind  tunnel  used,  we  utilized  a  half-span  model.  This  half¬ 
span  model  was  a  replica  of  one  used  at  NASA  in  their  Langley  Wind  Tunnel  Facility.  The 
model  differs  slightly  in  several  linear  dimensions  from  NASA's  model  because  of  construc¬ 
tion  problems.  The  dimensions  of  the  model  with  the  canard  in  the  low  position  are  shown 
in  Figures  6  and  7. 


k- - 4 

0  " 


Figure  6.  Canard  Model  (top  view) 
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57.0  " 


Figure  7.  Canard  Model  (side  view) 


B.  Probe 

We  used  the  seven-hole  probe  built  by  Hollenbaugh  and  reported  by  Gallington 
(Ref.  L).  This  probe  was  calibrated  last  year  in  the  subsonic  wind  tunnel.  The  hole 


C.  Traverse  Mechanism 

We  used  a  previously-designed  three-dimensional  traverse  mechanism  and  control 
unit  developed  by  the  Frank  J.  Seiler  Research  Laboratory  for  the  subsonic  wind  tunnel. 
The  axes  of  the  traverse  are  shown  in  Figure  9.  The  traverse  is  computer-driven  in  the 
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Y  and  Z  direction  and  manually  controlled  in  the  X  direction. 

D.  Wind  Tunnel 

The  experiment  used  the  large  closed-circuit  wind  tunnel  which  has  a  two  foot 
by  three  foot  test  section  and  a  velocity  range  from  50  ft/sec  to  400  ft/sec  at  atmos¬ 
pheric  pressure.  The  tunnel  is  driven  by  a  200-horsepower  synchronous  motor  coupled 
with  a  fluid  drive  unit  which  is  used  to  control  the  RPM  of  the  four-bladed  compressor. 

The  tunnel  does  not  have  a  heat  exchanger,  which  caused  tunnel  velocity  and  temperature 
to  increase  throughout  each  run. 

E.  Data  Acquisition  System 

Seven  Statham  PM6TC  ±1-350  transducers  were  used  to  convert  the  probe  pressures 
to  voltages.  Ten  Ectron  Model  560  amplifiers  were  used.  A  PDP-11/45  computer  was  used 
to  drive  the  traverse  and  collect  and  reduce  the  data.  The  following  peripherals  and 
software  were  used  in  the  data  acquisition/reduction  process: 

1)  A  teletypewriter  used  as  the  system  list  device 
. 2)  A  control/graphic  display  terminal  used  to  display  plots 

3)  A  fast  access  40-megabyte  disk  mass  storage  device 

4)  A  system  which  converts  analog  signals  from  the  transducers  and  traverse  into 
digital  form  suitable  for  computer  use 

5)  Standard  graphics  software  for  the  display  terminal 

6)  A  high  quality,  flat  bed  pen  plotter 

7)  A  thermal  hard  copy  device  hooked  into  the  display  terminal 
IV,  Procedure 

The  actual  test  procedure  can  be  considered  the  sequential  execution  of  four  compu¬ 
ter  programs.  They  were  VOLCAL.FTN,  TOPWNG.FTN,  C0EF7.FTN,  and  SCALAR. FTN.  See  Appen¬ 
dix  A  for  user  instructions  for  each  program.  The  apparatus  had  to  be  installed  in  the 
wind  tunnel  to  use  the  first  two  programs,  but  the  latter  two  required  only  the  PDP- 
11  /4  5  computer. 

The  first  step  after  the  model,  traverse,  and  probe  were  installed  in  the  wind  tun¬ 
nel  was  to  run  VOLCAL.FTN  on  the  computer.  VOLCAL  used  a  tri-linear  scheme  to  define 
the  measurement  region  as  a  variable  number  of  points  on  a  variable  number  of  planes 
in  a  linear  prismatoid.  VOLCAL  stored  the  three-dimensional  coordinates  of  all  the  points 
in  each  plane  in  a  data  file,  VOLUME.DAT. 

The  next  program,  TOPWNG.FTN,  drove  the  probe  to  each  point  defined  in  VOLUME.DAT. 

It  measured  the  pressures  on  each  of  the  probe  holes  along  with  the  tunnel  total  pressure, 
static  pressure,  and  temperature.  After  completely  measuring  all  the  points  in  a  given 
plane,  the  program  wrote  the  data  out  to  a  file  called  WNGMAP.DAT. 

C0EF7.FTN  reduced  all  the  data  stored  in  WNGMAP.DAT.  Seven  sets  of  calibration 
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coefficients  had  been  previously  determined  based  on  the  highest  of  the  seven  pressures. 
Two  angles  and  two  pressure  coefficients  were  determined  from  the  measured  pressures 
and  calibration  coefficients  which  described  the  local  total  and  static  pressures  and 
the  local  flow  angularity  (Ref.  1).  The  reduced  data  was  stored  in  a  file  called 
CONTUR.DAT. 

The  last  program,  SCALAR. FTN,  produced  the  graphs.  The  graphs  were  projected  on 
the  display  terminal  and  hard  copies  were  made  on  the  pen  plotter.  Both  axonometric 
and  contour  plots  were  created. 

V.  Discussion 

Data  were  collected  on  three  different  canard  wakes  using  the  parametric  model  de¬ 
scribed  above.  Canard  positions  of  0  inches,  1.125  inches,  and  2.75  inches  above  the 
plane  of  the  aft-swept  wing  were  investigated.  All  were  at  11  degrees  angle  of  attack. 

In  addition,  oil  flow  visualizations  were  performed  on  the  high  canard  model  at  angles  of 
attack  of  11  and  22  degrees.  Eight  planes  of  data  perpendicular  to  the  free  stream  were 
collected  for  each  wake.  The  oil  flows  will  be  discussed  first,  followed  by  a  discussion 
of  the  wakes. 

A  sketch  of  the  oil  flow  at  11  degrees  angle  of  attack  is  shown  in  Figure  10  and 
it  shows  the  main  features  of  the  flow  over  the  canard  and  wing.  On  the  canard  we  see 
a  separation  line  and  a  reattachment  line.  Between  the  two  the  oil  shows  the  existence 
of  a  large  degree  of  spanwise  flow  caused  by  a  vortex.  This  vortex  will  be  called  the 
canard  leading  edge  separation  vortex.  It  exists  at  the  low  angle  of  attack  of  11  degrees 
because  the  leadiiig  edge  radius  of  the  canard  is  near  zero.  Outside  the  separation  line 


Figure  10.  Oil  Flow  at  Alpha  =11  Degrees 


the  oil  has  no  streaks,  indicating  the  flow  is  severely  separated.  To  describe  this 
region  by  potential  flow  would  be  very  difficult. 

On  the  wing  a  similar  separated  flow  exists.  One  different  feature  is  the  occurrence 
of  a  weak  secondary  reattachment  line  outboard  of  the  separation  line.  The  small  second- 
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ary  separation  vortex  which  causes  this  phenomenon  is  well-documented  for  highly  swept 
wings  (Ref.  4). 

The  oil  flow  inside  the  reattachment  line  shows  the  existence  of  a  second,  very 
weak  vortex  emanating  from  the  non-faired  canard-fuselage  junction.  The  oil  flow  at 
22  degrees  shows  the  same  results,  differing  only  in  the  amount  of  separated  flow.  Since 
one  of  the  reasons  for  mapping  these  wakes  was  to  aid  in  their  correct  potential  flow 
modeling,  we  chose  to  concentrate  our  efforts  at  11  degrees  angle  of  attack. 

The  middle-positioned  canard  wake  will  be  discussed  in  detail  first.  Comparisons 
with  the  low-  and  high-positioned  canard  wakes  will  follow. 

Contour  maps  showing  lines  of  constant  local  total  pressure  in  eight  planes  of  data 
taken  perpendicular  to  the  free  stream  are  presented  in  Figure  11.  The  X-axis  is  parallel 
to  the  free  stream  and  therefore  is  at  an  11-degree  angle  with  the  longitudinal  axis 
of  the  fuselage.  The  Y-axis  of  the  canard  is  increasing  in  the  inboard  direction  and 
its  origin  is  at  the  trailing  edge  of  the  canard  tip. 

The  first  plane  of  data  at  X  =  .3  inches  behind  the  canard  tip  shows  a  relatively 
concentrated  region  of  low  total  pressure  inboard  and  above  the  canard  tip  (Figure  11a). 
The  point  of  minimum  pressure  marks  the  center  of  the  canard  leading  edge  separation 
vortex.  Outboard  of  this  is  another  low-pressure  region  located  above  the  canard  tip, 
marking  the  center  of  the  canard  tip  vortex.  These  two  vortices  in  such  close  proximity 
create  a  saddle-like  structure  when  plotted  in  three  dimensions  (total  pressure  plotted 
orthogonally  to  the  Y  and  Z  spacial  dimensions.  Both  vortices  are  rotating  in  a  clock¬ 
wise  sense,  when  looking  into  the  free  stream.  The  low-pressure  regions  located  seven 
inches  inboard  are  the  result  of  fuselage-canard  interactions  because  of  the  lack  of 
well-designed  fairing. 

As  we  move  aft  in  Figures  lib  through  1 lh  we  see  the  local  total  pressure  and 
coefficient  increasing  as  viscous  effects  reduce  the  magnitude  of  velocity  gradients 
in  the  flow.  The  integrated  total  pressure  deficit,  of  course,  remains  the  same. 

The  relative  strengths  of  the  canard  tip  and  separation  vortices  are  evidenced  by 
the  magnitudes  of  the  pressure  coefficient  and  also  by  the  relative  movements 

of  the  vortex  centers.  The  movement  of  the  canard  tip  vortex  in  a  clockwise  fashion 
about  a  point  close  to  the  separation  vortex  center  shows  the  separation  vortex  is  much 
stronger,  inducing  more  movement  in  the  tip  vortex  than  it  can  induce  in  the  separation 
vortex. 

In  Figures  lie  through  l lh  we  see  a  new  structure  forming  at  the  bottom  and  moving 
outboard  as  we  go  aft.  This  is  the  wing  leading  edge  separation  vortex.  The  entire 
structure  could  not  be  mapped  due  to  geometric  limitations  on  the  traverse  mechanism. 
Figures  ltg  and  1 1 h  show  the  wing  separation  vortex  is  assymmetrical  with  a  small  flow 
structure  located  inboard.  We  do  not  know  the  cause  of  this  small  structure  which  con¬ 
sistently  appears.  We  suggest  that  the  leading  edge  separation  vortex  on  a  highly  swept 
wing  needs  to  be  studied  in  greater  detail. 
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Nowhere  in  Figures  11a  through  llg  do  we  see  the  effects  of  the  trailing  vorticity 
shed  from  the  canard  trailing  edge  as  a  result  of  changes  in  the  spanwise  load  distribu¬ 
tion.  It  may  be  that  the  expected  feature  has  a  characteristic  dimension  that  is  much 
smaller  than  the  resolution  limit  of  the  probe.  This  also  needs  to  be  studied  in  more 
detail. 

The  approximate  locations  of  the  centers  of  the  canard  tip  and  separation  vortices 
versus  X  are  given  in  Table  1.  As  we  go  aft  the  canard  tip  vortex  goes  up  and  inboard. 


Table  1 

CANARD  VORTEX  CENTER  LOCATION 


MIDDLE  CANARD  ALPHA  =11° 

TIP 

SEPARATION 

X 

Y 

Z 

Y 

Z 

.3 

.5 

.7 

2.7 

1.2 

3.1 

.4 

1.0 

2.7 

1.0 

5.9 

.2 

1.3 

2.7 

.9 

8.7 

.2 

1.7 

2.8 

.9 

11.5 

.3 

2.1 

2.9 

.6 

14.3 

- 

- 

3.1 

.4 

17.2 

- 

3.  1 

0.0 

20.0 

- 

- 

3.  1 

-.5 

while  the  separation  vortex  goes  down.  This  downward  movement  is  approximately  half 
the  angle  of  attack. 

The  local  static  pressure  and  static  pressure  coefficient  are  inversely  related 
to  the  magnitude  of  local  velocity.  Examining  Figure  12a,  we  see  a  low-pressure  and 
thus  high-velocity  flow  structure  at  the  Y  and  Z  coordinates  (3.5,  1.9).  Just  below 
and  inboard  exists  a  relatively  high-pressure,  low-velocity  structure  forming  a  "static 
pressure  doublet."  One  possible  explanation  for  this  structure  is  that  the  low-velocity 
region  represents  the  vortex  core  with  its  region  of  rigid  body  rotation.  If  the  high 
velocity  represents  the  edge  of  the  viscous-potential  flow  boundary,  by  theory  it  should 
encircle  the  vortex  core.  The  fact  that  it  does  not  suggests  the  answer  may  be  that  the 
characteristic  of  the  flow  structure  is  less  than  the  resolution  limit  of  the  probe.  This 
may  also  be  causing  the  discrepancy  between  the  center  of  the  low-velocity  region  and 
the  center  of  the  total  pressure  loss.  The  resolution  limit  needs  to  be  quantified. 

An  interesting  feature  of  the  static  pressure  map  is  the  linear,  high-pressure/ 
low-velocity  ridge  located  between  points  (0,  .3)  and  (3.5,  .3).  This  appears  to  be 
the  wake  of  the  shed  vorticity  from  the  canard  trailing  edge  because  of  changes  in  the 
spanwise  load  distribution.  If  we  examine  static  pressures  on  a  line  perpendicular  to 
this  structure  we  see  a  low-pressure/high-velocity  region,  a  high-pressure/ low-velocity 
wake,  and  a  low-pressure/high-velocity  region  below.  This  classic  wake  structure  is 
quickly  dissipated  as  we  go  aft  because  of  viscous  effects. 
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Figure  13.  Dynamic  Pressure  Coefficient 
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Examining  the  rest  of  Figure  12,  we  see  that  the  major  features  of  the  flow  are 
the  canard  leading  edge  separation  vortex  in  the  center  and  the  wing  separation  vortex 
below  it.  A  third  feature  that  deserves  comment  is  the  shape  of  the  static  pressure 
contours  outboard  of  the  vortices,  which  is  essentially  a  series  of  parallel  lines.  This 
means  the  flow  velocity  inboard  is  higher  than  outboard  in  a  regular  fashion.  This  is 
due  to  the  induced  flow  field  of  the  60-degree  swept  wing.  The  inboard  portion  of  each 
figure  is  above  the  wing  surface  while  the  outboard  portion  is  in  front  of  the  leading 
edge.  We  expect  the  velocity  above  the  wing  to  be  higher  than  the  induced  velocity  in 
front  of  the  leading  edge. 

The  local  dynamic  pressure  is  presented  in  Figure  13  as  contour  plots  of  the  dynamic 
pressure  coefficient.  Since  this  coefficient  is  simply  the  difference  of  the  total  and 
static  pressure  coefficients,  it  shows  the  features  of  both.  Figure  13a  illustrates 
the  effect  of  the  linear  canard  wake  ridge  observed  in  Figure  12a.  The  high  and  low 
"static  pressure  doublet"  causes  only  very  slight  modification  to  the  upper  right  side 
of  the  total  pressure  structure  of  Figure  11a.  Figures  13b  through  13h  behave  as  ex¬ 
pected. 

Two  other  canard  wakes  were  mapped  with  the  canard  at  0  inches  and  2.75  inches  above 
the  plane  of  the  wing.  Since  these  wakes  gave  very  similar  results  to  the  middle  canard, 
we  will  summarize  the  total  pressure  data  by  giving  the  coordinates  of  the  canard  lead¬ 
ing  edge  separation  vortex  and  tip  vortex  in  Tables  2  and  3.  The  origin  of  the  coordi¬ 
nate  system  for  all  three  tables  is  the  same. 
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Table  3 

CANARD  VORTEX  CENTER  LOCATION 


HIGH  CANARD  ALPHA  =  11 

O 

TIP 

SEPARATION 

X 

Y 

z 

Y 

z 

.2 

-1.0 

2.3 

2.1 

2.6 

3. 1 

-  .4 

2.6 

2.5 

2.4 

5.9 

0.0 

2.8 

2.8 

2.1 

8.8 

- 

- 

3.0 

1.8 

11.7 

- 

- 

3.0 

1.6 

14.6 

- 

- 

3.0 

1.6 

17.4 

- 

- 

3.0 

1.3 

20.3 

- 

- 

3.2 

1.0 

VI .  Conclusions 

We  have  documented  the  development  and  use  of  a  new  technique  of  conducting  flow 
field  surveys  using  a  seven-hole  probe.  These  surveys  allow  the  calculation  of  most 
of  the  important  flow  field  parameters,  including  local  static  and  total  pressures.  This 
technique  was  used  to  survey  the  canard  wake  of  three  separate  models  at  subsonic  incom¬ 
pressible  speeds.  Our  results  allow  us  to  locate  the  canard  wakes  passing  above  the 
wing  surface.  They  also  indicate  a  possible  discrepancy  of  location  of  wakes  depending 
on  the  type  of  data  examined,  that  is,  static,  total,  or  dynamic  pressure.  The  relative 
accuracy  of  the  method  is  undetermined  as  is  the  resolution  limit  of  the  probe  and  the 
frequency  response  of  the  apparatus.  The  tests  were  run  only  at  incompressible  flow 
velocities.  Despite  these  criticisms,  a  valuable  educational  and  research  technique 
for  quantitative  flow  field  visualization  has  been  developed. 
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Appendix  A 

This  appendix  contains  the  user  instructions  for  VOLCAL.FTN,  TOPWNG.FTN,  C0EF7.FTN, 
and  SCALAR. FTN.  All  the  left- justified  capitalized  lines  are  either  responses  to  or 
replies  from  the  PDP-11/45  computer. 

VOLCAL.FTN  USER'S  MANUAL 

ENTER  A  HEADING  FOR  THE  GRAPHS  GENERATED  BY  SCALAR. FTN.  DESCRIBE  THE  MODEL  CONFIGURATION, 
THE  ANGLE  OF  ATTACK,  AND  ANY  OTHER  DESIRED  INFORMATION.  MAX  LENGTH  IS  40  CHARACTERS. 

The  user  now  enters  a  heading  such  as 

ALPHA  =  11  MIDDLE  CANARD,  SWEPT  BACK  WING 

The  computer  then  prints  out  what  is  entered  and  asks  if  it  is  correct.  If  correct, 

Y  is  entered.  If  incorrect,  N  is  entered  and  the  computer  again  requests  the  title 
information. 

ENTER  XA,  YA,  ZA 

The  computer  asks  for  the  X,  Y,  and  Z  coordinates  for  corner  point  A.  Figure  A-l 
shows  the  location  of  points  A  through  H.  The  dimensions  of  each  coordinate  are 


Figure  A-l.  Measurement  Volume 

in  inches  and  are  referenced  to  the  (0,0,0)  point  at  the  outboard  tip  of  the  middle 
canard.  Point  A  is  the  lower  right  corner  of  the  volume  as  referenced  to  the  front 
of  the  model. 

ENTER  XB,  YB,  ZB 


ENTER  XH,  YH,  ZH 

ARE  THESE  CORRECT?  (Y  OR  N) 

If  N  is  entered,  the  computer  asks  for  the  coordinates  again. 

ENTER  THE  DESIRED  NUMBER  OF  PLANES  OF  DATA  IN  THE  X  DIRECTION,  THE  NUMBER  OF  POINTS 
DESIRED  IN  THE  Y  DIRECTION  (OUT  THE  CANARD  TIP),  AND  THE  DESIRED  NUMBER  OF  POINTS  IN 
THE  Z  DIRECTION,  RESPECTIVELY.  USE  REAL  FORMAT. 

We  used  four  planes  of  data  in  the  X  direction,  20  points  in  the  Y  direction,  and 
14  points  in  the  Z  direction.  The  input  was  thus: 
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4,20,  14 

The  computer  then  prints  out  the  values  entered  above  and  asks  if  they  are  correct. 
If  N  is  entered,  then  it  will  go  back  and  ask  the  above  question  again.  The  program 
now  executes  to  completion. 


TOPWNG.FTN  USER'S  MANUAL 

ENTER  VOLUME , DAT  VERSION  NUMBER 

The  user  now  enters  the  version  number  of  VOLUME.DAT  that  corresponds  to  the  model 
configuration  being  tested. 

ENTER  TEMP,  PATM  (USE  REAL  FORMAT) 

The  user  enters  the  room  temperature  (°F)  and  the  room  pressure  (in.  Hg) .  After 
the  temperature  and  pressure  are  entered,  the  computer  prints  them  out  and  asks 
if  they  are  correct. 

START  TUNNEL  AND  ENTER  GO  WHEN  READY 

The  user  starts  the  tunnel  and  enters  GO  when  the  air  is  stabilized  at  the  desired 
velocity. 

ENTER  WHAT  VELOCITY  THE  RUN  IS  AT  AND  ANY  OTHER  INFORMATION  UP  TO  A  TOTAL  OF  20  CHARAC¬ 
TERS  LONG  YOU  WISH  TO  HAVE  IN  THE  GRAPH  TITLE 

The  user  enters  something  similar  to 

VEL  =  +100  FPS 

The  computer  prints  out  the  above  and  checks  if  it  is  correct.  If  incorrect, 

the  user  enters  N  and  the  computer  asks  for  it  again. 

THE  NEXT  PLANE  THAT  IS  RUN  IS  AT  X  =  .3.  DO  YOU  WANT  TO  RUN  THIS  PLANE  OF  DATA?  IF 
YOU  DO,  MOVE  THE  PROBE  TO  THE  DISPLAYED  X  VALUE  AND  ENTER  GO.  IF  YOU  WANT  ANOTHER  PLANE, 
THEN  ENTER  NO. 

The  user  acts  accordingly.  As  soon  as  GO  is  entered,  the  program  will  start  taking 

data.  When  a  plane  is  completed,  the  following  message  is  printed  out  on  the  screen 

if  there  is  another  plane  of  data  to  be  taken. 

MOVE  TRAVERSE  MECHANISM  TO  NEXT  X  VALUE,  X  =  3.1. 

ENTER  GO  WHEN  READY  TO  CONTINUE  RUN. 

As  soon  as  GO  is  entered,  the  program  takes  another  plane  of  data.  The  computer 
carriage  returns  when  it  completes  the  last  plane  of  data.  t 

C0EF7.FTN  USER'S  MANUAL 

ENTER  WNGMAP.DAT  VERSION  NUMBER 

The  user  enters  the  version  number  of  the  WNGMAP.DAT  data  file  to  be  reduced. 

SCALAR. FTN  USER'S  MANUAL 

DO  YOU  WANT  AN  HP-4662  HARDCOPY? 

If  the  user  wants  a  hard  copy,  he  sets  up  the  4662  printer  with  the  paper,  defines 
the  borders  (if  needed),  and  enters  Y.  N  is  entered  if  the  user  wants  the  graph 
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displayed  on  the  screen. 

ENTER  CONTUR.DAT  VERSION  NUMBER 

The  user  enters  the  corresponding  version  number  of  the  plane  of  data  to  be  dis¬ 
played  . 

ENTER  VALUE  FOR  2  COORDINATE  DATA 

1  CA 

2  CB 


13  CDYN 

The  user  enters  the  number  of  the  data  to  be  displayed. 

ENTER  A1,A2,A3,S1,S2,S3 

The  user  enters  0,  120,  90,  .4,  .4,  .4,  1  if  an  axonometric  projection  is  desired 
or  0,  90,  0,  .4,  .4,  .4  if  a  contour  projection  is  desired.  The  three  .4's  are 
the  scaling  factors.  Changing  these  will  change  the  size  of  the  graph. 

If  a  contour  projection  is  selected,  the  computer  will  print  out: 

ZMIN  ZMAX  INC 

-1.235  .243  .215 

DO  YOU  WISH  TO  CHANGE  LIMITS? 

If  more  or  less  intervals  are  desired,  the  user  enters  Y.  If  the  given  intervals 
are  satisfactory,  N  is  entered.  If  Y  is  entered,  the  computer  will  print  out: 
ENTER  ZLOW,  ZMAX,  ZINC 

The  user  now  enters  the  required  information  such  as 
-1.2,  .2,  .3 

The  program  now  plots  the  graph  either  or,  the  terminal  screen  or  on  the  4662  flat 
bed  pen  plotter. 
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EXAMINING  A  RULE  OF  THUMB  FOR  THE  RELATION  BETWEEN 
CAMBER  AND  ZERO-LIFT  ANGLE  OF  ATTACK 

E..S.  Jumper* 

Abstract 

This  paper  examines  the  rule  of  thumb  that  states  that  an  airfoil's  zero-lift  angle 
of  attack  is  roughly  equal  in  degrees  to  minus  the  magnitude  of  the  percent  camber. 

The  examination  consists  of  looking  at  coef f icient-of-lif t  versus  angle-of-attack  curves 
for  selected  NACA  airfoils  and  comparing  this  experimentally-determined  zero-lift  angle- 
i  f-attack  data  with  that  predicted  by  the  rule  of  thumb.  Thin-airfoil  theory  is  also 
employed  to  examine  the  validity  of  the  rule.  Finally,  the  paper  presents  an  example 
of  how  the  rule  can  be  used  to  help  interpret  data  by  relating  the  rule's  use  in  inter¬ 
preting  wind  tunnel  data  obtained  for  a  K1 ine-Fogleman  airfoil. 

1.  Introduction 

In  engineering  disciplines  there  are  certain  relationships  that  hold  in  such  a 
general  way  that  they  are  referred  to  as  rules  of  thumb.  As  Webster  would  define  it, 
a  rule  of  thumb  is  a  general  principle  regarded  as  roughly  correct  and  helpful  but  not 
intended  to  be  scientifically  accurate.  In  this  paper  I  am  going  to  examine  a  rule 
of  thumb  I  found  in  an  introductory  aerodynamics  text  by  Clancy  (Ref.  1).  This  rule 
of  -thumb  states  that  an  airfoil's  zero-lift  angle  of  attack  is  roughly  equal  in  degrees 
to  minus  the  magnitude  of  the  percent  camber. 

Figure  l  shows  a  typical  coef f ic ient-of-lif t ,  Cf  ,  versus  angle-of-attack,  t ,  curve 
for  an  airfoil  with  positive  camber.  Notice  that  the  airfoil  must  \-r  olaced  at  a  nega¬ 
tive  a  before  C /  is  equal  to  zero.  The  angle  at  which  rv  ',^als  zeio  i-  called  the 
zero-lift  angle  of  attack,  Referring  to  Figure  2,  camber  is  the  maximum  distance 

from  the  mean  camber  line  to  the  chord  line.  This  distance  is  indicated  by  z  in  Figure 

2.  The  camber  is  usually  specified  as  a  percent  of  the  length  of  the  chord  line,  c. 
Thus,  a  camber  of  four  percent,  for  example,  would  mean  that  zc  was  equal  to  0.04c  and 
would  be  considered  positive  (concave  up)  unless  otherwise  stated.  Using  our  rule  of 
thumb  and  the  example  of  the  four-percent  camber  airfoil,  we  would  expect  to  find 

for  this  airfoil  roughly  equal  to  -4  degrees. 

A  cursory  examination  of  Cj  /a  data  from  Abbott  and  Doenhoff  (Ref.  2)  for  various 
airfoils  seems  to  verify  the  rule,  but  how  far  can  we  push  it?  This  paper  examines 
the  applicability  of  the  rule  by  taking  a  closer  look  at  the  Abbott  and  Doenhoff  airfoil 
data  and  by  examining  a  theoretical  treatment  of  the  effect  of  camber  on  a  q,  using 
thin-airfoil  theory.  Finally,  the  paper  presents  an  example  of  using  the  rule  as  a 
tool  for  interpreting  experimental  data.  The  example  involves  looking  at  a- so-called 
Kline-Fogleman  airfoil. 

*Major,  USAF,  Associate  Professor  of  Aeronautics,  DFAN 
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ANGLE  OF  ATTACK 


'Angle  at  which  the  coefficient  of  lift 
is  zero  (zero-lift  ongle  of  attack,  aL0) 


Figure  1.  Typical  Coef f icient-of-Lif t  Versus 
Angle-of-Attack  Curve  for  an  Airfoil 
with  Positive  Camber 


MEAN  CAMBER  LINE 


MAXIMUM 

THICKNESS 


CHORD  LINE 


Figure  2.  Airfoil  Nomenclature 
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The  most  accessible  airfoil  data  for  examining  the  extent  of  applicability  of  the 
rule  of  thumb  relating  camber  and  a,,  is  found  in  Ref.  2.  In  this  reference  is  a  large 
sample  of  NACA  airfoil  data  assembled  by  Abbott  and  Doenhoff.  Of  the  set  of  data,  the 
NACA  four-  and  five-digit  series  are  easiest  to  interpret  in  terms  of  determining  the 
camber.  In  the  four-digit  NACA  series,  the  first  digit  gives  the  maximum  camber  of 
the  airfoil  in  percent  of  the  chord  (percent  camber) ,  the  second  number  gives  the  loca¬ 
tion  of  the  maximum  camber  in  tenths  of  chord,  and  the  last  two  numbers  give  the  maxi¬ 
mum  thickness  in  percent  of  chord.  Thus  an  NACA  4415  designated  airfoil  describes  a 
wing  section  that  has  a  maximum  camber  of  0.04c,  located  at  0.4c  measured  from  the  lead¬ 
ing  edge  of  the  airfoil,  and  a  maximum  thickness  of  0.15c,  where  c  is  the  chord  length. 
Similarly,  in  the  NACA  five-digit  series,  the  first  digit  again  indicates  the  percent 
camber,  but  the  next  two  numbers  represent  twice  the  percent  of  chord  of  the  location 
of  maximum  camber,  and  the  last  two  digits  again  indicate  the  thickness.  Thus,  an  NACA 
23012  designated  airfoil  describes  a  wing  section  which  has  a  maximum  camber  of  0.02c, 
located  at  0.15c  (.30/2)  from  the  leading  edge  with  a  maximum  thickness  of  0.12c. 

If  we  now  look  at  some  coef f icient-of-lif t  (C[)  versus  angle-of-attack  (a)  curves 
(Cf /a  curves)  gathered  from  wind  tunnel  tests,  we  can  examine  the  applicability  of  the 
rule  of  thumb.  Let  us  first  look  at  the  case  of  a  symmetrical  airfoil,  that  is,  one 
that  has  no  camber.  Figure  3  shows  two  examples  of  symmetrical  airfoils  -  the  NACA 
0006  and  the  NACA  0012  (Ref.  2).  As  both  the  rule  and  common  sense  would  predict,  the 
zero  percent  cambered  airfoil  has  zero  degrees  of  shift  in  the  zero-lift  angle  of  attack. 
Now  look  at  Figures  4,  5,  and  6  which  show  Cj /a  data  for  airfoils  that  have  one,  two, 
and  four  percent  camber  respectively  (Ref.  2).  Note  that  our  rule  of  thumb  correctly 
predicts  over  a  wide  range  of  percent  thickness  (ranging  from  8  percent  to  24  per¬ 
cent).  One  should  also  note,  however,  that  in  all  of  these  cases  the  location  of  the 
maximum  camber  is  at  0.4c. 

If  we  examine  data  from  the  same  source  (Ref.  2)  at  maximum  camber  locations  other 
than  0.4c,  we  begin  to  note  some  disagreement  between  the  ^  predicted  by  the  rule 
of  thumb  and  the  data.  As  the  maximum  camber  location  moves  further  forward  of  the 
0.4c  location,  we  see  that  the  rule  overpredicts  the  magnitude  of  |aLQ|.  For  example, 
Figure  7  shows  the  effect  of  moving  the  maximum  camber  forward  by  0.25c  to  0.15c.  Here 
the  rule  overpredicts  the  magnitude  of  |c*^ by  25  percent.  For  Two-percent  camber, 
the  curve  shifts  to  the  left  only  about  1.5  degrees,  or  about  -0.75  degrees/percent, 
instead  of  -1.0  degrees/percent  predicted  by  the  rule. 

To  examine  the  effect  of  moving  the  maximum  camber  aft  of  the  0.4c  -.ocation,  we 
must  go  to  airfoil  designations  which  do  not  indicate  camber  directly  by  their  number 
designation.  Of  the  NACA-designated  airfoils  reported  by  Abbott  and  Doenhoff,  those 
with  a  camber  line  designation  of  NACA  a  =  0  and  0.1,  maximum  camber  occurs  at  0.3c 
for  a  **  0.4  at  0.4c,  for  a  «  0.5  and  0.6  at  0.45c,  and  for  a  “  0.7,  0.8,  O.d,  and  1.0 
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SECTION  ANGLE  OF  ATTACK,  a0  (degrees) 

Figure  5.  NACA  2408  and  NACA  2424 
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SECTION  ANGLE  OF  ATTACK,  a0 


SECTION  ANGLE  OF  ATTACK,  a0  (degrees) 

Figure  7.  NACA  23012  and  NACA  23024 
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at  0.5c.  In  the  series  which  have  an  "a"  designation,  consisting  of  a  two-digit,  sub¬ 
script-,  three-digit  number,  we  have  an  indirect  indicator  of  the  percent  camber.  Of 

the  last  three  digits,  the  first  indicates  the  design  coefficient  of  lift  (C /  .)  in  tenth 

so  that  an  NACA  65j  -  212  (a  =  0.6)  airfoil  has  a  design  lift  coefficient,  C j  of  0.2 

(the  last  two  numbers  are  the  percent  maximum  thickness).  The  "a"  designator  (i.e., 
a  =  0.6)  refers  to  a  specific  mean  camber  line  distribution  which  is  enumerated  in  Ref. 

2,  based  on  a  design  coefficient  of  1.0.  As  mentioned  before  in  the  a  =  0.6  case,  maxi¬ 

mum  camber  occurs  at  0.45c,  but  further,  the  camber  at  that  point  is  0.0737c  for  C j 
for  a  design  coefficient  of  1.0.  To  get  the  maximum  camber  for  Cj  /  s  different  than 
l .0 ,  one  need  only  multiply  the  maximum  camber  for  the  C j  =  1.0  case  by  the  C j  in 
question.  Thus,  for  the  eyample  where  C/^  =  0.2,  the  maximum  camber  would  be  0.01474c 
(i.e.,  0.2  x  0.0737)  or  1.474  percent. 

Using  the  determination  criteria  above  for  camber  and  location  of  maximum  camber, 
we  can  now  examine  the  effect  of  maximum  camber  located  aft  of  the  0.4c  location.  As 
Figure  8  shows,  for  zero  camber  there  is  zero  shift  in  a^.  Figure  9  shows  data  for 
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Figure  8.  NACA  652  -  018  (a  =  1.0) 
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SECTION  ANGLE  OF  ATTACK,  a0  (degrees) 


Figure  11.  NACA  653  -  418  (a  =  1.0) 

1.474  percent  camber  at  0.45c.  Here  the  rule  seems  to  work.  Figure  10  shows  data  for 
4.422  percent  camber  at  0.45c.  Again,  the  rule  applies.  Figure  11  shows  that  for  maxi¬ 
mum  camber  at  0.5c  and  camber  of  2.948  percent,  the  rule  still  holds. 

So,  according  to  data  for  standard  NACA  wing  sections,  the  rule  seems  quite  accurate 
for  maximum  camber  locations  between  0.4c  and  0.5c  from  the  leading  edge,  and  it  appears 
to  be  less  accurate  as  the  maximum  camber  point  nears  the  leading  edge. 

Examination  of  the  data  for  modified  NACA  airfoils  (Ref.  3)  and  special  airfoils 
like  the  GAW(l)  (Ref.  4)  seem  to  support  the  rule  of  thumb.  Further,  since  these  air¬ 
foils  have  maximum  camber  locations  near  the  mid-chord  location,  the  rule  of  thumb  seems 
to  generally  apply  for  wing  sections  whose  maximum  camber  is  near  the  mid-chord  location. 
Implications  of  maximum  camber  locations  further  aft  are  addressed  in  the  next  section. 
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III.  Theoretical  Examination  of  the  Effect  of  Camber  on 

A  theoretical  examination  of  the  rule  of  thumb  can  be  made  by  using  thin-airfoil 
theory.  In  this  theory  a  wing  is  represented  by  a  camber  line  whose  distribution  is 
given  by  z,  the  distance  from  the  chord  line  to  the  camber  line,  which  is  a  function 
of  x,  the  distance  along  the  chord  line  measured  from  the  leading  edge  (see  Figure  12). 


THICKNESS 


Figure  12.  Airfoil  Nomenclature  and  Same  Airfoil 
in  Thin-Airfoil  Theory 

This  camber  line  represents  the  wing,  and  the  angle  of  attack,  ci,  is  the  same  angle 
between  the  chord  line  and  the  relative  wing  as  for  the  actual  airfoil.  If  we  further 
represent  the  distance  x  in  terms  of  a  geometric  function  of  a  convenient  angle  6  so 
t  hat 


x  =  he  (1  -  cos  0)  ,  (1) 

where  c  is  again  the  chord  length;  then,  according  to  thin-airfoil  theory,  the  coeffi¬ 
cient  of  lift,  C/ ,  is  given  by  (Ref.  5): 

Cl  =  2*  [a  +  -J-  fo  -g-  (cose  -  1)  de]  ,  (2) 

where  dz/dx  is  the  slope  of  the  camber  line  at  the  x  location  corresponding  to  0.  Clear¬ 
ly  Eqn.  (2)  tells  us  that  the  zero-lift  angle  of  attack,  (i.e.,  the  angle  of  attack 
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at  which  Ci  =  0)  ,  must  be 

a  - - —  f"  -4 —  (cosO  -  1)  d<-  (3) 

LO  tt  J  o  dx 

Thus,  if  we  know  z  as  a  function  of  x  (thus  •)  we  can  cither  calculate  »  ^  analytically 
from  Eqn.  (3)  or,  if  we  are  given  z  in  tabular  form,  we  can  calculate  ^  using  an  appro¬ 
priate  numerical  approximation  of  Eqn.  (3). 

The  simplest  case  to  examine  is  the  case  for  a  parabolic  camber  distribution  centered 
at  xc  equal  to  0.5c  (the  analysis  is  reproduced  here  from  Ref.  5  for  tutorial  purposes) 
so  that 

2  =  4  zc[t  -  (f)2]  (4) 

Eqn. (4)  is  represented  graphically  in  Figure  13.  By  changing  z c  we  can  represent  all 


percent  cambers  of  this  camber  distribution.  The  slope  of  the  camber  line  is  obtained 
by  differentiating  Eqn.  (4): 


dz 

dx 


(1  - 


c 


) 


or  by  substituting  x  in  terms  of  0  by  using  Eqn.  (1): 


dz 

dx 


=  4 


zc 

— —  cosO 
c 


If  we  then  substitute  Eqn.  (6)  into  Eqn.  (3),  we  have: 


a,  _  =  -  — —  f  4  — —  cosh  (cos6  -1)  do 

IAJ  TT  Jo  C 


(5) 


(6) 


(7) 


Finally,  when  we  integrate  Eqn.  (7),  we  obtain  a  as  a  function  of  zr  so  that 

LU  u 
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a 


LO 


(8) 


We  should  keep  in  mind  Chat  Eqn.  (8)  gives  u^g  in  radians,  so  that  in  the  generic  sense 
Eqn.  (8)  tells  us  that  o^g  is  -1.146  degrees/percent  of  camber.  For  this  camber  distri¬ 
bution  then,  the  rule  of  thumb  is  off  by  only  15  percent. 

We  can  also  generically  examine  the  rule  of  thumb  for  the  NACA  four-digit  series 
(on  which  no  wind  tunnel  data  is  readily  available)  by  using  a  crude  numerical  integra¬ 
tion  of  Eqn.  (3).  If  we  assume  that  the  camber  line  is  made  up  of  straight  segments 
arranged  according  to  tabular  camber  distributions  given  in  Ref.  2,  then  we  can  write 
Eqn.  (3)  as  follows: 


a 


LO 


-r  /I  "air  (cose  ‘  1}  de 


(cost*  -  1)  d6 


(9) 


or 


(10) 


Since  Ref.  2  gives  the  tabular  distributions  for  six-percent  camber  (all  the  others 
are  direct-ratio  multiples,  i.e.,  four  percent  would  require  multiplying  all  camber 
values  by  4/6),  1  used  these  as  the  generic  cases.  Using  Table  I,  1  computed  Uj  0  accord¬ 
ing  to  Eqn.  (10)  for  the  62 ,  63 ,  64 ,  65 ,  66 ,  and  67 distributions.  The  results 

were  -0.834  degrees/percent,  -0.8933  degrees/percent,  -0.9633  degrees/percent,  -1.039  de¬ 
grees/percent,  -1.138  degrees/percent,  and  -1.204  degrees/percent  for  maximum  camber  at 
0.2c,  0.3c,  0.4c,  0.5c,  0.6c,  and  0.7c  respectively.  Notice  that  the  general  trend  we 
observed  earlier  when  examining  experimental  data  of  decreasing  accuracy  of  the  rule 
as  the  maximum  camber  location  moved  forward  is  preserved  in  the  theoretical  calculations. 
Further,  we  can  now  see  that  as  the  location  of  maximum  camber  moves  aft  of  the  mid¬ 
chord,  the  rule  again  begins  to  err.  The  rule  is  most  accurate  between  maximum  camber 
locations  of  0.4c  to  0.6c.  Since  this  is  the  location  of  maximum  camber  for  most  wings 
(Ref.  2),  it  seems  the  rule  is  good  for  most  wings. 


IV .  Applying  the  Rule  of  Thumb  to  the  Interpretation  of  Experimental  Data 

Some  years  ago  several  of  my  students  examined  the  aerodynamic  characteristics  of 
a  cusped  airfoil  (Ref.  6  and  7).  At  that  time  an  airfoil  known  as  a  Kline-Fogleman  had 
appeared  in  the  news  as  having  exceptional  aerodynamic  characteristics.  The  Kline- 
Fogleman  airfoil  apparently  is  any  airfoil  which  has  a  discontinuous  lower  surface  (Ref. 
8);  the  airfoil  the  students  tested  was  of  this  classification.  The  findings  were  not 
particularly  spectacular  or  even  important,  but  they  did  serve  to  demonstrate  the  power 
of  using  the  rule  of  thumb  to  interpret  trends  in  data  that  at  first  appear  somewhat  odd 
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The  wing  section  shown  in  figure  14  served  as  a  model  for  several  series  of  tests. 
The  section  was  constructed  by  modifying  an  NACA  0015  with  a  vertical  cut  located  at 
the  quarter-chord  point  and  going  from  the  lower  surface  to  the  chord  line.  This  cut 


Figure  14.  Cusped  Airfoil  Formed  By  Modifying 
an  NACA  0015  Wing  Section 


was  joined  by  a  cut  from  the  trailing  edge  along  the  chord  line.  In  this  paper  I  will 
refer  to  this  section  configuration  as  a  cusped  airfoil.  All  tests  were  accompanied 
by  a  test  of  an  unmodified  NACA  0015  airfoil  for  comparison  purposes. 

In  these  series  of  t-«sts  two  types  of  planforms  were  used.  The  first  was  a  finite 
span  on  which  endplates  could  be  attached.  This  planform  was  mounted  on  a  force  balance 
so  that  aerodynamic  forces  could  be  measured  directly.  Figure  15  shows  selected  data 
representative  of  the  coefficient  of  lift  versus  angle  of  attack  for  this  planform. 

d 


Figure  15.  Data  for  NACA  0013  and  Cusped  Airfoil  with  Endplates  (Ref.  6) 


36 


USAFA-TR-8 1 -4 


The  second  planform  spanned  the  entire  width  of  the  wind  tunnel  so  that  the  data 
obtained  could  be  interpreted  as  two-dimensional.  Pressure  ports  located  at  various 
points  distributed  along  the  chord  direction  on  the  upper  and  lower  surfaces  of  the  wing 
provided  tne  information  necessary  to  plot  the  coefficient  of  pressure,  Cp,  versus  chord- 
wise  location  as  distance  from  leading  edge,  x,  divided  by  chord  length,  c,  or  x/c .  By 
appropriate  integration  of  these  plots  one  of  the  results  we  obtained  was  coefficient- 
of-lift  versus  angle-of-attack  plots.  Figure  16  shows  a  plot  of  selected  data  represent¬ 
ative  of  these  data. 

O 


Figure  16.  Two-Dimensional  Data  for  NACA  0015 
and  Cusped  Airfoil  (Ref.  7) 


Rather  than  interpret  Figures  15  and  16  as  absolute  or  definitive,  it  is  instruc¬ 
tive  to  look  at  what  appears  to  be  odd  trends,  which  at  first  glance  seem  to  defy  a  sim¬ 
ple  explanation.  If  we  examine  Figure  15  we  see  that  the  cusped  airfoil  has  an  which 
is  negative.  From  our  own  rule  of  thumb  we  realize  that  this  implies  camber  (unless 
otherwise  stated,  wre  mean  positive  camber  when  we  say  camber).  In  fact,  we  might  even 
be  tempted  to  say  the  camber  is  about  2.7  percent;  but  for  the  present  let  us  not  specu¬ 
late  on  the  amount  of  camber.  The  supposition  that  the  wing  has  camber  comes  as  no  sur¬ 
prise  since  the  wing  actually  does  have  camber.  From  geometric  considerations  only,  the 
amount  of  camber  is  about  4.5  percent.  What  is  interesting  is  that  the  slope  of 
the  cusped  wing  is  steeper  than  that  of  the  NACA  0015.  One  might  also  note  that  the 
slope  of  the  cusped  wing  is  not  a  straight  line  but  rather  decreases  as  a  increases  (this 
can  be  seen  most  easily  by  placing  a  straight-edge  along  the  curve).  The  slope  is  only 

17 
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4.3  per  radian,  which  is  not  unusual  when  compared  with  the  two-dimensional  limit  of 
5.7  per  radian  (2tt  Is  the  theoretical  two-dimensional  slope,  but  in  practice  we  find 
that  the  slope  is  only  5.7).  I  say  not  unusual  in  the  sense  that  we  might  think  that 
the  cusped  configuration  somehow  interfers  with  the  performance  of  the  endplates. 

Things  become  really  interesting  when  we  look  at  Figure  16.  First,  as  a  disclaimer 
of  sorts,  let  me  start  by  admitting  the  difficulty  of  obtaining  data  by  integrating 
Cp  versus  x/c  plots,  as  evidenced  by  the  value  of  Cj  at  a  =  0  for  the  NACA  0015  data. 

So  there  is  bound  to  be  some  error  bar  which  I  have  not  included  on  the  plot.  However, 
for  the  sake  of  instilling  some  confidence  in  the  data,  let  me  point  out  a  few  positive 
aspects.  First,  the  plot  of  the  NACA  0015  duplicates  the  data  from  Ref.  2  for  the  same 
section  (which  includes  the  C/^  slope  of  5.7  per  radian).  Second,  the  zero-lift  inter¬ 
sect,  a^Q»  f°r  the  cusped  airfoil  is  close  to  that  of  Figure  15.  Finally,  the  decreas¬ 
ing  slope  with  increasing  angle  of  attack,  noted  for  the  data  of  Figure  15,  is  again 
duplicated  in  the  data  of  Figure  16.  Now  let  me  point  out  that  the  slope  of  the  cusped 
airfoil  exceeds  the  two-dimensional  limit!  (It  is,  in  fact,  as  high  as  8.2  per  radian 
at  a  =  -4  degrees,  which  exceeds  even  the  theoretical  limit  of  2ir  per  radian.)  One  might 
be  tempted  to  claim  we  have  stumbled  on  a  wing  which  defies  our  present  understanding 
of  aerodynamics. 

Let  us  see  how  our  rule  of  thumb  might  be  helpful  in  interpreting  what  is  actually 
going  on.  Figure  17,  which  is  an  enlarged  portion  of  Figure  16  with  lines  added  at  vari- 

U_ 
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Figure  17.  Enlarged  Portion  of  Figure  16  with  Lines  of  Constant  Slope  Added 
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ious  locations,  is  helpful  in  this  regard.  Each  of  the  additional  lines  has  a  slope 
5.7  per  radian,  which  is  the  practical  two-dimensional  limit.  Examining  the  data 
presented  in  Figure  16  in  light  of  our  rule  of  thumb  seems  to  imply  that  the  cusped  wing 
section  is  acting  as  a  variable  camber  airfoil,  with  camber  increasing  with  angle  of 
attack. 

Based  on  this  interpretation,  we  went  back  to  the  wind  tunnel  and  located  flow  re¬ 
attachment  points  by  attaching  tufts  along  the  chord-wise  direction  in  the  cusped  region 
of  the  wing.  Our  conclusion  was  that  the  attachment  point  moved  forward  with  increasing 
angle  of  attack.  If  we  could  interpret  the  streamline  of  the  attachment  point  as  a  vir¬ 
tual  wing  surface,  then  the  wing  was  indeed  acting  as  a  variable  camber  airfoil.  This 
explanation  seems  much  more  satisfying,  and  we  have  a  way  of  checking  the  validity 
of  the  virtual  surface  interpretation. 

Figure  18  shows  the  attachment  point  for  a  “  0  degrees.  Note  also  that  a  line 
is  faired  in  from  the  attachment  point  to  the  lower  surface  of  the  airfoil  just  at  the 
cut  forming  the  cusp.  The  actual  surface,  coupled  with  the  faired  line,  should  approxi¬ 
mate  the  virtual  wing  section  representing  the  airfoil  the  potential  flow  field  thinks 
it  is  encountering.  This  new  shape  allows  us  to  locate  the  virtual  camber  line  also 
shown  in  Figure  18.  The  maximum  camber  for  the  virtual  section  of  Figure  18  is  approxi¬ 
mately  3.5  percent  at  0.6c.  According  to  our  rule  of  thumb*,  this  would  correspond  to  an 
of  approximately  -4.0  degrees,  which  in  turn  would  imply  a  Cj  of  approximately  0.4  at 
a-  0  degrees.  This  compares  extremely  well  with  our  measured  C|  at  a  =  0  degrees.  Thus  the 
case  for  assuming  a  virtual  wing  section  seems  well  justified.  (This  experimental  justifica¬ 
tion  of  using  a  potential  flow  solution  around  a  virtual  wing  section  induced  by  the  separated 
flow  region  in  which  we  have  trapped  vortices  tends  to  lend  credibility  to  modeling  tech¬ 
niques  by  Rossow  -  Ref.  9.) 

CAMBER  LINE  OF 


Maximum  camber  of  virtual  wing  ~  3.5  percent  located  at  approximately  0.6c 

Figure  18.  Cusped  Wing  Section  Showing  Attachment  Point  for  a  «  0  Degrees 
with  Virtual  Wing  Section  Faired  In 


♦Recall  from  Section  TIT  that  for  camber  located  at  0.6c,  the  more  exact  rule  is 

-1.138  degrees/percent,  which  gives  3  -4  degrees. 
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We  have  examined  the  rule  of  thumb  which  says  that  the  zero-lift  angle  of  attack 
of  a  cambered  airfoii  shifts  to  the  left  by  approximately  the  magnitude  in  degrees  of 
the  percent  camber.  In  doing  so,  we  found  that  the  rule  can  be  applied  with  little  error 
to  maximum  cam1''1*-  locations  of  0.4c  to  0.6c  measured  from  the  leading  edge.  Further, 
we  can  extend  tfu  usefulness  of  the  rule  by  noting  that  if  u^q  exceeds  (in  magnitude) 
the  rule,  then  we  can  assume  the  maximum  camber  is  located  aft  of  0.6c;  and  if  is 

less  (in  magnitude)  than  that  predicted  by  the  rule,  we  can  assume  the  maximum  camber 
location  is  forward  of  0.4c.  I  think  it  is  clear  from  the  example  presented  here  that 
using  the  rule  of  thumb  in  the  interpretation  of  causes  for  unexpected  (or  expected) 
experimental  results  can  be  of  great  benefit. 
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EXPERIMENTAL  AERODYNAMIC  EFFECTS  OF  THREE-DIMENSIONAL 
SQUARE  CROSS-SECTION  MISSILES  AT  MODERATE  ANCLES  OF  ATTACK 

T.R.  Yechout*,  C.J.  Zollars**,  and  D.C.  Daniel*** 

Abstract 

Subsonic  aerodynamic  force  and  moment  characteristics  have  b*-**:i  determined  for  a 
parametric  family  of  missiles  with  square  cross-sections.  The  g.  iriables  con¬ 
sidered  in  this  experimental  study  were  body  cross-section  cornet  ,  nose  planform 

and  fin  planform.  Data  were  obtained  for  missiles  with  and  without  i.  for  angles  of 
attack  from  0  to  30  degrees  at  roll  angles  of  0,  22.5,  and  45  degrees.  Surface  oil- 
flow  patterns  were  also  recorded.  The  results  clearly  show  the  strong  influence  of  body 
cross-section  corner  radius  on  the  aerodynamic  forces  and  moments. 

I .  Introduction 

Recently  the  possibility  of  using  rectangular  or  square  cross-section  bodies  as 
missiles  or  submunition  dispensers  has  received  increasing  attention.  This  interest 
is  due  primarily  to  the  increased  packing  efficiency  that  results  from  the  rectangular 
shape  as  compared  to  missiles  with  circular  cross-sections.  This  packing  efficiency 
gain  is  due  to  two  factors:  (1)  the  ease  of  packing  rectangular  modular  components  and 
(2)  greater  useable  volume  for  a  given  frontal  area.  What  must  be  determined,  however, 
is  whether  this  well-recognized  packing  advantage  is  out-weighed  by  possible  detrimental 
aerodynamic  effects.  Unfortunately,  only  a  limited  amount  of  aerodynamic  information 
exists  which  could  support  preliminary  design  or  analyses  of  missiles  with  square  cross- 
sections,  especially  those  requiring  moderate  amounts  of  maneuverability.  These  limited 
studies  include  the  two-dimensional  work  by  Pohlhamus  (Ref.  1  and  2)  and  the  more  recent, 
very  high  angle  of  attack  work  by  Clarkson  et  al  (Ref.  3). 

To  overcome  this  data  deficiency,  the  Air  Force  Armament  Laboratorv  and  the  Air 
Force  Academy  are  conducting  a  two-phase  research  program  to  analyze  the  aerodynamic 
characteristics  of  missiles  having  square  cross-sections.  The  first  phase  of  the  program 
is  concerned  with  measuring  the  experimental  forces  and  moments  on  square  missile  shapes 
for  a  wide  range  of  geometries,  orientations,  and  flow  conditions.  Pliotogr  .t|h  u  documen¬ 
tation  of  surface  oil-flow  patterns  which  define  the  shear  stress  distribution  along 
the  missile  surface  are  also  recorded.  The  second  phase  of  the  program  is  concerned 
with  detailed  flow-field  measurements.  This  paper  presents  partial  results  from  the 
first  phase  of  the  above  program.  Included  is  a  brief  discussion  of  the  wind  tunnel 
facility  used  in  the  investigation,  information  on  the  wind  tunnel  models,  major  force 
and  moment  results,  representative  oil-flow  results,  and  detailed  discussion  of  the  results. 

II .  Wind  Tunnel  Facility 

The  data  for  this  investigation  were  obtained  in  the  subsonic  wind  tunnel  of  the 

*Major,  USAF,  Assistant  Professor  of  Aeronautics,  DFAN 
**Captain,  USAF,  Instructor  of  Aeronautics,  DFAN 
***Assoc iate  Fellow  AIAA,  Research  Manager,  Air  Force  Armament  Laboratory 
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United  States  Air  Force  Academy  (Ref.  4).  This  facility  is  a  continuous  flow,  closed-circuit 
facility  which  has  a  test  section  of  two  feet  by  three  feet  and  is  capable  of  operation 
at  a  Mach  number  range  of  0.04  <  M„  <  0.35  at  atmospheric  pressure.  At  the  maximum  oper¬ 
ating  condition,  the  tunnel  is  capable  of  obtaining  a  unit  Reynolds  number  of  1.6  million 
per  foot.  The  facility  is  used  for  research  projects  as  well  as  educational  experiments 
and  demonstrations. 

III.  Wind  Tunnel  Models 

The  various  model  components  used  in  this  investigation  are  shown  in  Figure  1.  This 
figure  depicts  the  four  bodies,  each  with  a  different  cross-section  corner  radius;  the 


two  noses,  a  blunted-tangent  ogive,  and  a  pointed-tangent  ogive;  and  one  each  of 
the  three  sets  of  fins  that  were  used.  As  illustrated  by  the  figure,  the  four  bodies 
used  ranged  from  perfectly  round  (B4)  to  perfectly  square  (Bl).  All  components  were 
made  of  aluminum. 

The  normalized  body  corner  radii  (r/b)  values  investigated  were  0.0,  0.1,  0.2,  and 
0.5.  These  bodies  are  depicted  in  Figure  1  as  Bl,  B2,  B3,  and  B4  respectively.  Each 
body  section  was  12  inches  long  and  had  a  width  (diameter)  of  2  inches. 

The  blunted-tangent  ogive  nose  (Nl)  shown  in  Figure  1  was  3  inches  in  length.  The 
pointed-tangent  ogive  nose  (N2)  was  4  inches  long. 

The  three  sets  of  fins  were  designated  as  FI,  F2,  and  F3.  They  had  aspect  ratios 
(based  on  exposed  semi-span)  of  0.47,  0.67,  and  0.34  respectively.  Fins  FI  and  F3  had 
an  exposed  semi-span  of  1.25  inches  and  root  lengths  of  4  inches  and  6  inches  respectively. 
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Fin  F2  had  an  exposed  semi-span  of  1.5  inches  and  a  root  length  of  3  inches. 

The  four  bodies,  two  noses,  and  three  sets  of  fins  were  used  to  make  various  config¬ 
urations  of  the  missile  model.  Each  configuration  was  then  labeled  according  to  the  compo 
nent  nomenclature  shown  in  Figure  1.  For  instance,  the  square-body,  blunt-nose,  and 
three-inch  length  fin  configuration  was  noted  as  B1N1F2.  An  assembled  model,  mounted 
in  the  wind  tunnel,  is  shown  in  Figures  2a  and  2b.  As  can  be  seen  from  Figure  2b,  the 
fins  are  mounted  in  a  standard  cruciform  arrangement  at  the  corners  of  the  body. 

(a)  Finless  (b)  Finned 


Figure  2.  Finless  and  Finned  Wind  Tunnel  Model 


In  addition  to  the  models  described  above,  a  smaller  half-scale  version  of  model 
N1B3F1  was  also  fabricated.  This  model  was  used  for  scale  comparison  experiments  in 
the  subsonic  wind  tunnel,  and  it  also  will  be  used  in  future  experiments  to  investigate 
Reynolds  number  and  Mach  number  effects. 

IV.  Results 

The  primary  results  obtained  were  six  component  forces  and  moments  which  were  taken 
using  an  internally-mounted  strain  gauge  balance.  The  accuracy  of  the  balance  was  plus 
or  minus  0.5  percent  of  full  load  output.  It  was  mounted  6.97  inches  forward  of  the 
rear  of  the  model  (3.428  inches  forward  in  the  half-scale  model).  The  aerodynamic  force 
and  moment  coefficients  were  calculated  using  the  diameter  and  area  of  the  circular  cross 
section  body  (B4). 

The  six-component  data  were  taken  in  two-degree  increments  from  0  to  30  degrees 
total  angle  of  attack  at  roll  orientations  of  0,  22.5,  and  45  degrees.  The  free  stream 
Mach  number  was  approximately  0.3  for  all  tests;  the  unit  Reynolds  number  was  approximate 
ly  1.3  million  per  foot  for  all  tests. 

The  aerodynamic  forces  and  moments  on  the  model  without  fins  were  obtained  for  all 
nose,  body,  and  orientation  combinations.  For  the  model  with  fins,  data  were  obtained 
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for  configurations  using  nose  N1  and  fin  FI  with  all  bodies.  All  three-fin  configurations 
were  also  tested  using  nose  N1  and  body  B3  at  all  orientation  combinations. 

In  addition  to  the  force  and  moment  data,  surface  oil-flow  patterns  were  obtained 
for  various  body-fin  combinations.  These  data  were  taken  for  total  angles  of  attack 
of  15  and  30  degrees  at  roll  orientations  of  0,  22.5,  and  45  degrees. 

V.  Discussion  of  Results 

The  results  of  these  wind  tunnel  tests  are  presented  and  discussed  in  the  following 
paragraphs . 


A.  Aerodynamic  Force  and  Moment  Results 

Missile  Without  Fins:  These  force  and  moment  results  are  shown  in  Figures  3 
through  10.  Figure  3  shows  the  variation  of  normal-force  coefficient  with  angle  of  at¬ 
tack  for  each  of  the  four  bodies  with  the  blunt  nose.  As  can  be  seen  from  this  figure, . 


Figure  3.  Variation  of  Normal-Force  Coefficient 
with  Angle  of  Attack  for  Finless  Body 


the  effects  of  body  cross-section  corner  radius  are  very  noticeable,  particularly  above 
10  degrees  angle  of  attack.  The  normal-force  coefficient  for  the  body  with  the  perfectly 
square  cross-section  (Bl)  is  approximately  125  percent  higher  than  that  for  the  circular 
body  (B4)  at  28  degrees  angle  of  attack.  It  is  also  interesting  to  note  that  at  the 
higher  angles  of  attack  the  normal-force  coefficients  of  the  body  with  the  normalized 
corner  radii  of  0.2  (B3)  are  close  in  value  to  the  circular  body  (B4) ,  whereas  the  values 
for  the  body  with  normalized  corner  radii  of  0.1  (B2)  are  closer  to  the  perfectly  square 
body  (Bl).  This  nonlinear  behavior  is  depicted  further  in  Figure  4,  which  shows  the 
variation  of  normal-force  coefficient  with  body  corner  radius  for  discrete  angles  of 
attack  of  10,  12,  and  26  degrees. 
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Figure  4.  Variation  of  Normal-Force  Coefficient 
with  Body  Corner  Radius  for  Finless  Body 

The  body  cross-section  corner  radius  also  has  a  slight  effect  on  axial-force  coeffi¬ 
cient  as  shown  in  Figure  5.  Here,  the  axial-force  coefficient  decreases  with  increasing 
corner  radius.  This  is  due  primarily  to  the  fact  that  the  frontal  area  is  also  decreas¬ 
ing  with  increasing  cross-section  corner  radius. 


Figure  5.  Variation  of  Axial-Force  Coefficient 
with  Body  Corner  Radius  for  Finless  Body 

The  variation  of  side-force  coefficient  with  angle  of  attack  is  shown  in  Figure  6. 

As  would  be  expected  for  circuLar  cross-section  bodies  of  the  fineness  ratio  investigated 
here,  there  is  little  evidence  of  side-force  build-up  at  the  higher  angles  of  attack 
with  the  exception  of  body  NIB1.  This  is  the  perfectly  square  body  (r/b  =  0.0).  It 
has  a  side-force  coefficient  of  approximately  0.5  at  30  degrees  angle  of  attack.  This 
value  is  approximately  10  percent  of  the  normal-force  coefficient  generated  by  this  body 
at  that  angle  of  attack  and  is  felt  to  be  of  little  significance. 

An  interesting  scale  effect  result  is  given  in  Figure  7.  This  figure  depicts  the 
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Ct  (DEGREES! 

Figure  6.  Variation  of  Side-Force  Coefficient 
with  Angle  of  Attack  for  Finless  Body 


N1B3  (H«I1  Sul*) 


a  (DEGREES! 

Figure  7.  Scale  Effects  on  Normal-Force  Coefficient 
for  Finless  Body 


variation  of  finless  missile  normal-force  coefficient  with  angle  of  attack  for  body  NIB3 
and  the  half-scale  equivalent  of  this  body.  The  normal-force  coefficient  for  the  half¬ 
scale  model  is  increasingly  larger  from  approximately  6  degrees  angle  of  attack  upward. 
This  is  probably  due  to  the  fact  that  the  separating  boundary  layer  does  not  scale  geo¬ 
metrically  as  the  models  do. 

Some  effects  of  nose  planform  on  the  aerodynamics  are  shown  in  Figures  8,  9,  and 
10.  Figure  8  shows  the  variation  in  normal-force  coefficient  with  angle  of  attack  for 
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Figure  8.  Nose  Planform  Effects  on  Normal-Force  Coefficients 
for  Finless  Body 


a  (DEGREES) 

Figure  9.  Nose  Planform  Effects  on  Axial-Force  Coefficient 
for  Finless  Body 

the  blunted-  and  pointed-tangent  ogive  noses  affixed  to  body  B3.  The  longer  pointed 
nose  results  in  a  slightly  higher  normal-force  at  the  higher  angles  of  attack.  Figure 
9  shows  the  variation  in  axial-force  coefficient  with  angle  of  attack  for  the  two  noses. 
The  trends  are  as  expected,  with  the  sharp-nosed  body  having  a  lower  value  at  zero  angle 
of  attack,  and  the  difference  between  them  diminishing  as  angle  of  attack  increases. 
Finally,  the  variation  in  side-force  coefficient  with  angle  of  attack  for  the  two  noses 
is  shown  in  Figure  10.  There  is  essentially  no  difference  in  the  two  results. 
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Figure  10.  Nose  Planform  Effects  on  Side-Force  Coefficient 
for  Finless  Body 

Missile  With  Fins:  The  force  and  moment  test  results  are  shown  in  Figures  11  through 
14.  Figure  11  shows  the  variation  of  normal-force  coefficient  wich  angle  of  attack  for 
each  ‘‘in  combination  affixed  to  the  same  body.  The  finless  missile  values  of  normal- 
force  coefficient  are  also  shown  for  reference.  As  can  be  seen  from  the  results  in  Fig¬ 
ure  11,  the  finned  combination  NIB3F3  produces  the  highest  value  of  normal-force  coeffi¬ 
cient  for  the  configurations  presented  in  this  figure.  This  is  due  to  the  fact  that 
fin  F3  has  the  largest  surface  area.  Fins  FI  and  F2,  having  essentially  the  same  area 
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Figure  11.  Variation  of  Normal-Force  Coefficient 
with  Angle  of  Attack  for  Body  B3  with  Fins 
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as  each  other,  produce  essentially  the  same  results,  albeit  lower  than  F3. 

Some  interesting  effects  of  the  body  cross-section  corner  radius  for  the  missile- 
with-fins  combinations  are  shown  in  Figure  12.  In  this  case,  results  are  presented  for 


a  (DEGREES) 

Figure  12.  Effect  of  Corne  Radius  on  Normal-Force  Coefficient 
for  Bodies  with  Fin  Set  FI 

configurations  having  the  same  nose  (N 1 )  and  fins  (FI),  but  affixed  to  bodies  of  differ¬ 
ent  cross-section  corner  radius.  As  was  the  case  for  the  finless  missile  results,  the 
perfectly  square  body  (Bl)  produces  higher  values  of  normal-force  coefficient  at  the 
higher  angles  of  attack  than  the  other  bodies.  It  is  very  interesting  to  note,  however, 
that  the  variation  from  body  to  body  is  not  nearly  as  pronounced  as  it  was  in  the  finless 
missile  case  (Figure  3).  The  fins,  being  a  major  lift-producing  device,  tend  to  reduce 
significantly  the  effects  of  corner  radius  on  normal-force  coefficient  for  body-fin 


Figure  13.  Variation  of  Normal-Force  Coefficient 
with  Corner  Radius  for  Bodies  with  Fins 
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geometries.  This  is  demonstrated  further  in  Figure  13,  where  variations  of  finned  mis¬ 
sile  normal-force  with  body  cross-section  corner  radius  are  shown  for  three  discrete 
angles  of  attack. 

Rolling-moment  test  results  are  presented  in  Figures  14a,  b,  and  c.  These  results 
show  the  variation  in  rolling-moment  coefficient  with  total  angle  of  attack  for  fixed 
roll  angles  of  0,  22.5,  and  45  degrees.  The  vehicle  geometries  considered  are  again 
the  nose  N1  and  fins  FI  affixed  to  bodies  of  different  cross-section  corner  radius.  As 
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Figure  14a.  Variation  of  Rolling -Moment  Coefficient 
with  Total  Angle  of  Attack  (4>  =  0.0°)  for  Bodies  with  Fin  Set  FI 

can  be  seen  from  these  results,  the  corner  radii  have  some  small  effect  on  the  rolling 
moment  at  the  higher  angles  of  attack  when  the  vehicle  has  a  symmetric  roll  orientation, 
i.e.,  0  or  45  degrees.  The  effect  for  the  asymmetric  roll  orientation  (<J>  =  22.5  degrees), 
as  shown  in  Figure  14b,  is  pronounced.  For  this  condition,  the  rolling  moment  for  the 
configurations  with  the  more  square  bodies  (N1B1F1  and  N1B2F1)  is  large  at  the  high  angles 
of  attack  and  has  a  positive  sign.  The  rolling-moment  coefficient  for  configuration 
NIB3F1,  which  has  a  normalized  body  cross-section  corner  radius  of  0.20,  stays  near  zero 
throughout  the  angle  of  attack  range.  The  configuration  with  the  circular  body  cross- 
section,  N1B4F1,  has  a  rolling  moment  completely  unlike  the  others;  at  the  higher  angles 
of  attack  it  has  large  negative  values.  This  interesting  behavior  is  probably  related 
to  flow  separation  from  the  model  at  these  high  angles  of  attack  and  by  the  immersion 
of  the  fins  in  the  leeside  body  vortices.  In  fact,  as  will  be  shown  in  the  next  section, 
the  magnitude  of  the  corner  radius  does  strongly  affect  the  flow  separation  line(s)  on 
the  body,  and  hence  the  circumferential  location  of  the  body  vortices. 
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Figure  14b.  Variation  of  Rolling-Moment  Coefficient 
with  Total  Angle  of  Attack  (  <p  =  22.5°)  for  Bodies  with  Fin  Set  FI 


Figure  14c.  Variation  of  Rolling-Moment  Coefficient 
with  Total  Angle  of  Attack  (  <J)  =  45°)  for  Bodies  with  Fin  et  FI 

B.  Surface  Oil-Flow  Results 

Surface  oil-flow  patterns  were  obtained  for  various  body-fin  configurations 
during  the  course  of  the  experiments  described  in  this  paper.  Some  of  the  more  signifi- 


cant  of  these  patterns  are  shown  in  Figures  15  and  !.'>. 

Figure  15  shows  Che  oil-flow  patterns  that  were  obtained  for  the  configurations 
with  nose  Nl  and  fins  FI  affixed  to  the  four  bodies  of  different  cross-section  corner 
radius.  The  view  of  this  figure  is  of  the  side  of  the  model  that  would  have  been  toward 
the  observer  when  the  wind  tunnel  test  was  run.  The  angle  of  attack  for  these  runs  was 
30  degrees.  Review  of  this  figure  reveals  that  the  separation  patterns  for  configura¬ 
tions  NIB1F1  and  N1B2F1  (Figure  14a,  b)  are  very  similar  forward  of  the  fins.  The  flow 

(a)  (b) 


Figure  15.  Oil  Flow  Patterns.  (a)  tt  =  30°,  p  =  0°  (N1B1F1); 
(b)  <t  =  30°,  p  =  0°  (N1B2F1);  (c)  rit  =  30°,  $  =  0°  (NIB3F1); 
(d)  at  =  30°,  4>  =  0°  (N1B4F1) 


seems  to  separate  immediately  as  it  comes  around  the  bottom  corners;  at  that  point  there 
is  a  clearly  defined  counter-flow  region  on  the  side  of  each  configuration,  followed 
by  a  reattachment  line  which  moves  diagonally  down  the  body  from  the  bott-om  surface, 
near  the  shoulder  of  the  nose,  to  the  top  surface  just  in  front  of  the  fins.  The  flow 
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along  the  side  of  configuration  N1BJKI  is  not  at  all  like  that  of  the  first  two  configur¬ 
ations.  in  this  case,  the  flow  seems  to  separate  and  reattach  almost  immediately  as 
it  comes  around  the  bottom  corner.  The  separation  pattern  for  the  configuration  with 
the  circular  cross-section,  N1B4F1,  is  as  expected. 

Figure  16  shows  oil-flow  patterns  that  were  obtained  for  the  same  configurations 
and  total  angle  of  attack  as  above,  but  lor  a  roll  angle  of  22.5  degrees.  The  side  of 
the  model  that  is  shown  in  this  figure  was  photographed  at  the  completion  of  the  test 
after  the  model  had  been  rolled  157.5  degrees  to  facilitate  the  photography.  To  help 
eliminate  confusion,  a  schematic  of  the  cross-flow  plane  during  and  after  the  experiment 
is  shown  in  Figure  17.  it  can  be  seen  from  Figure  16a  and  b  that  these  two  bodies  again 
create  very  similar  separation  patterns.  From  this  one  would  expect  the  rolling  moments 
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Figure  17.  Roll  Orientation  for  Oil  Flow  Patterns  ( <(>  =  22.5°) 

is  the  case.  The  separation  pattern  for  configuration  N1B3F1,  as  shown  in  Figure  16c, 
is  quite  different  from  the  first  two  configurations.  The  separation  line  is  not  nearly 
as  clearly  defined  and  appears  to  have  moved  further  around  the  body.  The  net  result 
as  shown  in  Figure  14b  is  a  near-zero  rolling  moment.  This  could  be  due  to  the  fact 
that  the  body  vortex  essentially  passes  between  fins  1  and  2  as  shown  in  Figure  17.  The 
final  oil-flow  results  are  shown  in  Figure  16d  for  the  configuration  with  the  circular 
cross-section,  ri^re  the  typical  twin  separation  lines  appear  for  bodies  of  this  type. 

This  totally  different  separation  pattern  has  resulted  in  completely  different  rolling 
moment  characteristics  for  this  body  at  this  particular  orientation.  This  is  no  doubt  due 
to  completely  different  fin  loadings  that  resulted  from  the  shift  in  the  separation  lines. 

VI.  Conclusion 

Aerodynamic  force,  moment,  and  oil-flow  results  have  been  presented  for  a  parametric 
family  of  missiles  with  square  cross-sections.  Particular  attention  has  been  given  to 
the  effects  of  body  cross-section  corner  radius.  This  parameter  has  been  shown  to  have 
a  very  strong  effect  on  finless  missile  normal-force  coefficient  and  a  lesser,  but  still 
noticeable  effect  on  finned  missile  normal-force  coefficient.  Body  cross-section  corner 
radius  was  also  shown  to  have  a  strong  effect  on  the  circumferential  location  of  the 
flow  separation  line  on  the  body  at  combined-roll  and  high  angle-of-attack  orientations. 
The  location  of  this  separation  line  can  in  turn  affect  both  the  magnitude  and  sense 
of  the  rolling-moment  coefficient. 

The  results  presented  in  this  paper  are  part  of  a  continuing  research  program  de¬ 
voted  to  aerodynamics  of  missiles  with  square  cross-sections.  Future  research  will  em¬ 
phasize  Reynolds  number  effects  and  flow-field  measurements. 
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Symbols 

b  Body  Cross-Section  Width  (Square  Missile) 

Body  Cross-Section  Diameter  (Round  Missile) 

CA  Axial-Force  Coefficient 

Cfj  Normal-Force  Coefficient 

Cy  Side-Force  Coefficient 

Cj^  Rolling-Moment  Coefficient 

C,,,  Pitching-Moment  Coefficient 

Cn  Yawing-Moment  Coefficient 

M^,  Free  Stream  Mach  Number 

r  Body  Cross-Section  Corner  Radius 

V,,,  Free  Stream  Velocity 

a  Angle  of  Attack 

at  Total  Angle  of  Attack 

<f>  Roll  Angle 
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CALIBRATION  OF  SEVEN-HOLE  PROBES 
SUITABLE  FOR  HIGH  ANGLES  IN  SUBSONIC  COMPRESSIBLE  FLOWS 

A. A.  Gerner  and  C.L.  Maurer* 

Abstract 

This  paper  illustrates,  by  example,  a  method  for  calibrating  seven-hole  probes  to 
measure  local  total  and  static  pressures  and  relative  flow  angles  of  up  to  70  degrees 
in  subsonic  compressible  flows.  To  conserve  air  in  our  blowdown  wind  tunnel,  we  used 
the  method  of  Latin  Squares  to  statistically  sample  a  large  and  otherwise  unmanageable 
data  set,  thereby  reducing  to  a  minimum  the  number  of  data  points  required  to  construct 
a  polynomial  fit  to  the  data.  The  three-variable  third  order  polynomials  found  to  repre¬ 
sent  the  probe  calibration  permit  all  the  desired  output  quantities  to  be  found  explicit¬ 
ly  from  pressures  measured  on  the  probe  in  an  unknown  flow  field.  This  method  determines 
the  flow  angle  to  within  ±2  degrees  with  95  percent  certainty. 

I.  Introduction 

Many  present  and  future  aircraft  designs  are  beginning  to  employ  such  devices  as 
leading-edge  strakes,  forward-swept  wings,  and  canards.  These  devices  have  demon¬ 
strated  the  potential  for  enhancing  aircraft  maneuverability  and  control  by  producing 
strong  vortices.  Some  modern  aircraft,  such  as  the  Concorde,  actually  rely  on  vortices 
and  complex  flows  to  create  lift.  However,  in  some  instances,  primary  and  secondary 
vortices  can  interact  unfavorably,  causing  separation  and  loss  of  lift  on  portions  of 
a  wing.  To  be  sure,  these  vortex-laden  flow  fields  ate  quite  intricate  and  difficult 
to  analyze.  Flow  visualization  techniques  offer  a  way  to  gain  insight  into  vortex  inter¬ 
actions,  but  suffer  from  their  inability  to  provide  quantitative  information.  To  over¬ 
come  this  limitation,  small  probes  can  be  inserted  directly  into  the  flow  stream  to  gather 
meaningful  pressure  information.  Historically,  non-nulling  five-hole  probes  have  been 
used  to  determine  local  total  and  static  pressures  at  a  particular  point  in  a  flow,  as 
well  as  flow  directions  up  to  40  degrees  relative  to  their  axis.  Nevertheless,  it  is 
not  inconceivable  for  the  local  flow  angles  of  strong  vortices  to  exceed  60  degrees. 

For  this  reason,  the  Air  Force  Academy,  under  a  grant  from  NASA-Ames,  has  developed  a 
unique  seven-hole  probe.  In  addition  to  local  total  and  static  pressure  measurements, 
these  probes  have  demonstrated  the  ability  to  determine  flow  angles  up  to  80  degrees 
relative  to  their  axis.  When  combined  with  a  computerized  data  acquisition  system,  they 
are  capable  of  taking  data  at  a  rate  of  nearly  two  data  points  per  second,  much  faster 
than  nulling  probes  which  require  considerable  time  to  balance  probe  tip  pressures  be¬ 
fore  each  pressure  measurement  can  be  taken. 

In  addition,  these  probes  are  very  small  (about  one-tenth  of  an  inch  in  diameter), 
so  they  do  not  significantly  disturb  the  flow  they  are  measuring.  But  because  of  this 
small  size,  they  suffer  from  inherent  manufacturing  defects.  As  a  result,  each  probe 
must  be  calibrated  before  it  can  be  a  useful  measuring  device.  Gallington  describes 
such  a  calibration  procedure  (Ref.  1)  which  is  both  fast  and  effective,  but  is  only  valid 
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for  Incompressible  flows.  This  constraint  effectively  restricts  the  probes  to  surveying 
flows  with  free-stream  Mach  numbers  of  0.3  and  below  (Ref.  2).  However,  the  mere  presence 
of  a  body  in  an  airflow  causes  the  flow  to  accelerate  in  certain  regions,  resulting  in 
local  velocities  greater  than  the  free-stream  velocity.  In  addition,  local  flow  velocities 
in  vortex  fields  are  likely  to  exceed  those  of  the  free  stream  too.  As  a  result,  even 
though  the  free-stream  conditions  are  slow  enough  to  justify  an  incompressible  flow  as¬ 
sumption,  the  local  flow  conditions  might  in  fact  be  compressible.  The  purpose  of  this 
report  then  is  to  develop  a  power  series  calibration  scheme  which  accurately  determines 
the  actual  flow  conditions  from  pressures  measured  on  the  probe  in  compressible  flows. 

Because  of  the  similarities  in  compressible  and  incompressible  theory,  we  begin 
our  discussion  by  developing  fully  the  incompressible  calibration  theory.  Then,  by  anal¬ 
ogy,  we  expand  this  theory  to  obtain  the  desired  form  of  the  calibration  suitable  for 
compressible  flow.  Next,  we  describe  the  apparatus  and  procedures  we  used  to  calibrate 
a  probe  in  subsonic  compressible  flow,  and  finally,  we  discuss  the  results  of  that  cali¬ 
bration. 

II.  Incompressible  Flow  Theory 

To  simplify  our  discussion  on  the  calibration  of  seven-hole  probes  in  compressible 
flow,  we  will  begin  with  a  background  in  the  theory  of  incompressible  flow  calibrations. 
Much  of  the  information  presented  in  the  following  sections  is  drawn  from  the  work  of 
Gallington  (Ref.  1).  Readers  wishing  further  information  on  the  theory  of  incompressible 
probe  calibrations  are  highly  encouraged  to  consult  this  reference. 

Axis  System  for  Low  Flow  Angles 

To  begin  our  discussion  of  calibration  theory,  we  restrict  our  treatment  to  low 
flow  angles;  typically,  those  values  for  which  the  angle  between  the  velocity  vector 
of  the  flow  stream  and  the  probe's  axis  are  less  than  30  degrees.  The  more  familiar 
reference  system  for  low  flow  angles  measures  velocity  vectors  in  terms  of  the  angle 
of  attack,  a,  and  the  angle  of  sideslip,  8.  However,  in  choosing  our  reference  system, 
we  adopted  the  tangential  reference  system  illustrated  in  Figure  1.  In  this  system, 
is  taken  to  be  the  projection  on  the  vertical  plane  of  the  angle  between  the  velocity 
vector  and  the  probe's  axis.  And  to  preserve  symmetry,  is  defined  as  the  projection 

on  the  horizontal  plane  of  the  angle  between  the  probe's  axis  and  the  relative  wind. 

For  this  reason,  the  tangential  reference  system  differs  slightly  from  the  u-8  system, 
but  will  be  used  to  evaluate  the  low  angle  flow  properties. 

Pressure  Coefficients  for  Low  Flow  Angles 

For  low  flow  angles  it  is  desirable  to  define  dimensionless  pressure  coefficients 
which  utilize  all  seven  measured  probe  pressures  and  are  sensitive  to  changes  in  flow 
angularity  with  respect  to  the  probe's  x-axis.  From  Figure  2,  one  such  pressure  coeffi- 
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Figure  1.  Low  Angle  Reference  System 


cient  sensitive  to  changes  in  angle  of  attack  in  the  x-z  planes  is  defined  as: 

P4  -  Pi 


“l 


P7  -  Pl-6 


(1) 


where  the  numerator  measures  changes  in  flow  angularity  based  on  the  differences  in  oppo¬ 
site  port  pressures,  and  the  denominator  nondimensionalizes  the  term  with  the  apparent 
dynamic  pressure.  This  pseudo-dynamic  pressure  is  obtained  from  the  difference  between 
the  central  port  pressure,  P7,  which  approximates  the  total  pressure  at  low  angles,  and 
the  average  of  the  six  surrounding  pressures,  Pi_6 ,  which  collectively  approximates  the 
static  pressure.  From  the  definition  of  this  pressure  coefficient,  it  is  easy  to  see 
that  two  other  possibilities  also  exist:  one  which  measures  the  pressure  differential 
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Figure  2.  Port  Numbering  Convention  and  Principal  Axes 
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between  ports  three  and  six,  and  the  other,  which  measures  the  pu.„  ;ure  differential 
between  ports  two  and  five.  The  complete  set  of  these  pressure  coefiicients  include: 


C 

>n 


-  Pi 

P7  -  Pl-6 


c 

02 


Pi  -  P(, 

P?  -  Pl-6 


P?  -  Pi 

P 7  -  Pl-6 


(2) 


But  before  these  coefficients  can  be  of  any  use  to  us,  they  must  be  resolved  into  the 
a  -  l3.p  reference  system.  This  is  done  by  weighing  the  contribution  of  each  coefficient 
in  Eqns.  (2)  along  the  respective  axis,  which  results  in  the  following  two  equations: 


C  =  4-  (2C  +  C  -  C  ) 

Cl  3  CXI  0t2  013 


c,  =  — (C  +  C  ) 

0  /3-  “2 


(3) 


The  first  equation  defining  Cu  contains  all  three  coefficients  of  Eqns.  (2).  In  particu¬ 
lar,  C..(  has  the  greatest  significance,  which  only  makes  sense  since  it  lies  directly 
along  the  axis  of  interest.  The  equation  defining  Cg  only  takes  into  account  the  last 
two  coefficients  of  Eqns.  (2),  assigning  to  each  an  equal  significance.  The  fact  that 
C.ti  is  not  included  in  this  equation  again  makes  sense,  since  it  is  directly  aligned 
with  the  a^-direction,  ideally  making  it  insensitive  to  changes  in  the  perpendicular  B  - 
direction.  In  summary,  the  procedure  of  obtaining  Ca  and  C,,  requires  two  tasks.  First, 
determine  Cal ,  C,x , ,  and  Ca,  from  the  seven  measured  pressures  using  Eqns.  (2)  and  then 
substitute  these  intermediate  quantities  into  Eqns.  (3)  for  the  desired  coefficients. 

Having  defined  the  two  angular  pressure  coefficients,  it  is  now  appropriate  to  dis¬ 
cuss  the  remaining  low  angle  pressure  coefficients,  0o  and  C^,  defined  as: 
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P7  -  Pl-6 


(4) 


C0  is  the  apparent  total  pressure  coefficient  and  functions  as  a  correction  factor 
to  convert  actual  pressures  measured  by  the  probe  to  accurate  values  of  local  total  pres¬ 
sure.  From  the  numerator,  it  is  seen  that  the  coefficient  measures  the  difference  between 
the  pseudo-total  pressure  measured  by  the  probe,  P7 ,  and  the  actual  total  pressure.  Just 
as  with  Eqn.  (1),  the  coefficient  is  nondimens ionalized  by  the  denominator,  which  is 
a  measure  of  the  apparent  dynamic  pressure  as  previously  described. 

Cq  functions  much  like  C0,  but  instead  of  correcting  probe  pressures  to  total  pres¬ 
sure,  Cq  relates  these  pi  assures  to  the  actual  dynamic  pressure.  In  this  coefficient, 
the  numerator  represents  the  probe's  approximation  of  dynamic  pressure  while  the  denomi¬ 
nator  consists  of  the  actual  dynamic  pressure. 
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Axis  System  for  High  Flow  Angles 

Up  to  this  point,  our  discussion  has  been  limited  to  a  description  of  the  pressure 
coefficients  used  for  flow  angles  below  30  degrees.  Yet,  Lhe  real  advantage  of  using 
a  seven-hole  probe  lies  in  its  ability  to  determine  flow  angles  as  high  as  80  degrees 
to  the  probe's  x-axis. 

A  reference  system  better  suited  to  high  angle  measurement  than  the  tangential  sys¬ 
tem  is  the  polar  reference  system,  which  measures  flow  angularities  in  terms  of  0  and 
and  is  shown  in  Figure  3.  In  this  system,  ",  the  pitch  angle,  is  the  angle  the  velocity 
vector  makes  with  respect  to  the  probe's  x-axis;  and  ;,  the  roil  angle,  describes  the 
azmuthal  orientation  of  the  velocity  vector  in  the  y-z  plane,  measured  counterclockwise 
from  the  negative  z-axis  as  viewed  from  the  front.  Although  a  singularity  exists  direct¬ 
ly  along  the  x-axis,  this  does  not  represent  any  problems  in  high  angle  measurement  and 
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Figure  3.  High  Angle  Reference  System 

is  avoided  entirely  by  switching  to  the  tangential  reference  system  at  low  angles. 

Pressure  Coefficients  for  High  Flow  Angles 

At  low  angles  of  attack,  all  seven  of  the  measured  probe  pressures  are  used  to  form 
the  pressure  coefficients.  However,  at  high  angles  of  attack,  as  illustrated  in  Figure 
4,  the  flow  tends  to  detach  over  the  downstream  portions  of  the  probe.  Pressure  ports 
lying  in  this  separated  region  are  insensitive  to  changes  in  flow  angularity;  consequent¬ 
ly,  it  is  not  feasible  to  use  their  pressures  in  a  meaningful  coefficient.  As  a  result, 
the  pressure  coefficients  for  high  angle  measurements  must  be  defined  so  that  they  include 
only  the  pressures  from  ports  in  attached  flow. 

Typically,  the  separation  points  of  a  cylinder  in  turbulent  flow  are  over  100  degrees 
from  the  frontal  stagnation  point  (Ref.  3).  And  for  a  conical  body,  such  as  the  probe  tip, 
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Figure  4.  Flow  Over  Probe  at  High  Angle  of  Attack 

the  flow  is  likely  to  remain  attached  longer.  In  addition,  the  u-velocity  component 
is  also  likely  to  extend  the  separation  points.  Thus,  for  the  condition  depicted  in 
Figure  4,  pressure  ports  three,  four,  five,  and  seven  lie  in  reliably  attached  flow; 
port  one  is  in  separated  flow;  and  the  disposition  of  ports  two  and  six  is  uncertain. 
Using  only  those  ports  in  attached  flow,  a  coefficient  sensitive  to  the  angle  of  pitch 
is  defined  as: 


C 


P-  -  P? 


P>, 


P7  +  I'o 

2 


(5) 


Following  the  same  rationale  as  in  the  low  angle  case,  the  numerator  measures  changes 
in  9  based  on  the  differences  in  opposing  port  pressures.  In  this  example,  P7,  the  smal¬ 
ler  of  the  two  pressures,  is  subtracted  from  P4.  Here  again,  the  coefficient  is  nondimen- 
sionalized  by  dividing  through  with  the  apparent  dynamic  pressure.  This  pseudo-dynamic 
pressure  is  determined  from  the  difference  between  the  peripheral  port  pressure,  P4 , 
which  at  high  angles  approximates  the  total  pressure,  and  the  average  of  P3  and  P5 ,  which 
when  taken  together  are  relatively  independent  to  changes  in  roll  and  approximate  the 
static  pressure. 

Using  a  similar  argument,  a  coefficient  which  changes  in  proportion  to  roll  angle 
is  appropriately  defined: 


P3  +  P5  (6) 

2 


The  numerator  of  Eqn.  (6)  is  sensitive  to  changes  in  <P,  in  that  ,s  the  velocity  vector 
rolls  in  either  direction,  the  windward  pressure  rises  and  the  leeward  pressure  falls. 

In  this  way,  the  difference  between  the  two  pressures  varies  significantly  for  variations 
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in  roll,  yet  the  average  of  their  sums  remains  relatively  constant.  Once  more,  the  coef¬ 
ficient  is  nondimens ionalized  with  the  same  denominator  as  in  the  previous  case. 

Obviously,  the  above  two  coefficients  are  only  valid  for  a  narrow  range  in  roll 
angle  about  port  four.  That  is,  as  we  rotate  the  velocity  vector  to  either  side  of  port 
four,  the  region  of  separated  flow  approaches  either  port  five  or  three.  Therefore, 
to  insure  that  all  pressures  are  taken  from  ports  in  attached  flow,  we  restrict  Eqns.  (5) 
and  (6)  to  roughly  a  60-degree  pie-shaped  sector  centered  on  port  four.  In  this  way, 
six  pie-shaped  sectors  are  summarily  defined  for  high  angle  measurement,  such  that  each 
has  its  own  set  of  coefficients  based  on  the  pressures  in  attached  flow.  The  remaining 
angular  pressure  coefficients  are  defined  with  the  same  method  used  to  develop  Eqns. 

(5)  and  (6),  resulting  in  the  following  set  of  equations: 
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Similarly,  the  C0  and  Cq  coefficients  are  developed  with  the  same  rationale  used 
to  derive  their  low  angle  counterparts;  the  only  difference  resides  in  the  choice  of 
the  pressures  which  are  roughly  equivalent  to  total  and  static  pressures.  These  pressures, 
of  course,  vary  in  relation  to  the  sector  a  particular  coefficient  describes.  The  com¬ 
plete  set  of  these  coefficients  include: 
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(8) 


Division  of  Angular  Space 

Having  defined  a  host  of  coefficients  for  low  and  high  angles  and  for  various  sec¬ 
tors  around  the  probe,  a  question  arises  as  to  when  a  particular  set  of  coefficients 
should  be  used.  To  arbitrarily  assign  angular  cut-offs,  based  on  probe  symmetries,  would 
be  naive  since  actual  data  might  suggest  better  division  lines.  With  this  in  mind,  a 
better  scheme  for  locating  the  sector  division  lines  is  based  on  the  isobars  depicted 
in  Figure  5.  This  method  defines  seven  sectors,  the  central  low  angle  sector  and  the 
six  high  angle  periphery  sectors,  and  allocates  data  points  to  a  given  sector  based  on 
the  highest  port  pressure  measured  on  the  probe. 


Polynomial  Power  Series  Expansion 

Once  the  data  points  are  allocated  to  the  proper  sector  with  its  corresponding  pres¬ 
sure  coefficients,  a  fourth  order  polynomial  expansion  is  used  to  solve  for  the  desired 
quantities.  In  two  variables  (l.e.,  using  the  two  angular  pressure  coefficients)  this 
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SECTOR  1 


where  A  is  either  aT>  6T>  CQ  or  Cq  for  low  angles  and  0,  C0  or  Cq  for  high  angles, 
with  the  subscript  denoting  the  ith  such  quantity.  The  K's  are  the  calibration  coeffi¬ 
cients,  with  the  superscripts  denoting  the  quantity  to  which  a  particular  set  of 
K's  belong,  and  the  subscripts  identifying  the  coefficient  of  a  particular  term  in 
the  power  series  expansion.  Note  that  in  the  high  angle  case,  the  C^'s  and  Cg's  are 
replaced  by  Co's  and  C^'s  respectively.  In  matrix  notation  for  "n"  data  points  of  a 
particular  sector,  a  set  of  Eqn.  (9)'s  are  represented  as: 
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But  to  simplify  further  discussion,  Eqn.  (10)  is  rewritten  as: 

[A]  -  [C][K]  (11) 

In  this  form,  the  n  x  1  A-matrix  contains  n  values  of  the  flow  parameters  of  interest, 
the  n  x  15  C-matrix  contains  the  expanded  angular  pressure  coefficients  for  each  of  the 
n  flow  parameters,  and  the  15  x  1  K-matrix  contains  the  calibration  coefficients  for 
the  particular  flow  parameter  of  interest. 


Determining  the  Calibration  Coefficients 

During  the  calibration  process,  the  quantities  within  the  A-matrix  are  set  up  by 
the  experimenter  in  terms  of  the  known  tunnel  conditions,  and  the  terms  within  the  C- 
matrix  are  determined  from  the  measured  probe  pressures.  A  calibration  procedure,  there¬ 
fore,  involves  the  calculation  of  the  unknown  K-matrix.  This  calculation  is  performed 
by  rearranging  Eqn.  (11)  to  solve  for  the  unknown  calibration  coefficients.  With  matrix 
algebra,  this  is  performed  according  to  the  procedure  outlined  by  Netter  and  Wasserman 
(Ref.  4): 

First  multiply  each  side  of  Eqn.  (11)  by  the  transpose  of  the  C-matrix: 

[C]T  [A]  -  [C]T  [C][KJ  -  [CTC][K]  (12) 

Next  multiply  each  side  by  the  inverse  of  the  recently  created  C ^C-matrix: 

[CTC]-1  [C]T[A]  -  [CTcrl[CTCl[K]  (13) 

Realizing  that  the  product  of  a  matrix  and  its  inverse  results  in  the  identity  matrix, 
Eqn.  (13)  simplifies  to  yield  a  solution  for  the  unknown  K-matrix  in  terms  of  the  known 
C-  and  A-matrices: 


t 
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[K]  =  [CTC]_1[C]T[A] 

This  technique  determines  the  calibration  coefficients  by  a  least  squares  curve  fit  to 
the  experimental  data. 


Determining  the  Desired  Flow  Properties 

Once  the  calibration  coefficients  are  determined,  the  calibration  process  is  complete. 
The  probe  is  ready  to  be  inserted  in  an  unknown  flow  field  and  the  desired  flow  properties 
determined.  Once  in  the  flow  field,  the  probe's  measured  pressure  readings  allow  us 
to  determine  the  angular  pressure  coefficients;  these  coefficients  are  then  manipulated 
to  fill  the  C-matrix  of  Eqn.  (10).  Since  the  K-matrix  is  already  known,  the  desired 
flow  properties  in  the  A-matrix  are  then  determined  explicitly.  For  a  particular  condi¬ 
tion,  the  solutions  for  the  desired  flow  properties  take  on  the  following  functional  forms: 


Inner  Sector  (low  flow  angles) 
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Outer  Sectors  (high  flow  angles) 


e  '  f<W-K?  t,tSce„*K?V  •  • 


+  Ki sCi 
^n 


*  f(C6„'V*  K?  *  ^C9n+  *Kn+  ‘  •  • 


C°n  "  f(CV  <V  =  K,  n  +  K2  nC8n  +  K3  %n  + 


Cqn  . 


V 


Cqn  -  f(C6n.  C$n)=  Ka  +  K2  "C9n  +  K3  "C^n  + 


+  ^=c;n 


°n  u 
+  K]  5  Ca 


Mn  u 

+  Kl5  C5n 


(16) 


Once  the  local  total  and  dynamic  pressure  coefficients  are  specified,  it  is  possible  to 
determine  the  local  total  and  dynamic  pressures.  This  is  accomplished  by  rearranging  Eqns. 
(4)  in  the  case  of  the  inner  sector,  and  Eqns.  (8)  in  the  case  of  the  outer  sectors, 
thereby  solving  for  the  desired  pressures.  This  calculation  involves  the  polynomial 
result  for  CG  or  C^  and  the  seven  measured  probe  pressures.  As  an  example,  the  inner 
sector  equations  for  local  total  and  dynamic  pressures  are  derived  as  follows: 
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Recalling  that 
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the  local  total  and  dynamic  pressures  are  then  solved  by  manipulating  Eqns.  (4);  thus: 


P  ,  -  P?  -  C  (P7  -  P 1  — 6  ) 
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The  same  procedure  is  extended  to  determine  the  local  total  and  dynamic  pressures  in 
the  outer  sectors. 


III.  Extension  to  Compressible  Flow 

Up  to  this  point  we  have  discussed  the  form  of  a  calibration  procedure  which  is 
only  valid  for  incompressible  flow.  This  is  so  because  the  calibration  coefficients 
depend  directly  upon  the  angular  pressure  coefficients,  which  as  pressure  coefficients 
are  themselves  dependent  upon  Mach  number  (Ref.  5).  As  such,  seven-hole  probes  with 
the  present  method  of  calibration  are  limited  to  surveying  flows  within  the  incompres¬ 
sible  regime.  This  restriction  is  lifted  easily  enough  with  the  inclusion  of  an  addi¬ 
tional  pressure  coefficient  representative  of  compressibility.  Consequently,  the  number 
of  variables  in  the  polynomial  curve  fit  is  increased  from  two  (for  incompressible  flow) 
to  three  (for  compressible  flow). 

Determining  a  Coefficient  of  Compressibility 

Just  like  the  angular  pressure  coefficients  (C  and  C  or  C  and  C  ),  the  compres- 

Ct  p  W  <p 

sibility  coefficient  (hereafter  denoted  C^)  must  be  nondimensional  and  determined  strict¬ 
ly  from  pressures  measured  on  the  probe.  In  addition,  C^  must  have  the  feature  that 
it  approaches  zero  at  very  low  Mach  numbers.  That  is,  at  very  low  Mach  numbers  all  terms 
bearing  a  compressibility  coefficient  should  be  insignificant  in  the  power  series  expan¬ 
sion.  This,  in  effect,  brings  us  back  to  an  incompressible  power  series  expansion  in 
two  variables,  where  the  significant  terms  only  contain  the  angular  pressure  coefficients. 
A  further  constraint  requires  the  compressibility  coefficient  to  approach  a  finite  value 
in  the  hypersonic  limit.  In  other  words,  for  very  high  Mach  numbers,  large  changes  in 
Mach  number  should  have  a  neglible  effect  on  the  compressibility  coefficient.  This  re¬ 
flects  a  limitation  shared  by  all  pressure  probe  methods  in  that  Mach  number  becomes 
indeterminate  at  hypersonic  speeds  (Ref.  1).  These  requirements  are  satisfied  by  model¬ 
ing  the  compressibility  coefficient  after  the  dynamic  to  total  pressure  ratio  depicted 
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in  Figure  6.  This  figure  is  idealized  for  isentropic  flow;  as  such,  it  is  not  useful 
for  flow  speeds  much  beyond  Mach  one.  However,  since  our  calibration  is  limited  to  sub¬ 
sonic  speeds,  the  isentropic  idealization  is  valid  through  the  sonic  limit. 


Figure  6.  Compressibility  Effects  as  a  Function  of  Mach  Number 


To  develop  a  compressibility  coefficient  from  probe  measured  pressures,  we  need 
to  represent  the  total  and  dynamic  pressures  from  probe  measured  pressures.  For  the 
low  angle  case,  a  pseudo-total  pressure  is  approximated  by  P7  with  a  pseudo-dynamic  pres¬ 
sure  approximated  by  P7  -  P1-6.  Thus,  for  the  inner  sector,  a  compressibility  coefficient 
in  terms  of  the  dynamic  to  total  pressure  ratio  is  modeled  by: 


C 


M? 


P.7  -  P 1  —6 


P? 


(18) 


The  compressibility 
manner;  these  outer 


coefficients  for  each  of  the  outer  sectors  are  modeled  in  a  similar 
sector  coefficients  include: 

Pe  +  P2 


Pi 


CMi 
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M2 


Pi  +  Pj 


P2 


(19) 
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Selection  of  Data  Points 

Typical  incompressible  probe  calibrations  take  approximately  80  data  points  in  two 
variables  (C^  and  Cg)  for  each  of  the  seven  sectors  (Ref.  1).  This  results  in  a  total 
of  about  560  data  points  for  a  complete  calibration.  To  extend  this  present  scheme  into 
yet  another  dimension  (Mach  number)  would  create  a  data  set  of  intractable  proportions. 
Consequently,  it  is  necessary  to  represent  the  data  set  with  a  sample  of  more  manageable 
proportions.  In  addition,  this  sample  must  be  chosen  such  that  the  density  of  chosen 
data  points  throughout  the  data  set  is  homogenous.  In  other  words,  the  sample  must  be 
an  accurate  representation  of  the  data  set,  otherwise  the  calibration  routine  will  not 
offer  consistent  accuracy  throughout  the  range  of  data. 

A  method  of  ensuring  a  homogenous,  yet  random,  sample  of  a  three-dimensional  parame¬ 
ter  space  is  suggested  by  Cochran  and  Cox  (Ref.  6).  The  technique  is  known  as  the  method 
of  Latin  Squares,  an  example  of  which  is  shown  in  Figure  7.  The  plan  depicted  in  this 
figure  is  a  3  x  3  square  which  actually  represents  a  three-dimensional  parameter  space; 
one  variable  along  the  vertical  axis,  a  second  variable  along  the  horizontal  axis,  and 
the  third  variable  denoted  by  the  letters  along  the  axis  going  into  the  page.  This 
square  is  better  visualized  in  Figure  8,  where  it  is  shown  in  three  dimensions  instead 
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Figure  7.  3x3  Latin  Square  Plan 
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Figure  8.  3x3  Latin  Square  Shown  in  Three  Dimensions 

of  the  two-dimensional  rendition  of  Figure  7.  The  3x3x3  cube  represents  the  entire 
parameter  space,  and  each  sub-cube  represents  a  discrete  data  point  within  the  entire 
data  set.  In  Figure  7,  the  A' s  denote  the  data  points  to  be  selected  in  the  A-plane 
(i.e.,  the  plane  in  which  the  third  variable,  Mach  number,  is  held  constant),  the  B' s 
denote  the  points  selected  in  the  B-plane,  and  so  on.  In  this  way,  the  entire  data 
set  is  sampled  by  the  points  which  appear  as  the  shaded  cubes  in  Figure  8.  The  unique 
feature  of  this  sampling  technique  is  that  no  matter  from  which  axis  direction  the  cube 
is  viewed,  the  entire  frontal  space  will  appear  covered  by  data  points,  one  point  deep; 
that  is,  no  two  points  will  appear  to  overlap  each  other.  Furthermore,  in  any  given 
plane  where  one  variable  Is  held  constant,  every  value  of  the  remaining  two  variables 
is  sampled  exactly  once.  Consequently,  the  method  of  Latin  Squares  guarantees  a  homo¬ 
genous  density  of  data  points  comprising  the  sample  of  the  representative  data  set.  Of 
course,  the  technique  of  Latin  Squares  is  not  limited  to  a  3  x  3  plan,  but  can  be  ex¬ 
panded  to  a  12  x  12  plan  if  desired;  however,  the  more  common  plans  range  between  the 
5x5  and  8x8  squares  (Ref.  6).  This  is  true  because  as  plan  size  increases,  the  ratio 
of  points  sampled  to  total  points  in  the  data  set  decreases.  To  illustrate,  the  3x3 
square  of  Figure  7  samples  a  27-point  data  set  with  nine  points,  whereas  a  10  x  10  square 
uses  100  points  to  sample  a  1000-point  set.  Accordingly,  the  3x3  Latin  Square  samples 
33  percent  of  the  entire  data  set,  while  the  10  x  10  Latin  Square  samples  only  10  percent 
of  the  complete  set.  Nevertheless,  a  3  x  3  square  describes  a  given  parameter  space 
with  only  27  data  points,  while  a  10  x  10  square  divides  the  same  parameter  space  into 
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1000  discrete  data  points.  Thus,  a  larger  Latin  Square  represents  a  given  parameter 
space  with  greater  resolution,  but  samples  the  resulting  data  set  with  a  smaller  percen¬ 
tage  of  points.  The  experimenter,  therefore,  is  faced  with  a  compromise  in  choosing 
the  Latin  Square  best  suited  to  his  needs. 

For  reasons  to  be  discussed  later,  we  selected  a  6  x  6  Latin  Square  (shown  in  Figure 
9)  for  each  sector.  This  gives  us  a  data  set  with  six  different  values  for  each  of  the 
three  variables  0,  <J>,  and  M,  for  a  total  of  216  data  points  per  sector.  Using  the  sampling 
technique  of  Latin  Squares,  each  216-point  data  set  is  represented  by  a  sample  of  36 
points.  As  such,  a  total  of  256  data  points  will  be  experimentally  tested  in  all  seven 
sectors  for  a  complete  calibration  of  the  probe. 


D 

B 

C 

D 

E 

D 

B 

D 

D 

C 

D 

E 

C 

D 

E 

F 

B 

□ 

D 

A 

D 

E 

C 

B 

E 

C 

D 

B 

n 

D 

D 

E 

B 

D 

D 

C 

Figure  9.  6x6  Latin  Square  Plan 

Polynomial  Power  Series  Expansion  in  Three  Variables 

In  the  incompressible  calibration,  a  fourth  order  polynomial  expansion  in  two  vari¬ 
ables  was  used  for  a  total  of  15  terms  with  15  corresponding  calibration  coefficients. 

By  adding  a  third  variable  to  the  calibration,  the  number  of  terms  in  the  fourth  order 
expansion  jumps  from  15  to  35  with  35  corresponding  calibration  coefficients.  In  order 
to  obtain  a  valid  estimation  of  standard  deviation,  a  surplus  of  about  20  data  points 
over  the  number  of  calibration  coefficients  is  necessary  (Ref.  1).  This  sets  the  required 
number  of  data  points  for  a  given  sector  at  55;  however,  our  sample  size  has  already 
been  constrained  to  36  points  due  to  the  selection  of  a  6  x  6  Latin  Square.  It  is  possi¬ 
ble  to  take  more  than  one  36-point  Latin  Square  sampling  within  the  216-point  data  set, 
but  in  the  interest  of  keeping  the  total  number  of  points  required  for  a  complete  cali¬ 
bration  down  to  a  manageable  number,  so  that  time  spent  in  the  wind  tunnel  can  be  mini¬ 
mized,  we  elect  to  limit  each  of  the  seven  samples  to  36  data  points.  As  a  result,  a 
fourth  order  curve  fit  is  no  longer  feasible;  consequently,  the  polynomial  expansion  is 
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reduced  to  the  next  lowest  order. 

A  third  order  polynomial  expansion  in  three  variables  requires  20  calibration  coef¬ 
ficients.  This  leaves  us  with  a  16  data  point  surplus,  which  is  sufficient  to  calculate 
a  valid  standard  deviation.  Using  the  same  format  at  Eqn.  (9),  a  third  order  expansion 
in  three  variables  takes  on  the  following  form: 


A 


i 


Iff  + 


Order  of  Terms 
0th 


K2Cai  +  K^C6i  +  K4CMi  +  1st 

+  KtC^  +  KsCaiCBi  +  K9Ca.  CMi  +  KioC^C^  +  2nd 


KiiC^  +  KfzCj^  +  +  KisCj^  + 

K*6CaiC6i  +  K^C8iCMi  +  ^CaiC^  +  K^Cg^  +  K^oC^Cg^ 


where  A  is  either  8^,  0o,  or  Cq  for  the  inner  sector  and  6,  d> ,  CQ,  or  Cq  tor  the 
outer  sectors,  with  the  subscript  denoting  the  ith  such  quantity.  The  K's  are  the  cali¬ 
bration  coefficients,  where  the  superscript  denotes  the  quantity  to  which  a  particular 
set  of  K's  belong,  and  the  subscripts  identify  the  coefficient  of  a  particular  term  in 
Che  power  series  expansion.  In  the  case  of  an  outer  sector,  the  Ca's  and  Cg’s  are  re¬ 
placed  by  CQ1  s  and  C^'s,  respectively. 


Determination  of  Mach  Number 

As  stated  previously,  this  experiment  limits  its  scope  to  surveying  Mach  numbers 
slightly  below  sonic  flow  on  down  to  incompressible  flow.  Consequently,  there  are  no 
shocks  ahead  of  the  probe  and  the  isentropic  flow  relation  applies: 
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00 
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(21) 


Rearranging  Eqn.  (21)  and  deriving  Mach  number  in  terms  of  the  dynamic  to  total  pressure 
yields : 


(22) 


And  since  we  know  CD,  Cq,  and  the  seven  probe  pressures,  we  can  explicitly  determine 
the  local  Mach  number  from  the  dynamic  to  total  pressure  ratio.  For  the  inner  sector, 
this  ratio  is  determined  from  Eqns.  (4)  and  the  seven  probe  pressures: 
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(23) 


Analagous  equations  for  the  outer  sectors  are  developed  in  the  same  way,  except  that 
Eqns.  (8)  and  the  appropriate  probe  pressures  are  used  instead.  During  calibration, 
the  actual  values  of  Eqn.  (23)  and  its  outer  sector  counterparts  are  known  from  the  mea¬ 
sured  tunnel  conditions.  Solving  Eqn.  (21)  in  terms  of  the  dynamic  to  total  pressure 
ratio  as  a  function  of  Mach  number  gives  us: 
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(24) 


Thus,  if  the  total  and  static  pressures  cannot  be  measured  directly,  we  can  still  deter¬ 
mine  the  dynamic  to  total  pressure  ratio  if  only  the  Mach  number  is  known.  By  noting 
the  difference  between  Eqns.  (23)  and  (24),  we  can  estimate  how  accurately  our  polynomial 
curve  fit  determines  the  Mach  number  for  the  inner  sector;  similar  arguments  are  also 
extended  to  the  outer  sectors. 


Estimating  Accuracy  of  the  Curve  Fit 

Providing  there  are  approximately  twenty  more  data  points  than  the  number  of  calibra¬ 
tion  coefficients,  it  is  statistically  feasible  to  calculate  a  global  estimate  of  the 
accuracy  of  the  curve  fit  for  each  of  the  flow  parameters.  This  is  done  by  computing 
the  standard  deviation  of  the  difference  between  the  experimental  data  and  the  polynomial 
prediction  of  that  data.  For  flow  angles  and  Mach  number,  the  following  relation  applies 
(Ref.  7): 


a(A) 


■\/  n  £  (AEXPi  ‘  AP0LY1) 


(25) 


where  o  is  the  standard  deviation;  n,  the  total  number  of  data  points;  and  A,  the  desired 
flow  parameter. 

Even  though  the  total  and  dynamic  pressures  are  determined  from  C0  and.Cq,  the  ac¬ 
curacies  of  these  pressures  are  not  representative  of  the  accuracies  obtained  for  CQ 
and  Cq.  To  estimate  the  uncertainty  in  determining  the  total  and  dynamic  pressures  (for 
the  inner  sector,  for  example)  requires  the  defining  equations,  Eqns.  (17),  and  applying 
to  them  the  method  of  Kline  and  McClintock  (Ref.  8).  Then,  by  taking  the  standard  devia¬ 
tion  of  these  uncertainties  and  nondimensionalizing  them  with  the  dynamic  pressure,  we 
arrive  at  the  following  estimates  for  the  accuracies  of  these  pressures: 
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where  is  the  average  value  of  this  coefficient  for  a  given  sector  denoted  by  "n". 

Providing  the  errors  between  the  actual  values  and  the  polynomial  predictions  of 
those  values  are  normally  distributed,  there  is  a  68.3  percent  probability  that  the  poly¬ 
nomial  prediction  will  fall  within  one  standard  deviation  of  the  actual  value.  And  at 
1.96<?,  this  certainty  is  increased  to  95  percent. 


IV.  Apparatus  and  Calibration  Procedure 
Probe  Geometry 

The  seven-hole  probe  was  constructed  at  the  Air  Force  Academy  by  packing  seven 
properly-sized  stainless  steel  tubes  into  a  larger  stainless  steel  tube  as  shown  in  Figure 
10.  Although  the  inner  seven  tubes  have  an  outside  diameter  of  only  0.028  inches  with 
a  .005-inch  wall  thickness,  accurate  alignment  is  reasonably  insured  since  the  tubes 
can  only  pack  in  one  unique  way.  Once  assembled  in  this  arrangement,  the  tubes  are  then 
soldered  together  and  machined  to  provide  the  25-degree  half  angle  at  the  tip. 


Figure  10.  Probe  Geometry 


Probe  Mounting  Hardware 

The  test  facility  we  used  has  a  variable  position  sector  capable  of  traversing  25 
degrees  in  either  direction  for  a  total  sweep  of  50  degrees.  Yet,  for  a  complete  cali¬ 
bration  of  the  probe,  a  range  of  0  through  80  degrees  is  necessary.  Consequently,  we 
constructed  two  stings  (see  Figure  11):  a  15-degree  bent  sting  for  low  angle  measurements 
from  -10  to  40  degrees  angle  of  pitch,  and  a  55-degree  bent  sting  for  high  angle  measure¬ 
ments  from  30  to  80  degrees  angle  of  pitch.  In  addition,  both  stings  were  designed  to 
permit  the  tip  of  the  probe  to  pivot  about  a  fixed  point  in  the  center  of  the  tunnel. 
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Figure  11.  Sting  Geometry 


This  insured  uniform  flow  over  the  tip,  despite  changes  in  the  angle  of  pitch. 

Referring  to  Figure  12,  each  sting  has  36  holes  drilled  in  the  front  lace.  These 
holes  are  evenly  spaced  at  10-degree  intervals  and  circumscribe  a  complete  circle  in 
roll.  In  this  way,  roll  angle  is  accurately  set  by  engaging  the  alignment  pin  on  the 
probe  holder  with  the  alignment  hole  on  the  sting's  face. 


I 
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Data  Point  Selection 

Since  there  are  a  total  of  36  discrete  roll  angles,  we  are  allowed  to  test  six  dif¬ 
ferent  roll  angles  in  each  outer  sector.  This  is  the  primary  reason  why  we  chose  the 
6x6  Latin  Square.  Although  the  technique  of  Latin  Squares  allows  us  to  conveniently 
sample  large  data  sets,  there  are  some  drawbacks  associated  with  seven-hole  probe  cali¬ 
bration.  Primarily,  we  can  no  longer  test  the  entire  range  of  data  and  then  allocate 
the  data  to  a  given  sector  based  on  the  highest  port  pressure.  Instead,  we  must  deter¬ 
mine  beforehand  which  points  will  be  tested  and  the  sector  they  belong  to.  As  a  result, 
we  are  forced  to  draw  the  angular  cutoffs  for  each  sector  without  knowing  where  they 
actually  lie. 

In  terms  of  angle  of  pitch,  experience  with  incompressible  calibrations  has  shown 
that  30  degrees  is  the  smallest  angle  at  which  almost  all  points  still  fall  in  the  outer 
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Figure  12.  Probe,  Holder,  and  Sting  Assembly 

sectors.  Consequently,  we  selected  30  degrees  as  the  lower  limit  of  high  angle  measure¬ 
ment.  And  to  match  the  number  of  roll  measurements  with  an  equal  number  of  pitch  settings, 
we  measured  pitch  from  30  to  80  degrees  in  10-degree  increments  for  the  outer  sectors. 
Specifying  the  roll  angles  for  each  sector,  however,  presented  some  difficulties.  Assum¬ 
ing  a  perfectly  constructed  probe,  the  outer  sector  boundaries,  formed  by  the  isobars 
between  two  adjacent  peripheral  ports,  are  coincident  wxth  a  roll  angle  on  which  data 
is  taken.  This  predicament  is  illustrated  in  Figure  13a  for  sector  four.  In  this  exam¬ 
ple,  it  is  uncertain  whether  the  points  along  the  <)>  =  330  degrees  line  belong  in  sector 
four  or  five,  or  whether  the  points  on  the  <j>  =  30  degrees  line  belong  in  sector  three 
or  four.  To  resolve  this  problem,  we  rotated  the  probe  5  degrees  and  then  permanently 
fixed  it  to  the  holder  in  this  position.  From  Figure  13b  we  see  that  the  isobars  have 
shifted  off  the  points  and  that  each  sector  has  six  clearly  defined  roll  angles.  Realis¬ 
tically,  the  probe  is  not  constructed  perfectly,  causing  the  isobars  to  deviate  from 
their  ideal  Positions.  But  even  so,  the  isobars  can  deviate  up  to  5  degrees  in  either 
direction  before  the  points  of  one  sector  fall  in  another  sector.  At  this  point,  one 
might  ask  if  the  probe  can  be  accurately  calibrated  with  the  given  offset.  In  actuality, 
the  offset  makes  no  difference.  Since  only  the  probe  was  rotated,  leaving  the  index 
mark  of  the  probe  unchanged,  the  calibration  process  will  interpret  the  pressures  in 
terms  .if  the  roll  angle  set  by  the  index  mark.  An  alternative  to  the  offset  technique 
would  use  the  roll  angles  between  those  presently  measured  (i.e.,  3,  15,  25  .  .  .  355 
degrees);  however,  this  Is  an  option  we  ch.  not  to  take. 

Now  that  we've  determined  the  angular  coordinates  of  the  data  points  to  be  taken  in 
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Since  the  angular  coordinates  of  all  data  points  are  known,  the  last  task  is  to 
determine  the  Mach  number  at  which  each  is  sampled.  The  Mach  numbers  for  each  of  the 
six  letters  in  the  Latin  Square  of  Figure  9  are:  A,  M  =  0.37;  B,  M  =  0.45;  C,  M  =  0.53; 

D,  M  =  0.66;  E,  M  =  0.77;  F,  M  =  0.91.  The  resulting  distribution  of  data  points,  and 
the  Mach  number  each  was  tested  at,  can  be  inferred  from  Figure  14,  which  depicts  the 
inner  sector  and  a  typical  outer  sector. 

Test  Procedure  and  Software 

The  experiment  was  conducted  in  the  lxl  foot  blowdown  wind  tunnel  at  the  United 
States  Air  Force  Academy.  All  pressure  measurements  were  made  with  a  Scanivalve  Corpor¬ 
ation  Model  T  Scanivalve  and  pressure  transducer,  calibrated  by  a  Wallace  and  Tiernan 
sonar  manometer  digital  U-tube.  Data  were  collected  and  reduced  by  a  Digital  PDP  11/45 
computer  and  LPS-ll  Laboratory  Peripheral  System. 

Before  any  data  were  taken,  the  test  conditions  had  to  be  set  for  each  run.  Angle 
of  pitch  was  set  first  with  a  Gunner's  quadrant.  The  quadrant  was  clamped  to  the  sting's 
incidence  setting  block  (refer  to  Figure  12)  and  set  to  the  desired  angle  of  pitch  minus 
the  angle  of  the  elbow  on  the  sting  (either  15  or  50  degrees).  Next,  the  sector  was 
hydraulically  positioned  to  level  the  quadrant.  Since  the  sector  normally  has  a  tendency 
to  drift  between  runs,  we  froze  the  mechanism  by  driving  it  against  blocks.  Shims  were 
used  to  fine  tune  the  adjustments.  Using  this  technique,  the  pitch  was  set  to  within 
±0.0112  degrees  of  the  desired  angle,  and  without  drift.  Once  set,  the  pitch  was  not 
changed  until  all  points  at  that  angle  were  taken.  The  roll  angle  was  set  by  pulling 
the  spring-loaded  probe  forward  to  clear  the  alignment  pin,  then  rotating  it  to  the  de¬ 
sired  setting  and  engaging  it  into  the  new  position.  For  a  complete  360-degree  revolu¬ 
tion,  the  probe  was  first  rotated  counterclockwise  (as  viewed  from  the  front)  from  0 
to  180  degrees  in  10-degree  increments.  Next,  it  was  rotated  clockwise  to  the  190-degree 
position.  And  finally,  the  probe  was  rotated  counterclockwise  from  190  to  350  degrees 
in  10-degree  increments.  This  subjected  the  tubes  to  a  maximum  twist  of  180  degrees. 
After  setting  the  roll  angle,  the  tunnel  was  closed  by  a  hydraulic  ram  and  bolted  shut. 
Mach  number  was  then  set  by  manually  adjusting  the  tunnel's  inner  geometry  with  a  series 
of  hand  cranks  on  the  outside  of  the  tunnel. 

Prior  to  each  run,  the  above  variables  were  verified  and  input  into  the  computer 
by  a  program  titled  TRISHP.  Once  the  input  conditions  were  set,  air  was  blown,  and  after 
a  few  seconds  to  acquire  steady-state  conditions,  data  were  recorded.  Following  this, 
the  tunnel  was  reopened  and  the  process  repeated  for  a  new  data  point.  Typically,  each 
run  took  about  ten  minutes  to  set  up.  Consequently,  tests  ran  for  over  a  week  to  acquire 
all  256  data  points. 

After  all  data  were  taken,  a  second  program  titled  TRICAL  performed  the  matrix  opera¬ 
tions,  determined  the  calibration  coefficients,  and  estimated  the  accuracy  of  the  poly¬ 
nomial  expansions  in  fitting  the  known  data.  Having  done  this,  the  probe  calibration 
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is  complete. 

V.  Discussion 

The  primary  purpose  of  this  experiment  is  to  create  a  power  series  curve  fit  which 
accurately  determines  the  actual  flow  conditions  from  pressures  measured  on  the  probe. 

The  best  way  to  evaluate  the  performance  of  the  curve  fit  is  to  analyze  the  standard 
deviations  between  the  experimental  data  and  the  polynomial  determination  of  those  data. 
These  standard  deviations  are  presented  in  Table  1.  In  addition,  the  standard  deviations 
obtained  from  past  incompressible  calibrations  are  included  in  Table  2  for  comparison 
(Ref.  1). 

From  Table  1,  the  standard  deviations  of  and  8^  are  0.78  degrees  and  0.72  degrees 
respectively.  Although  these  values  appear  to  be  quite  good,  they  are  significantly 

Table  1 

STANDARD  DEVIATIONS 
COMPRESSIBLE  FLOW  CALIBRATION 


higher  than  those  obtained  from  the  incompressible  calibrations.  Specifically,  the  in¬ 
compressible  flow  calibrations  produced  standard  deviations  of  0.42  degrees  for  ot^  and 
0.36  degrees  for  B^,.  Similar  variations  are  also  noted  in  the  percent  errors  of  total 
and  dynamic  pressures.  Yet,  when  we  recall  that  the  incompressible  calibrations  used 
a  fourth  order  power  series  with  80  data  points  per  sector  and  the  compressible  calibra¬ 
tion  used  a  third  order  power  series  with  only  36  data  points  per  sector,  the  slight 
differences  are  quite  understandable.  Still  another  measure  of  curve-fit  accuracy 
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in  the  inner  sector  is  apparent  from  Figures  15  and  16.  Both  of  these  figures  depict 
isolines  of  Ca  and  Cg  plotted  against  uT  and  &T  with  CM  held  constant.  Figure  15  pre¬ 
sents  this  information  for  a  low  Mach  number,  with  Figure  16  representative  of  a  high 
Mach  number.  In  both  cases,  the  isolines  of  the  coefficients  are  nearly  straight  and 
relatively  orthogonal  to  each  other,  implying  a  linear  dependence  on  their  respective 
angle  and  independence  to  the  other  angle.  These  properties,  however,  begin  to  break 
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Figure  16.  Isolines  of  Ca  and  Cg  versus  aT  and  gT  for  Low  Angles  at  High  Mach  Number 

down  as  we  exceed  the  limits  of  our  experimental  data.  But  since  the  boundaries  of  the 
experimental  data  coincide  with  the  inner-outer  sector  interface  (see  Figure  17),  data 
points  lying  outside  the  inner  sector's  experimental  boundary  will  fall  into  the  outer 
sectors.  Consequently,  the  accuracy  in  determining  flow  properties  is  uniform  through¬ 
out  the  inner  sector. 

In  the  outer  sectors,  the  average  standard  deviation  of  the  error  in  calculating  if 
is  0.57  degrees  for  the  compressible  calibration  and  1.17  degrees  for  previous  incompres 
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Figure  17.  Ideal  and  Actual  Sector  Boundaries  Based  on 
the  Highest  Probe  Measured  Pressures 

sible  calibrations.  Judging  from  Figure  17,  one  might  expect  the  average  standard  devi¬ 
ation  of  <p  for  the  compressible  calibration  to  be  substantially  greater  than  that  for 
the  incompressible  calibrations,  because  the  actual  boundaries  of  sectors  one  through 
four  do  not  coincide  with  the  ideal  boundaries  for  which  the  calibration  was  made.  Al¬ 
though  this  phenomenon  is  entirely  consistent  with  the  manufacturing  anomalies  associated 
with  probe  construction,  the  actual  boundaries  cannot  be  determined  in  advance  and,  there¬ 
fore,  cannot  be  taken  into  account  in  a  calibration  scheme  using  the  method  of  Latin 
Squares.  Nevertheless,  despite  the  actual  locations  of  the  sector  boundaries,  the  stand¬ 
ard  deviation  in  <t>  for  the  compressible  calibration  agrees  favorably  with  its  incompres¬ 
sible  counterpart.  The  reason  for  this  is  shown  in  Figure  18.  Even  though  Figure  18 
is  based  on  data  taken  from  incompressible  calibrations,  it  is  representative  of  seven- 
hole  probes  in  general.  That  is,  the  coefficient  of  roll  continues  to  behave  linearly, 
or  ir.  a  way  that  can  be  represented  easily  by  a  polynomial,  slightly  beyond  the  isobaric 
sector  boundaries.  As  such,  the  effect  of  using  data  within  10  degrees  in  roll  from 
the  selected  boundaries,  as  Figure  18  suggests,  does  not  have  an  adverse  effect  on  fitting 
the  data.  Consequently,  the  disagreement  between  the  actual  and  ideal  sector  boundaries 
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Figure  18.  Incompressible  Data  Showing  C^,  versus  4>  in  Sector  1 
for  9  =  84  Degrees 


of  Figure  17  does  not  significantly  influence  the  determination  of  roll  angle. 

Despite  the  small  error  in  calculating  the  roll  angle,  the  standard  deviation  of 
the  error  in  calculating  pitch  angle  at  4.27  degrees  is  unacceptably  high,  especially 
when  compared  to  the  0.84-degree  standard  error  of  previous  incompressible  calibrations. 

The  actual  variation  between  the  experimental  data  and  its  polynomial  prediction  is  pre¬ 
sented  in  Figure  19,  which  depicts  a  typical  outer  sector.  From  this  illustration,  we  can 
see  that  the  greatest  error  in  determining  the  pitch  angle  occurs  at  high  angles  of  attack. 
This  occurs  as  a  result  of  the  polynomial's  inability  to  fit  the  actual  data.  Typically, 
the  shape  of  the  Cfi  versus  9  curve  looks  much  like  the  lift  curve  of  a  stalling  airfoil. 

Figure  20  illustrates  such  a  curve  for  sector  one,  based  on  the  extensive  data  available 

from  incompressible  probe  calibrations.  The  reason  why  the  curve  hooks  over  as  it  does 
is  evident  after  examining  the  two  pressures  comprising  the  numerator  of  the  coefficient 

of  pitch.  According  to  Figure  20,  the  center  port  pressure  decreases  with  increasing 

angle  of  pitch.  Beyond  some  point,  a  suction  develops  at  this  port,  causing  the  pressure 
there  to  dip  below  the  free  stream  static  pressure.  But  near  80  degrees  angle  of  pitch, 
the  suction  breaks  and  the  pressure  begins  to  increase.  As  this  occurs,  the  slope  of 
the  Cp7  curve  approaches  the  slope  of  the  Cpj  curve.  Once  the  two  slopes  are  equal , 
the  rate  of  change  of  the  numerator  is  zero,  causing  the  slope  of  the  coefficient  of 
pitch  also  to  be  zero.  No  calibration  may  be  made  beyond  this  point,  because  each  value 
of  C(5  then  corresponds  to  two  values  of  9.  And  since  9  is  a  function  of  Cfl ,  the  poly- 
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Figure  19.  Three-Dimensional  Data  Set  Showing  Correlation  Between 
Experimental  and  Polynomial  Data  for  Sector  5 
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nomial  can  onLy  calculate  one  unique  value  of  0  for  a  given  value  of  C<}.  Therefore, 
based  on  the  information  presented  in  Figure  20,  we  would  expect  to  be  capable  of  cali¬ 
brating  the  probe  out  to  80  degrees  angle  of  pitch.  This  was  indeed  the  case  for  the 
past  incompressible  calibrations,  but  not  the  case  in  our  compressible  calibration,  as 
Figure  19  vividly  points  out. 

The  failure  of  the  calibration  at  high  angles  of  pitch  is  most  likely  the  result 
of  probe  geometry.  Specifically,  the  half  angle  of  the  probe  used  in  our  compressible 
calibration  is  25  degrees  (see  Figure  10)  as  opposed  to  the  30-degree  half  angle  of  the 
probes  used  in  the  incompressible  calibrations.  This  steeper  half  angle  was  incorpor¬ 
ated  to  permit  a  closer  approximation  of  total  pressure  (as  measured  by  the  peripheral 
ports)  at  high  angles,  in  the  hope  of  extending  the  range  of  calibration  beyond  80  degrees 
of  pitch.  However,  this  reasoning  overlooked  the  effect  increasing  the  half  angle  would 
have  on  the  central  port.  That  is,  a  steeper  half  angle  requires  the  flow  passing  around 
the  probe  tip  to  turn  through  a  greater  angle,  causing  the  suction  over  the  central  port 
to  break  at  a  lower  angle  of  pitch.  Consequently,  the  slope  of  the  coefficient  of  pitch 
levels  off  earlier,  limiting  the  calibration  to  a  pitch  angle  below  80  degrees.  But 
even  as  this  theoretical  upper  limit  of  the  calibration  is  approached,  the  breakdown 
in  linearity  of  the  coefficient  of  pitch  with  increasing  angle  of  pitch  forces  the  poly¬ 
nomial  to  work  harder  to  fit  the  data.  Since  at  high  angles  of  attack  small  changes  in 
the  coefficient  of  pitch  result  in  large  changes  in  the  angle  of  pitch,  a  small  mismatch 
between  the  actual  data  and  the  polynomial  curve  fit  translates  into  a  large  error  as 
Figure  21  illustrates.  Thus,  as  a  result  of  the  difficulty  in  fitting  the  actual  data 
near  the  theoretical  limit  of  calibration  as  established  by  the  probe's  geometry,  the 
calibration  appears  to  be  limited  to  an  angle  of  pitch  below  80  degrees. 

To  confirm  these  suspicions,  we  ran  a  second  data  reduction,  excluding  all  80-degree 
pitch  data.  The  standard  deviations  of  the  truncated  data  set  are  displayed  in  Table  3 
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Figure  21.  limitations  of  Third  Order  Curve  Fit 
at  High  Angles  of  Attack 
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Table  3 

STANDARD  DEVIATIONS 
COMPRESSIBLE  FLOW  CALIBRATION 
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and  reflect  a  significant  decrease  in  the  error  associated  in  calculating  the  experimen¬ 
tal  quantities.  For  example,  the  standard  deviation  in  calculating  pitch  angle  went 
down  from  A. 27  degrees  to  0.79  degrees.  The  correlation  between  the  data  points  of  sec¬ 
tor  five  also  improved  dramatically  and  is  depicted  in  Figure  22.  However,  one  must 
keep  in  mind  that  since  the  80-degrec  points  were  removed  from  the  data  set,  the  data 
set  is  no  longer  as  accurately  represented.  That  is,  in  addition  to  losing  the  80-degre 
data,  we  also  lost  an  equal  number  of  associated  roll  angle  and  Mach  number  data.  Never 
theless,  the  greatly  reduced  standard  deviations  support  our  contention  that  for  our 
probe  in  the  vicinity  of  80  degrees,  the  third  order  curve  fit  is  incapable  of  accuratel 
fitting  the  data. 

VI .  Conclusions  and  Recommendations 

The  polynomial  expansion  in  three  variables  accurately  extends  the  calibration  of 
seven-hole  probes  into  the  compressible  regime.  Based  on  the  reasonably  close  correla¬ 
tions  between  the  standard  deviations  of  past  incompressible  calibrations  and  the  stand¬ 
ard  deviations  of  the  compressible  calibrations,  the  method  of  Latin  Squares  furnishes 
a  sample  space  which  accurately  represents  a  much  larger  three-dimensional  parameter 
space.  Finally,  the  third  order  curve  fit  accurately  represented  the  parameter  space 
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Figure  22.  Three-!) imensiona  1  Data  Set  Showing  Correlation  Between 
Experimental  and  Polynomial  Data  for  the  Truncated  Data  Set  of  Sector  5 
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out  to  an  angle  of  pitch  of  70  degrees,  but  fell  completely  apart  when  required  to  fit  the 
data  extending  to  80  degrees  angle  of  attack. 

It  is  recommended  that  an  additional  Latin  Square  be  run  in  an  outside  quadrant. 

A  set  of  calibration  coefficients  can  then  be  determined  using  that  Latin  Square  of  data 
only;  the  same  should  be  done  with  the  original  square.  The  results  can  then  be  examined 
to  see  how  well  the  calibration  coefficients  of  one  square  determine  the  values  of  the 
other  square  and  vice  versa.  Close  correlations  should  establish  beyond  a  reasonable 
doubt  the  validity  of  a  Latin  Square  in  representing  a  parameter  space.  Next,  both  Latin 
Squares  should  be  combined  to  determine  another  set  of  calibration  coefficients,  which 
in  turn  can  be  used  to  see  if  more  data  points  significantly  improve  the  degree  of  fit. 
Lastly,  using  the  combined  set  of  data  points,  extend  the  calibration  to  a  fourth  order 
polynomial  and  examine  it  to  see  if  the  addition  of  the  fourth  order  terms  significantly 
improve  the  degree  of  fit  up  to  80  degrees  angle  of  attack. 
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angle  between  probe  axis  and  velocity  vector  projected  on 
horizontal  plane  through  probe's  x-axis 

total  angle  between  velocity  vector  and  probe's  x-axis 

ratio  of  specific  heats 

angle  between  a  plane  containing  the  velocity  vector  and 
the  probe's  x-axis  and  a  vertical  plane  through  the  probe's 
x-axis  measured  positive  clockwise  from  port  number  four 
as  viewed  from  the  front 

standard  deviation 
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DEVELOPMENT  OF  DESIGN  CRITERIA  FOR  RAPID  PRE-EJECTION 
CREW  RESTRAINT  OF  UPPER  TORSO 

A.M.  Higgins*  and  G.  Kroh** 


Abstract 

This  paper  investigates  the  possibility  of  reducing  the  period  between  the  initia¬ 
tion  of  the  ejection  sequence  (i.e.,  lifting  ejection  initiation  handles,  pulling  the 
D-ring,  etc.)  and  the  first  motion  of  the  ejection  seat  out  of  the  aircraft.  The  pro¬ 
posed  method  of  reducing  this  time  period  is  to  provide  a  more  rapid  retraction  and 
restraint  of  the  crew  member's  upper  torso  during  the  ejection  process  than  is  presently 
used.  The  results  of  an  experimental  test  program  using  human  volunteers  are  presented 
in  this  report  and  demonstrate  that  human  subjects  can  easily  be  retracted  10  inches 
in  approximately  150  milliseconds.  Additional  tests  at  retraction  lengths  of  12  and 
14  inches  to  determine  if  similar  retraction  times  can  be  attained  for  these  retraction 
distances  are  recommended. 


I .  Introduction 

For  obvious  reasons,  an  aircraft  escape  system  must  function  as  rapidly  as  possible 
once  it  has  been  initiated.  This  aspect  of  escape  system  performance  has  become  increas¬ 
ingly  important  with  the  introduction  of  newer  aircraft  such  as  the  A-10  which  operate 
continuously  in  low-altitude,  high-speed  environments  where  milliseconds  lost  in  clearing 
the  aircraft  can  mean  the  difference  between  life  and  death  for  an  airman.  The  most 
obvious  point  in  the  ejection  sequence  where  time  can  be  reduced  is  during  the  period 
allowed  for  retraction  and  restraint  of  the  crew  member  prior  to  first  motion  of  the 
ejection  seat.  The  purpose  of  retracting  and  restraining  an  aircrew  member's  upper 
torso  is  to  align  the  crew  member's  spine  with  the  force  vector  of  the  ejection  catapult 
prior  to  actual  initiation  of  the  catapult  (Ref.  1).  It  should  be  noted  that  the  canopy 
of  the  aircraft  is  being  removed  to  clear  the  ejection  path  concurrently  with  this  re¬ 
traction.  The  present  Air  Force  escape  specification  (Ref.  2)  requires  that  the  powered 
inertia  reel  of  the  ejection  system  be  capable  of  positioning  an  aircraft  seat's  occupant 
through  0  to  18  inches  of  travel.  It  also  allows  a  retraction  time  of  0.3  seconds  at 
70  degrees  Farenheit.  This  300-millisecond  delay  requirement  is  significant,  and  its 
elimination  has  been  suggested.  In  fact,  the  Air  Force's  A-10  System  Program  Office 
at  Wright-Pat terson  AFB  at  one  time  actively  supported  the  removal  of  this  delay  from 
the  A-10's  ejection  seats.  What  is  really  being  suggested  is  a  trade-off  between  in¬ 
creased  injury  potential  (principally  to  the  spine)  that  would  exist  if  the  delay  for 
spinal  positioning  and  restraint  was  removed  and  decreased  mortality  due  to  a  more  rapid 
removal  of  the  aircrew  member  from  the  stricken  aircraft.  Unfortunately,  the  trade- 
ofL  is  extremely  difficult  to  quantify,  and  consequently  the  decision  to  retain  or  remove 
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the  time  delay  is  very  hard  to  make. 

If  one  looks  more  closely  at  this  time  delay  it  seems  to  have  been  based  on  require¬ 
ments  for  the  mechanical  inertia  reel  assembly  rather  than  consideration  of  human  phys¬ 
iological  tolerance.  The  question  a  designer  must  first  ask  then  is,  "How  fast  can 
a  crew  member  be  safely  retracted  and  restrained?"  If  it  can  be  done  safely  in  150 
milliseconds,  for  example,  we  have  gained  valuable  time  without  compromising  the  neces¬ 
sity  of  spinal  alignment  before  the  catapult  force  of  the  ejection  system  in  applied 
to  the  human  body.  In  fact,  a  requirement  for  an  inertia  reel  that  can  accomplish  this 
goal  was  recently  published  as  an  Air  Force  Technical  Need  (Ref.  3).  For  this  analysis 
we  assumed  that  the  mechanical  requirements  of  a  more  rapid  retraction  can  accomplished 
as  well  as  faster  removal  of  the  aircraft  canopy.  We  therefore  set  out  to  determine 
if  the  retraction  time  can  be  safely  reduced  to  something  significantly  less  than  300 
milliseconds. 

II .  Test  Approach 

We  felt  that  the  only  way  we  could  conclusively  demonstrate  actual  reduced  retrac¬ 
tion  time  was  to  physically  experiment  and  measure  the  retraction  of  human  bodies.  Of 
course,  one  cannot  "just  begin"  when  human  volunteer  subjects  are  involved.  The  testing 
had  to  be  approached  carefully  with  an  engineering  test  plan  and  a  medical  human-use 
protocol  that  was  carefully  thought  out  and  which  considered  the  safety  of  the  human 
test  subjects  as  paramount. 

Our  approach  was  to  first  analyze  each  part  of  the  retraction  and  restraint  sequence 
and  look  for  critical  parameters  that  characterize  the  test  severity.  For  example, 
the  initial  stage  of  retraction  involves  large  retraction  cable  forces  which  are  trans¬ 
mitted  to  the  shoulders  of  the  subject  through  the  restraint  harness.  One  critical 
parameter  for  this  phase  of  retraction  then  was  cable  force.  At  the  same  time,  a  care¬ 
ful  review  of  existing  literature  had  to  be  done  to  determine  the  tolerable  limits  for 
these  critical  parameters.  Here,  for  example,  we  had  to  determine  what  levels  of  cable 
force  (harness  loads)  are  allowable  and  what  levels  were  expected  in  the  tests. 

Unfortunately,  little  information  on  the  retraction  of  humans  in  ejection  seats 
was  available  (Ref.  I  and  4).  Therefore,  the  following  test  procedure  approach  was 
adopted  in  order  to  gather  information:  tests  would  begin  at  a  level  of  retraction  speed 
that  would  not  jeopardize  the  human  subjects,  and  they  would  then  progress  in  steps 
of  severity,  where  the  test  severity  was  judged  to  be  a  function  of  retraction  length 
as  well  as  the  speed  of  retraction.  Consequently,  we  decided  to  do  three  test  series. 

The  first  series  would  involve  only  10-inch  retractions  and  would  begin  at  a  retraction 
speed  where  previous  test  data  was  available  and  the  retractions  were  shown  to  be  toler¬ 
able.  The  second  and  third  series  would  involve  12-inch  and  14-inch  retractions  beginning 
at  the  initial  retraction  speeds  used  in  the  10-inch  tests. 

The  reason  we  did  not  plan  to  conduct  retraction  tests  for  retraction  lengths  greater 
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Chan  14  inches  is  that  although  operational  inertia  reels  generally  provide  18  inches 
of  strap  length,  the  actual  maximum  retraction  length  is  approximately  14  to  15  inches 
for  a  crew  member  who  is  grossly  out  of  position  in  his  seat.  This  situation  would 
most  likely  occur  in  a  multi-seat  aircraft  in  which  a  crew  member  is  unexpectedly  ejected 
from  an  aircraft  undergoing  high  "G"  maneuvers  or  in  an  accidental  ejection.  Each  test 
series  would  begin  at  a  harmless  level  of  speed  and  progress  in  increasing  increments 
of  severity. 

To  establish  the  expected  maximum  levels  of  the  selected  critical  retraction  param¬ 
eters  we  ran  a  pre-test  series  using  a  95th  percentile  anthropometric  dummy  and  recorded 
these  critical  retraction  parameters  for  each  retraction  phase.  We  then  took  the  abso¬ 
lute  maximums  of  these  critical  parameters  and  showed  that  they  were  below  human  toler¬ 
ance  levels  (Ref.  5  and  6).  During  the  human  test  program,  we  recorded  these  same  criti¬ 
cal  parameters  for  each  human  subject  and  compared  them  to  the  values  obtained  from 
the  tests  using  the  dummy.  With  this  technique  we  could  compare  the  trends  in  data 
as  we  progressed  to  higher  levels  of  test  severity.  This  technique  allowed  us  to  see 
if  the  values  obtained  from  the  human  tests  were  larger  than  those  obtained  from  the 
tests  using  the  dummy  or  if  the  trends  were  such  that  we  could  predict  more  severe  read¬ 
ings  for  the  humans  at  the  next  level  and  thus  stop  the  tests.  We  required  that  all 
subjects  be  tested  at  each  severity  level  before  we  advanced  to  the  next  level  in  order 
to  insure  all  inputs  were  available  before  we  made  a  decision  to  go  to  the  next  test 
level. 

III.  Experimental  Apparatus  and  Procedure 

The  primary  piece  of  test  equipment  was  the  Body  Positioning  and  Restraint  Device, 
BPRD  (Ref.  1),  located  in  the  Protection  Branch  of  the  Aerospace  Medical  Research  Labora¬ 
tory  (AMRL)  at  Wright-Patterson  Air  Force  Base.  This  device  is  essentially  a  mock-up 
of  an  ejection  seat  that  has  a  number  of  electronically-actuated,  hydraulically-driven 
pistons  and  a  control  system  housed  in  the  area  behind  the  seat.  Each  piston  has  a 
cable  attached  to  it,  and  through  pulley  arrangements  these  cables  are  routed  to  various 
points  where  they  can  be  connected  to  the  test  subject's  body  by  a  harness.  Only  the 
upper  torso  piston  and  cable  were  used  in  these  tests,  and  Figure  1  shows  a  schematic 
of  the  BPRD  test  set-up.  Figure  2  shows  a  subject  on  the  BPRD  in  three  positions  during 
a  typical  retraction  experiment.  The  positions  shown  are  full  forward,  an  intermediate 
position,  and  full  back.  Data  sensor  locations  are  shown  and  these  sensors  will  be 
discussed  in  detail  later.  The  harness-cable  arrangement  can  be  seen  passing  over  the 
subject's  shoulders  and  through  the  seat  back. 

To  retract  the  subject  a  switch  is  thrown  which  actuates  a  solenoid,  allowing  the 
piston  to  be  hydraulically  driven  and  thus  retracting  the  cable  and  the  subject.  Higher 
hydraulic  pressures  impart  higher  forces  to  the  piston  and  result  in  faster  retractions. 
Therefore,  in  these  tests  hydraulic  pressure  becomes  a  measure  of  retraction  severity. 


96 


I  SA!  A- IK-SI-A 


Figure  1.  Rodv  Positioning  and  Restraint  Device  (BPRD)  Test  Set-ip 


ieure  2.  Subject  on  BPRD  in  Three  Positions  Durin 
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The  test  procedure  used  in  all  tests  is  clearly  spelled  out  by  the  BPRD  Test  Sequence 
Checklist  which  appears  in  Figure  3.  This  checklist  was  approved  by  the  AMRL  safety 
officer  prior  to  testing  and  was  followed  exactly  for  all  human  tests.  It  essentially 
assures  that  (1)  the  sutject  is  appropriately  dressed  and  physically  ready  for  the  test, 
(2)  the  test  instrumentation  and  equipment  are  operating  properly,  (3)  all  initial  mea¬ 
surements  have  been  taken,  (4)  the  subject  is  properly  positioned  for  the  test,  and 
(5)  all  test  steps  have  been  taken  in  the  proper  order.  After  the  test,  it  assures  that 


BPRD 

TEST  SEQUENCE 

TEST  CONDUCTOR'S 

CHECKLIST 

1 . 

Subject  No. 

19. 

Piston  Position 

2. 

Test  No. 

20. 

Blood  Pressure 

3. 

Subject  Weight 

21 . 

Mouthpiece 

4. 

Check  ECG  Signal 

22. 

Helmet  Strap  Fastened 

5. 

Cameras  Loaded 

23. 

Accel  Cal/Zero 

a.  Side  Cam  No. 

24. 

Subject  Position  (Preload)* 

b.  Lense  No. 

25. 

Cable  Clearance 

c.  Leg  Cam  No. 

26. 

Body  Measurement 

d.  Lense  No. 

27. 

Pre-Test  Photo 

e.  Front  Cam  No. 

28. 

Bldg  40 

f.  Lense  No. 

29. 

Pressure  Check 

g.  Frame  Rate 

30. 

Piston  Length 

h.  f-Stop 

31. 

Safety 

6. 

Subject  Arrival* 

32. 

Clearance  to  Operator 

7. 

Blood  Pressure 

33. 

Photo  Ready 

8. 

Chest  Pack 

34. 

Inst  Ready 

9. 

Harness  Fit 

35. 

Med  Monitor  Ready 

10. 

Brief  Subject 

36. 

Subject  Ready 

11. 

Helmet  Fit 

37. 

Test* 

12. 

Subject  Seated 

38. 

Post  Test  Med 

13. 

Fiducial  Placement 

39. 

Post  Test  Photo 

14. 

Pistons  Retracted 

40. 

Date 

15. 

Harness  Adjusted 

41. 

Remarks 

16. 

Subject  Position 

17. 

Harness  Attached 

18. 

Lap  Belt  Attached 

*Record  Time 

Figure  3.  BPRD  Test  Sequence  Checklist 


the  subject's  condition  is  checked  and  all  post-test  data  are  recorded. 

As  part  of  this  procedure,  the  subject  was  asked  to  lean  into  the  harness  and  apply  a 
10-to  20-pound  preload  to  the  cable  to  simulate  the  normal  cable  tension  exerted  by 
the  inertia  reel  on  an  aircrew  member  in  the  cockpit.  The  actual  value  of-  the  tension 
was  displayed  to  the  subject  on  a  digital  meter  so  he  could  adjust  it.  It  should  also 
be  pointed  out  that  we  purposely  did  not  reproduce  one  aspect  of  operational  retractions 
with  our  experiments.  In  our  tests,  the  subject  was  fastened  to  the  retraction  cable 
with  the  piston  fully  retracted,  and  the  subject  then  adjusted  his  harness  and  seat 
belt  for  comfort.  This  harness  adjustment  technique  then  assures  that  no  significant 
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afterload  would  be  placed  on  the  subject  by  the  hydraulics  after  the  piston  is  retracted. 
In  the  cockpit  the  inertia  reel  usually  does  produce  an  afterload.  It  should  also 
be  noted  that  a  fully  instrumented  retraction  test  with  a  dummy  was  run  as  a  check  of 
the  BPRD  prior  to  the  first  human  test  of  the  day.  Also,  as  a  safety  precaution,  high¬ 
speed  motion  pictures  from  at  least  one  test  level  were  carefully  reviewed  prior  to 
exposing  subjects  to  the  next  test  condition. 

IV.  Data 

The  data  obtained  from  these  tests  can  be  divided  into  three  separate  categories: 
medical,  electronic,  and  photographic.  The  specific  measurements  made  in  each  of  these 
categories  are  listed  below, 

A.  Medical  Data 

1.  Pre-test  physical  examination  of  subject 

2.  EKG  and  pulse  before,  during,  and  after  the  test 

3.  Blood  pressure  before  and  after  the  test 

4.  Post-test  physical  examination  of  subject 

5.  Subjective  assessment  of  retraction  severity  by  subject 

B.  Electronic  Data 

1.  Chest  acceleration  in  all  three  axes,  resultant  and  chest  severity  index 

2.  Head  acceleration  in  all  three  axes  and  resultant 

3.  Retraction  cable  force 

4.  Piston  displacement 

5.  Piston  velocity 

6.  Piston  acceleration 

7.  Back  impact  force:  three  load  cells 

8.  Head  impact  force:  three  load  cells 

9.  Pelvis  impact  force:  two  load  cells 

10.  Seat-pan  force:  three  load  cells 

C.  Photographic  Data 

1.  Photometric  data 

2.  Video  recordings 

3.  Photographic  stills 

In  this  paper  the  only  medical  data  we  will  comment  on  is  the  subjective  assessment 
of  the  retraction  severity  made  by  the  subjects. 

We  should  note  that  summing  of  the  load  cell  outputs  attached  to  a  particular  surface 
will  yield  the  total  impact  force  on  that  surface,  and  cable  force  (harness  load)  is 
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measured  in  the  retraction  cable  by  means  of  a  load  cell  positioned  under  a  pulley  over 
which  the  cable  is  routed.  The  chest  acceleration  is  measured  by  a  triaxial  accelerometer 
pack  mounted  on  the  xiphoid  area  using  a  belt  arrangement.  Head  acceleration  was  meas¬ 
ured  using  a  newly  developed  accelerometer  mouth  mount.  The  triaxial  acceleration  pack¬ 
age  was  attached  to  a  metal  mouthpiece  as  shown  in  Figure  4  (Ref.  7).  Each  of  these 
mouthpieces  were  individually  cast  by  the  dental  clinic  at  Wright-Patterson  AFB  to  fit 
each  subject's  tooth  pattern.  The  weight  of  the  complete  assembly  (mouthpiece  and  ac¬ 
celerometer  pack)  is  50  grams. 


A.  Mouthpiece  Attached  to  Subject's  Teeth 


B.  Mounting  Block  Attachment  to  Mouthpiece  C.  Assembled  Accelerometer  Mouth  Mount 

(Accelerometers  and  plastic  mounting 
block  are  covered  with  an  electri¬ 
cally  insulating  material) 


Figure  4.  Accelerometer  Mouth  Mount 
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A  high-speed  photographic  film  (nominal  500  frames  per  second)  was  used  to  determine 
body  kinematics  by  individually  tracking  various  points  on  the  subject's  body.  This 
photometric  data  will  not  be  discussed  here.  The  video  recordings  (120  frames  per  second) 
were  obtained  with  a  television  camera  and  were  used  by  the  principal  investigator  and- 
medical  monitor  to  review  the  body  kinematics  of  the  subject  immediately  after  a  test 
and  prior  to  the  next  subject's  exposure. 

V.  Test  Subject  Preparation 

All  of  the  subjects  were  volunteers  from  the  AMRL  Hazardous  Duty  Panel  and  demon¬ 
strated  normal  skull,  spine,  and  chest  x-rays,  EEG,  pulmonary  function,  treadmill  exer¬ 
cise  testing,  and  general  physical  examination.  Use  of  these  volunteers  was  in  strict 
accordance  with  AFR  80-33  and  its  supplement. 

The  subjects  wore  leotards  (which  were  dyed  orange  to  enhance  photographic  coverage) , 
AF  flight  helmets  (HGU-26/P) ,  and  integrated  parachute  harnesses  (PCU-15/P  for  males 
and  PCU-16/P  for  females).  The  subjects  donned  their  helmets  and  harnesses  prior  to 
mounting  the  BPRD. 

Also,  prior  to  testing,  the  medical  investigator  precisely  identified  various  points 
on  the  subject's  body  for  photographic  tracking.  A  tape  fiducial  was  then  positioned 
over  the  spot.  This  procedure  assured  consistent  fiducial  location  from  subject  to 
subject,  and  for  the  same  subject  from  test  to  test. 

VI.  Experimental  Tests 

The  tests  described  here  were  completed  in  two  phases.  Phase  1  was  a  pilot  study 
done  at  low  hydraulic  retraction  pressures  (200,  300  PSI)  to  validate  measurement  tech¬ 
niques,  evaluate  potential  human/test  apparatus  measurement  techniques,  check  out  data 
programming  procedures,  and  to  serve  as  an  orientation  experience  for  the  volunteer 
subjects.  These  tests  essentially  repeated  earlier  tests  performed  at  AMRL  (Ref.  8). 

Phase  II  extended  the  studies  from  the  harmless  levels  of  the  Phase  I  tests  to  retraction 
speeds  comparable  to  those  experienced  operationally.  Maximum  retraction  pressure  reached 
in  Phase  II  was  900  PSI.  The  retraction  distance  was  10  inches  in  all  tests. 

VII.  Discussion  of  Results 

Eleven  human  subjects,  nine  males  and  two  females,  participated  in  these  tests, 
and  they  represent  a  wide  range  of  sitting  heights  and  body  weights  as  shown  in  Table  1. 

Each  of  these  subjects  were  tested  at  each  severity  level  (retraction  pressure 
level),  and  Figure  5  is  a  series  of  photographs  showing  a  typical  retraction  sequence. 
These  photographs  were  taken  from  selected  frames  of  the  high  speed  film  and  display 
both  the  side  view  and  the  front  view  of  the  retraction  event.  These  sequential  photo¬ 
graphs  show  subject  F  being  retracted  at  the  maximum  severity  level  (900  PSI),  and  the 
times  shown  with  the  photographs  correspond  to  an  important  event  in  the  retraction 
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Tab  1 e  1 

SUBJECT  SITTING  HEIGHT  and  weight  range 


MID-SHOULDER  SITTING 
HEIGHT  RANGE 

WEIGHT  RANGE 

59.3  cm  to  68.0  cm 

118  ibs  to  210  lbs 

5  percentile  60.2  cm 

95  percentile  69.2  cm 

_ 

5  percentile  140  lbs 

95  percentile  210  lbs 

sequence.  For  example,  Figures  5c  and  5f  show  the  subject's  position  at  the  point 
where  the  retraction  piston  has  just  been  fully  retracted.  This  retraction  time  can 
easily  be  determined  from  the  displacement  transducer  data  shown  in  Figure  6.  The  piston 
retraction  time,  then,  is  the  time  elapsed  from  the  beginning  of  piston  motion  to  the 
point  where  the  piston  is  full  down  (designated  as  Point  3).  Similarly,  the  other  event 

times  can  be  determined  from  the  load  cell  output,  as  shown  in  Figure  7.  It  is  inter¬ 

esting  to  note  here  that  the  peak  cable  force  coincides  with  the  peak  seat-pan  load, 
which  indicates  that  the  subject  is  not  only  pulled  back  into  the  seat,  but  is  also 
pulled  down  into  the  seat.  Notice  that  Point  4  is  selected  to  be  the  point  where  seat 
back  load  reaches  its  peak. 

The  chest  acceleration  for  this  test  is  shown  in  Figure  8.  This  figure  shows  that 

the  maximum  chest  acceleration  for  this  test  is  12  G's  and  occurs  during  the  onset  of 

cable  loading  when  the  chest  is  being  accelerated  and  again  when  the  torso  impacts  the 
seat  back  and  is  being  brought  to  rest. 

The  subject's  head  acceleration  for  this  test  is  shown  in  Figure  9.  Here  the  maxi¬ 
mum  resultant  acceleration  is  11.5  G's.  This  is  the  acceleration  value  that  was  used 
to  compare  head  acceleration  severity  between  the  various  subjects  for  the  composite 
results  as  discussed  below. 

Individual  results  such  as  those  just  discussed  were  recorded  for  each  subject 
at  every  test  level  and  were  then  plotted  in  comparison  to  the  test  data  obtained  using 
the  dummy  subject.  These  composite  results  for  piston  retraction  time,  for  example, 
are  shown  in  Figure  10,  where  the  abscissa  represents  the  hydraulic  pressure . level  or 
the  retraction  severity  level.  The  data  indicates  clearly  that  increasing  the  pressure, 
i.e.,  the  force  applied  to  the  piston,  results  in  a  faster  retraction,  as  we  expected. 

The  body  retraction  time  was  defined  as  the  time  from  the  onset  of  cable  force  to  the 
peak  of  seat  back  force.  This  appeared  from  our  tests  to  be  a  reasonable  criterion 
for  completion  of  retraction.  Cessation  of  head  motion  or  head  orientation  was  not 
used  as  a  measure  of  completion  since  some  subjects  had  continued  head  motion  past  this 
point,  ending  with  upright  head  posture,  while  others  maintained  a  head-forward  position. 

In  either  case,  head  motion  or  position  would  seem  to  be  an  inappropriate  criterion  for 

i 


Figure  8.  Chest  Acceleration  Recorded  for  Figure  9.  Head  Acceleration  Recorded  for 

Subject  F  at  Highest  Severity  bevel  (900  1’SI)  Subject  F  at  Highest  Severity  I.evel  (900  PSI) 
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Figure  10.  Composite  Piston  Retraction  Times 


measuring  sequence  completion.  For  example,  if  we  consider  the  retraction  situation 
where  the  subject's  torso  is  already  on  the  seat  back,  the  retraction  would  certainly 
end  when  the  harness  is  tightened  (the  peak  back  force  point),  and  head  motion  would  not 
be  a  consideration. 

Figure  10  also  shows  that  at  the  higher  pressure  levels  we  have  reached  a  point 
of  diminishing  returns.  Notice  also  that  the  human  is  always  retracted  more  rapidly 
than  the  dummy,  although  at  the  higher  pressures  the  difference  becomes  insignificant. 

The  female  subjects  had  lower  body  mass  and  were,  therefore,  retracted  more  rapidly. 
Their  piston  retraction  time  was  always  at  the  lower  end  of  the  range  shown  in  Figure 

10,  whereas  the  largest  male  was  always  near  the  top  of  the  retraction  time  range. 

The  maximum  cable  force  obtained  for  each  subject  at  each  level  is  shown  in  Figure 

11.  Here,  as  expected,  cable  force  increased  with  pressure  level.  Note  that  the  human 
peaks  are  always  less  than  the  data  for  the  dummy.  For  the  human  subjects,  the  large 
males  always  experienced  the  largest  cable  forces  and  the  females  the  smallest.  This 

is  consistent  with  results  of  the  piston  retraction  behavior  since  a  larger  mass  requires 
greater  force  to  accelerate  it. 

The  maximum  seat-pan  force  obtained  for  each  subject  and  the  dummy  are  shown  in 
Figure  12.  Three  load  cells  were  used  to  measure  the  seat-pan  load  and  the  values  plotted 
in  Figure  12  represent  the  sum  of  those  load  cells.  The  load  cells  were  set  to  zero 
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Figure  11.  Maximum  Cable  Force  at  Each  Pressure  Level 
for  All  Human  Subjects  and  Dummy 
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Figure  12.  Seat-Pan  Load  for  Each  Pressure  Level 
for  All  Human  Subjects  and  Dummy 
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prior  to  the  retractions  to  remove  the  effect  of  body  weight  so  that  the  values  obtained 
would  reflect  only  the  down-loading  on  the  spine  due  to  the  retraction  process.  Figure 
12  clearly  demonstrates  that  this  loading  does  occur  and  reaches  a  maximum  value  of 
approximately  130  pounds.  This  value,  however,  is  below  the  values  for  the  dummy  and 
well  within  human  tolerance. 

The  composite  chest  acceleration  results  are  shown  in  Figure  13  and  demonstrate 
that  the  humans  experienced  significantly  lower  values  of  maximum  acceleration  than 
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Figure  13.  Maximum  Chest  Accelerations  at  Each  Pressure  Level 
for  All  Human  Subjects  and  Dummy 


the  dummy.  This  data  also  demonstrates  the  nonlinear  trend  of  the  dummy's  chest  acceler¬ 
ation  which  we  were  afraid  we  might  encounter  in  the  tests  of  humans.  Our  concern, 
of  course,  was  that  although  all  test  results  from  one  test  level  would  be  well  within 
tolerance  limits,  a  large  increase  in  test  severity  could  occur  between  test  levels 
because  of  the  nonlinear  behavior  of  the  biological  system,  and  we  would  therefore  expose 
the  test  subject  to  possible  injury.  The  human  data,  however,  is  approximately  linear 
and  demonstrates  a  fairly  small  range  in  resultant  acceleration  peaks  at  each  level 
of  hydraulic  pressure. 

The  composite  head  acceleration  results  shown  in  Figure  14  demonstrate  a  similar 
trend.  Here  again,  the  dummy  experienced  much  higher  levels  of  acceleration  than  the 
human  subjects.  This  was  primarily  due  to  the  dummy's  head  impacting  the  headrest  of 
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the  seat.  These  impacts  began  at  retraction  pressure  levels  of  as  low  as  400  PSI,  and 
at  the  higher  levels  of  retraction  pressure  these  impacts  were  quite  severe.  In  fact, 
some  of  these  impacts  were  so  severe  that  helmet  fractures  occurred.  Our  risk  analysis 
indicated  that  the  active  muscular  effects  of  the  human  neck  would  significantly  reduce 
the  head  motion  and  thereby  provide  safe  exposure.  Our  human  tests  at  low  pressure 
levels  confirmed  that  the  acceleration  levels  of  the  dummy  head  were  significantly  higher 
than  those  of  the  human  volunteers. 

We  therefore  proceeded  with  our  tests  using  the  stepped  retraction  procedure  pre¬ 
viously  discussed.  To  our  surprise,  we  observed  no  head  impacts  of  the  humans  at  all, 
even  at  the  highest  pressure  levels.  As  Figure  14  indicates,  the  maximum  resultant 
human  head  acceleration  was  less  than  one-half  of  the  dummy's  at  the  maximum  severity 
level.  This  result  once  again  demonstrated  the  importance  of  continued  human  testing. 
Data  obtained  from  these  tests  can  be  used  to  complete  optimal  system  designs  that  are 
demonstrably  safe,  whereas  data  from  tests  utilizing  dummies,  human  cadavers,  or  animals 
would  indicate  injuries  at  these  levels  and  lead  to  system  designs  that  are  not  only 
overweight  and  difficult,  but  which  also  degrade  system  performance. 


VIII.  Summary  and  Conclusions 


The  results  of  these  tests  indicate  that  retraction  of  the  human  torso  in  a  one-G 


field  can  be  safely  performed  in  significantly  less  than  300  milliseconds.  The  means 
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and  standard  deviations  of  the  various  retraction  test  parameters  are  shown  in  Table  2. 
Remember  that  the  retraction  time  was  defined  as  the  elapsed  time  between  the  onset 
of  cable  force  to  the  point  where  maximum  seat  back  force  occurred. 


Table  2 

RETRACTION  TEST  PARAMETERS 
(Tabulated  values  are  experimental  means  with 
the  standard  deviation  in  parentheses) 


TEST  PARAMETERS 

RETRACTION  LENGTH 

10  INCHES 

600  PSI 

900  PSI 

Retraction  Time 
(milliseconds) 

148  (8.2) 

134  (11.5) 

Peak  Cable  Force  (lbf) 

185  (24.3) 

250  (28.5) 

Average  Reel  Retraction 

Velocity  (in/sec) 

77.9  (2.4) 

97.4  (2.7) 

Peak  Seat-Pan  Force  (lbf) 

75  (18) 

102  (23.4) 

Peak  Back  Force  (lbf) 

354  (71) 

466  (70.6) 

Peak  Chest  Acceleration  (C's) 
Retraction 

Seat  Back  Impact 

7.2  (1.4) 

7.5  (1.2) 

11.6  (1.5) 
11.9  (1.9) 

Peak  Head  Acceleration  (G's) 
Retraction 

Deceleration 

7.1  (1.8) 

5.9  (2.3) 

11.7  (3.0) 
10.2  (3.3) 

All  of  the  measured  parameters  are  well  within  human  tolerance  limits.  The  subjects 
tolerated  the  tests  well  at  all  severity  levels,  and  considered  the  tests  to  be  harmless 
even  at  the  highest  severity  level. 

We  conclude  then,  that  for  for  10-inch  retraction  distances,  times  much  less  than 
300  milliseconds  are  required  for  upper  torso  retraction. 

The  next  step  in  this  testing  process  will  be  to  complete  a  test  series  of  upper 
torso  retractions  at  longer  distances  of  up  to  14  inches  to  determine  if  equivalent 
retraction  time  reductions  can  be  obtained  for  retraction  distances  equalling  the  maxi- 
mums  encountered  in  operational  retractions.  If  so,  then  it  seems  possible  that  a  sig¬ 
nificant  improvement  ir.  crew  survival  can  be  realized  through  a  reduction  in  the  initial 
time  delay  that  occurs  between  ejection  system  initiation  and  actual  seat  motion. 
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SI  UNITS  (MASS’,  WEIGHT?) 

LET'S  DO  IT  RIGHT 

R.F.  Felton* 

A  few  years  ago  the  United  States  began  in  ernest  to  change  from  the  English  system 
of  units  to  the  new  International  System  of  Units  (SI  units).  Those  who  encouraged  the 
change  had  high  hopes  that  some  of  our  bad  habits  in  using  the  English  system  would  dis¬ 
appear  in  the  change,  but  alas,  it  seems  that  the  American  public  still  does  not  under¬ 
stand  the  difference  between  mass  and  weight.  This  paper  examines  the  concepts  of  mass 
and  weight  by  using  force  to  explain  the  origin  of  the  difficulty  in  understanding  the 
correct  relationship  between  mass  and  weight. 

I.  Introduction 

In  1960  the  United  States,  along  with  36  other  countries,  took  part  in  officially 
establishing  the  "Systeme  International  d'Unites"  or  the  International  System  of  Units 

(SI  units).  Although  the  United  States  government  intended  to  rapidly  switch  from  use 

of  the  English  system  of  pounds,  feet,  and  seconds  to  the  new  SI  units  of  kilograms, 
meters,  and  seconds,  the  change  has  been  slow  in  coming.  It  is  only  in  recent  years 
Chat  the  change  has  really  made  itself  felt.  Unfortunately,  some  of  our  bad  habits  with 
regard  to  our  concepts  of  mass  and  weight,  which  we  had  in  the  English  system  of  units, 

have  found  their  way  into  the  new  world  of  SI  units.  In  fact,  the  problem  is  even  more 

severe  in  SI,  since  people  do  not  yet  have  a  feel  for  the  new  units. 

This  confusion  between  the  concepts  of  mass  and  weight  results  from  the  fact  that 
most  people,  including  many  engineers  and  scientists,  do  not  really  understand  the  English 
system  of  units.  And  it  is  no  wonder,  since  there  is  more  than  one  English  system  of 
units.  Units  (i.e.,  man's  attempt  to  put  weights  and  measures  on  some  relative  scale) 
were  first  developed  by  people  in  order  to  communicate  and  conduct  their  business  with 
others.  There  are  only  a  few  fundamental  units  (sometimes  referred  to  as  standard  or 
base  units).  These  fundamental  units  are  physical  creations  of  man  upon  which  all  mea¬ 
surements  depend.  For  example,  the  unit  for  the  dimension  of  length  is  a  fundamental 
unit.  All  the  other  units  are  derived  from  fundamental  units  through  the  use  of  equations 
(relationships  that  explain  observed  phenomena).  For  example,  the  unit  for  velocity 
is  a  derived  unit  (velocity  =  length/time;  both  length  and  time  are  fundamental  dimensions 
and  as  such  have  defined  fundamental  units,  and  velocity  is  derived  from  this  relation¬ 
ship).  A  number  of  different  unit  systems  have  been  developed  over  the  years  (Ref.  1). 

In  all  these  systems,  it  was  necessary  for  someone  to  define  the  few  fundamental  units. 

The  actual  definitions  of  these  fundamental  units  are  not  crucial  to  the  problem  of  under¬ 
standing  the  difference  between  mass  and  weight;  the  important  thing  is  that  there  are 
units  that  are  fundamental  and  are  defined  to  be  a  certain  thing  by  agreement  between 
individuals . 


*Colonel,  USAF,  Tenure  Professor  and  Head  of  Aeronautics  Department 
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II.  Mass,  Weight,  and  Force 

A.  Mass 

What  is  mass?  In  all  unit  systems,  mass  has  been  taken  as  a  fundamental  or 
base  unit.  Mass  is  a  measure  of  the  inertia  of  a  body,  that  is,  a  body's  resistance 
to  acceleration.  In  a  sense,  mass  is  the  actual  body  itself,  a  characteristic  of  the 
body  that  is  independent  of  its  location.  Thus,  a  body  will  have  a  certain  constant 
mass,  regardless  of  whether  the  body  is  on  the  earth  or  on  the  moon.  Again,  by  agreement, 
a  particular  body,  a  platinum-iridium  cylinder,  was  defined  as  having  a  certain  mass, 
and  the  masses  of  all  other  bodies  are  determined  by  comparing  them  to  this  standard 
body.  The  most  common  method  for  measuring  the  mass  of  a  body  is  the  beam  balance.  By 
using  a  beam  balance,  the  mass  of  the  body  is  equal  to  the  sum  of  the  known  masses  used 
to  balance  the  beam  to  a  zero  position. 

B.  Weight 

What  is  weight?  Weight  is  a  measure  of  the  attraction  that  one  body  (for  in¬ 
stance,  the  earth)  has  for  the  mass  of  another  body.  The  defining  equations  that  allow 
one  to  determine  this  force  of  attraction  are  credited  to  Isaac  Newton.  Newton's  law 
of  gravitation  states  that  the  force  acting  on  a  body  of  mass  m  is  proportional  to  the 
mass  of  the  body  times  the  mass  of  the  second  body,  m' ,  and  inversely  proportional  to 
the  distance  between  their  gravitational  centers  squared,  r2 .  In  equation  form  this 
relationship  is  expressed  as  follows: 


F 

g 


=  G 


m  m' 


U) 


where  F^  is  the  symbol  for  the  force  acting  on  the  body,  and  G  is  a  proportionality  factor. 
The  numerical  values  of  F^,  G,  m,  and  r  depend  on  the  choice  of  a  system  of  measurement 
units.  For  example,  Fg  might  be  expressed  in  one  of  the  following  ways:  dyne,  newton, 
poundal,  or  pound-force.  Since  we  are  accustomed  to  measuring  weight  on  the  surface 
of  the  earth,  m'  and  r  remain  essentially  constant,  and  we  may  combine  G,  m'  and  r2  of 
Eqn.  (1)  into  a  simpler  form: 


Fg  =  kmg  (2) 

where  Fg  is  quite  often  given  the  symbol  W,  representing  weight,  k  is  a  constant  of  pro¬ 
portionality  dependent  on  the  system  of  units  used,  and  g  is  the  so-called  "acceleration 
of  gravity."  The  value  of  g  varies  slightly  on  the  earth,  and  different  values  of  g 
could  be  derived  from  Eqn.  (1)  for  the  moon  or  elsewhere.  Again,  the  numerical  values 
of  Fg,  k,  m,  and  g  depend  on  the  system  of  units  being  used. 
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C.  Force 

We  usually  think  of  force  in  terras  of  an  ability  to  accelerate  a  mass.  Again, 
we  are  indebted  to  Newton,  for  it  is  Newton's  second  law  that  expresses  this  relationship: 

F  =  kma  (3) 

where  F  is  the  symbol  for  the  force,  k  is  the  same  constant  of  proportionality  as  in 
Eqn.  (2)  (again  dependent  on  the  unit  system),  m  is  again  the  mass  on  which  the  force 
is  applied,  and  a  is  the  acceleration  or  change  in  velocity  of  the  body  with  time. 

Unit  systems  take  two  fundamentally  different  forms  based  on  whether  Eqn.  (2)  or 
Eqn.  (3)  is  chosen  to  define  the  units  of  force.  In  those  systems  in  which  Eqn.  (3) 
is  used  to  define  the  units  of  force,  mass  has  a  defined  unit,  as  do  length  and  time, 
but  the  unit  of  force  is  a  derived  unit,  that  is,  the  unit  of  force  is  defined  by  Eqn. 

(3).  Notice  that  if  we  choose  the  proportionality  factor,  k,  to  be  1,  then  F  =  ma,  and 
an  absolute  unit  of  force  is  defined  as  the  force  that,  when  acting  upon  one  unit  of 
mass,  will  produce  one  unit  of  acceleration.  Table  1  lists  a  few  systems  of  units  in 
which  the  unit  of  force  is  an  absolute  unit  and  is  derived  from  F  =  ma. 

There  are  systems  of  units  based  on  Eqn.  (2)  and  which  use  the  earth's  gravity  as  a 
standard.  The  force  of  attraction  of  the  earth  on  the  mass  of  a  body  on  or  near  the  earth's 
surface  is  easy  to  observe,  and  a  number  of  individuals  decided  that  it  would  be  conven¬ 
ient  to  have  a  system  of  units  that  made  use  of  this  observation.  Since  the  force  of 
attraction  for  a  given  mass  depended  on  the  location  (refer  to  m'  and  r  of  Eqn.  (1)), 
it  is  necessary  to  define  this  force  at  a  specific  location,  and  then  use  it  as  a  stand¬ 
ard.  It  was  decided  that  the  force  of  attraction  would  be  defined  at  sea  level  and  at 
45  degrees  north  latitude,  such  that  its  numerical  value  would  be  equal  to  the  numerical 
value  of  the  mass  of  the  body.  This  definition  for  force  is  referred  to  as  the  gravita¬ 
tional  unit  of  force.  Thus,  in  this  system  force,  mass,  length,  and  time  all  have  de¬ 
fined  base  units,  and  the  proportionality  factor,  k,  is  no  longer  1.  Table  2  lists  a 
few  systems  of  units  in  which  force  is  a  gravitational  unit. 

D.  Mass  and  Weight  .  .  .  There  Is  a  Difference 

Since  the  force  of  attraction  (force  of  gravity)  is  referred  to  as  weight,  W, 
then  W  =  mg  if  an  absolute  system  of  units  is  used,  and  W  =  kmg  if  a  gravitational  sys¬ 
tem  is  used.  Take  as  an  example  a  person  whose  mass  is  200  pound-mass  (lbm),  then  the 
weight  of  the  person  on  the  earth  would  be  6434  poundals  if  W  =  mg  is  used  (W  =  200  x 
32.17);  or  200  pound-force  (lbf)  if  W  =  kmg  is  used  (W  =  1/32.17  x  200  x  32.17).  Thus, 
when  the  gravitational  system  is  used,  one  finds  that  the  numerical  value  of  weight  can 
be  equal  to  the  numerical  value  of  mass  if  the  proper  units  are  used  (assuming,  of  course, 
that  one  neglects  the  small  variation  in  the  value  of  g  over  the  earth's  surface).  But, 
the  units  are  not  the  same;  that  is,  it  is  true  that  200  equals  200,  but  200  pound-force 
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Table  1 

UNIT  SYSTEMS  BASED  ON  AN  ABSOLUTE  FORCE  UNIT  (F  =  ma) 


SYSTEM 

BASE  UNITS 
(partial  list) 

DERIVED  UNITS 
(only  force  is  listed) 

NAME 

UNIT 

SYMBOL 

NAME 

UNIT 

SYMBOL 

cgs  Units 

mass 

length 

time 

gram 

centimeter 

second 

8 

cm 

s 

force 

dyne 

English 

Units 

mass 

length 

time 

pound-mass 

foot 

second 

lbm 

ft 

s 

force 

poundal 

English 

Engineering 

Units 

mass 

length 

time 

si 

ft 

s 

force 

pound-force 

lbf 

SI  Units 
( International 

System  of  Units) 

Note:  this  is  the 
new  worldwide  system 
of  units 

mass 

length 

time 

kilogram 

meter 

second 

kg 

m 

s 

force 

newton 

N 

Examples:  A  one  pound-force  (1  lbf)  is  derived  from  F  =  ma  to  be  that  force  that 
when  acting  on  a  mass  of  one  slug  (1  si)  will  produce  an  acceleration 
of  1  ft/s2. 


A  newton  is  derived  from  F  =  ma  to  be  that  force  that  when  acting 
on  a  mass  of  one  kilogram  (1  kg)  will  produce  an  acceleration 
of  1  m/s2 . 

is  not  equal  to  200  pound-mass.  Force  and  mass  are  different  things;  saying  that  200 
pound-force  equals  200  pound-mass  is  like  saying  200  apples  equals  200  oranges. 

Because  of  space  travel,  most  people  realize  that  there  is  a  difference  between 
mass  and  weight.  They  know  that  one's  weight  on  the  moon  is  approximately  one-sixth 
of  that  on  earth,  but  that  one's  mass  on  the  moon  is  the  same  as  on  the  earth.  However, 
most  people  just  state  that  their  weight  is  200  pounds  rather  than  200  pound-force.  This 
normally  doesn't  cause  any  problems  when  the  earth's  gravitational  system  is  used.  But 
with  an  absolute  system  one  must  make  a  distinction  between  mass  and  weight. 
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Table  2 

UNIT  SYSTEMS  BASED  ON  A  CRAVITAT TONAL  FORCE  UNIT  (F  =  kma) 


( 

BASE  UNITS 
partial  list) 

SYSTEM 

NAME 

UNIT 

SYMBOL 

mass 

pound-mass 

lbm 

English 

Engineering 

length 

foot 

ft 

Units 

time 

second 

s 

force 

pound-force 

lbf 

mass 

gram-mass 

cgs  Units 

length 

centimeter 

cm 

time 

second 

s 

force 

gram-force 

8f 

mks  Units 

mass 

kilogram- 

mass 

k8m 

Note:  this  is  not 
the  system  of  units 

length 

meter 

m 

that  the  world  is 

time 

second 

s 

changing  to 

force 

kilogram- 

force 

kgf 

VALUE  OF  k 


SYMBOL  NUMERICAL  VALUE  UNITS  OF  k 


lbf  •  s2 
lbm  •  ft 


980.62 


9.8062 


kg  •  m 


Example:  A  one  pound  force  (1  lbf)  has  been  defined  as  a  force  equal 
to  the  attraction  of  gravity  on  a  one  >oo*rd-mass  i  .’bm)  at 
sea  level  and  45  degrees  north  latiti;.  ■...  A  1  lbf  acting  on 
a  1  lbm  will  produce  an  acceleration  of  32.17  ft/s2. 

The  new  SI  units  is  an  absolute  system  of  units  (Ref.  2),  and  one  must  not  confuse 
mass  and  weight  (force).  In  SI  units,  the  unit  for  mass  is  the  kilogram  (kg),  and  the 
unit  for  force  is  the  newton  (N).  Thus,  for  the  person  whose  mass  is  200  pound-mass, 
the  mass  in  SI  units  is  approximately  90.7  kilograms  (see  Table  3).  The  weight  of  the 
person  would  be  889  newtons,  where  one  uses  W  =  mg  (W  =  90.7  x  9.8).  Hence,  there  is 
a  difference  between  mass  and  weight,  so  be  careful  and  correct.  You  can  talk  about 
your  mass  in  kilograms,  but  your  weight  is  in  newtons.  Most  people  will  undoubtedly 
talk  about  their  mass,  since  it  is  a  smaller  number;  however,  if  you  are  on  a  diet,  it 
will  be  harder  to  lose  mass  than  it  would  be  to  lose  weight.  For  example,  for  the  person 
whose  mass  was  90.7  kilograms  and  whose  weight  was  889  newtons,  it  will  be  easier  to 
lose  1  newton  than  to  lose  l  kilogram. 


III.  Cone lusion 

All  of  the  unit  systems  mentioned  in  this  article,  plus  several  more,  ar=  currently 
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Table  3 

HELPFUL  CONVERSION  INFORMATION 
(approximate  values) 


MASS 

FORCE 

1  kg  '  1000  g 

1  g  =  0.0022  lbm 

1  slug  *  32.17  lbm 

1  g  =  0.0353  ounce* 

1  lbf  =  32.17  poundals 

1  N  =  0.02248  lbf 

*avoirdupois 


in  use  in  the  United  States.  Needless  to  say,  these  units  have  caused  a  great  deal  of 
confusion,  and  have  caused  people  to  make  many  errors.  The  desirability  of  having  only 
one  unit  system  is  obvious.  The  new  SI  unit  system  is  superior  to  any  other  system  of 
units,  and  many  European  countries  are  converting  to  it.  Often  referred  to  as  the  metric 
system,  SI  is  a  new  system,  even  though  it  does  contain  a  number  of  aspects  of  the  old 
metric  system.  It  just  happened  that  the  U.S.  was  late  in  realizing  that  if  toe  wanted 
to  communicate  with  the  rest  of  the  world,  then  we  had  to  convert  to  SI  units.  In  fact, 
the  U.S.  cannot  afford  not  to  convert  to  SI  units  if  we  want  to  conduct  business  on  a 
worldwide  basis.  But,  let  us  do  it  right  and  try  to  avoid  embarrassments  like  the  one 
shown  in  Figure  1,  which  appear  1  on  the  back  of  a  ruler  from  the  Bureau  of  Standards 
(of  all  places). 


COM  PAX  I  Mimic  AND  CUSTOM  AIT  UN  ITT 

Unglh  Volume  Weight 


Metric  Meter 

Ci/tKJmory  Inch 

For* 

»vrd 

Rod 

Mile 


liter  _ Kilegrvm 

Teotpoon  Groin 

Tafcletpoon  Ounce 

Cup  Reend 

Pint  Ton 

Overt 
Gollon 


Metric  Unit*  or#  in  10't 

1000  millimeter*  100  centimeter*  •  I  mete* 
1000  meter*  «  I  kilometer 
1  meter  *  I  I  yorX 
I  liter  ■  I  I  quo't* 

I  kilogrom  *  7  7  pound* 

f  e  oppro*’  motel 


U  S.  DCRMTMINT  Of  COMMIftCf 
MeHenet  twreou  •* 

Wo*tH«e*on.  0  C  20734 

fcMMTO*  tuMiceN^i  )7* 

>im4  Detmbe  ttt} 


to-  teW  h  **  inpriMfiWe  <W 
Owmre,  0  4  Go— nnmm  fruiting 
Other  Woilwigm.  0  C  2<M02 
lO*4r*  b»  SO  CeW*  Me  Cl  J  103741 
hi»  IS  cent 


Figure  1.  Back  of  a  Ruler  from  the  National  Bureau  of  Standards 
Expressing  Weight  in  Kilograms 
(Notice  this  is  not  consistent  with  SI  units) 
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THE  EARLY  HISTORY  OF  THE  JET  ENGINE 
Sit  Frank  Whittle* 

Editor's  Note 

Sir  Frank  Whittle  independently  developed  the  first  successful  turbojet  engine  for 
the  Allies  during  World  War  II.  His  interest  in  jet  engines  began  while  a  cadet  at  Cran- 
well  in  1928,  and  he  patented  his  early  ideas  in  1930.  Through  his  persistence  these 
ideas,  first  rejected  as  impractical,  were  eventually  accepted  by  the  British  Air  Ministry. 
With  Whittle  as  its  head,  the  subsequent  program  lead  to  not  only  an  operating  turbojet 
engine,  but  also  to  a  successful  jet-powered  aircraft,  the  Gloster-Whittle  E-28/39.  The 
paper  that  follows  is  Sir  Frank  Whittle's  account  of  those  pioneer  days  as  he  presented 
it  to  the  faculty  and  cadets  of  the  United  States  Air  Force  Academy  on  March  13,  1975. 

Well,  gentlemen,  as  a  graduate  of  Cranwell  it's  a  great  pleasure  to  be  speaking 
at  the  Air  Force  Academy,  and  I  compliment  you  on  your  smartness.  1  don't  think  we  can 
do  better  at  Cranwell,  [but]  I  wouldn't  guarantee  it. 

Well,  all  this  jet  business  started  at  Cranwell,  though  you  young  folks,  of  course, 
probably  can't  believe  that  there  was  a  day  when  the  jet  engine  didn't  exist,  and  you 
must  wonder  how  it  was  that  it  took  so  long  to  think  of  such  a  simple  idea  and  why  there 
was  so  much  trouble  about  getting  it  going.  Well,  I'll  try  and  explain  all  that. 

I  didn't  actually  think  of  the  turbojet  while  I  was  still  at  Cranwell,  but  it  had 
its  origin  there  because  each  term  we  had  to  write  a  thesis.  In  my  first  three  terms, 

I  wrote  a  three-part  thing  on  chemistry  in  the  service  of  the  RAF.  It  dealt  with  bombs 
and  high  explosives  and  poison  gases,  and  in  my  fourth  term  I  chose  the  rather  ambitious 
subject  of  future  developments  in  aircraft  design.  .  .1  worked  out  a  formula  for  range 
and  speed,  connecting  range  and  speed,  and  engine  efficiency  and  lift-drag  ratio  of 
the  airplane,  which  I  afterwards  found  out  was  quite  a  well-known  formula,  known  as  the 
Breguet  formula.  And  I  was  very  disappointed  to  find  out  that  whatever  you  did  for  a 
given  engine  efficiency  and  a  given  lift-drag  ratio,  you  couldn't  increase  the  range. 

But  it  did  show  that  if  you  wanted  to  fly  both  fast  and  far,  you  had  to  fly  very  high, 
in  the  stratosphere.  So  I  set  out  to  try  and  find  a  power  plant  which  would  work  in 
the  stratosphere,  which  the  piston  engine  would  not  do.  And  at  a  speed  of  500  miles 
and  hour,  the  propeller  wasn't  much  good  because  of  the  high  Mach  numbers. 

Well,  I  thought  about  rockets,  and  the  fuel  consumption,  or  course,  is  much  too 
high  for  an  airplane.  And  I  took  a  look  at  gas  turbine-driven  propellers,  but,  of  course, 
there  the  trouble  was,  as  I  just  mentioned,  the  propeller.  So  I  kept  on  working  away 
after  I  graduated  from  Cranwell.  While  I  was  a  pilot  In  111  Squadron,  1  thought  of  a 
piston  engine  inside  a  hollow  fuselage,  driving  a  low  pressure  compressor  producing  a 
jet  -  the  same  idea  that  the  Italians  produced  in  1940;  I  think  it  was  the  Caprono  - 
Campini's  scheme.  But,  when  I  worked  out  what  the  fuel  consumption  was  and  so  on,  it 
didn't  seem  to  me  to  be  of  much  good. 


*Air  Commodore,  RAF  (ret) 
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Then,  when  I  was  undergoing  training  as  a  flying  instructor  at  the  Central  Flying 
School,  I  suddenly  felt,  "Well,  why  not  throw  the  piston  engine  away  and  substitute  a 
turbine  and  drive  a  high  pressure  compressor  instead  of  the  low  pressure  compressor," 
and  that  was  the  turbojet. 

Now  I  should,  in  all  fairness,  say  that  I  was  not  the  first  to  think  of  the  turbo¬ 
jet  as  an  engine.  There  was  a  Frenchman  named  Maxime  Guillaume  who  took  out  a  patent 
in  1921.  Of  course,  it  was  long  before  its  time. 

Whatever  happened  to  Maxime  Guillaume  I  never  knew,  ana  whatever  he  did  about 
it,  1  don't  know.  But  anyway,  1  got  a  patent  in  Britain  in  1930,  after  the  Air  Ministry 
had  turned  the  thing  down.  They  had  already  made  up  their  minds  that  the  gas  turbine 
was  a  dead  duck,  and  all  the  firms  I  tried  to  sell  the  idea  to  decided  that  the  gas 
turbine  wouldn't  work. 

I  took  out  the  patent  and  got  it,  and  the  first  [figure]  shows  that  (Figure  1). 

I  am  showing  this  because  ...  it  isn't  very  representative  of  what  a  jet  engine 
really  looked  like,  but  unfortunately  I  haven't  got  with  me  a  picture  which  shows 
a  section  of  the  early  engines.  But  this  patent  drawing  does  give  some  clue  as  to  what 
form  the  engine  was  going  to  take. 

At  the  top  is  the  pressure  volume  diagram  for  the  stationary  condition,  and 
the  curve  -  it's  off  the  picture,  I'm  afraid  -  the  Curve  D  to  C,  I  think  it  is, 
represents  the  way  the  pressure  is  increased  and  volume  reduced  in  the  compression 
process,  and  from  C  to  E  you  see  the  increase  in  volume  due  to  heating  in  the  combus¬ 
tion  chamber  at  constant  pressure.  And  then  from  E  to  F  is  the  expansion,  the  drop 
in  pressure  through  the  turbine,  enough  to  drive  the  compressor,  and  then  that  leaves 
F  to  G  as  the  balance  of  energy  for  producing  the  energy  of  the  propelling  jet. 

The  machinery  for  doing  the  job  is  shown  here  very  crudely.  I  was  sitting  on 
fence  where  the  compressor  was  concerned.  Two  stages  of  axial-flow  compressors 
were  followed  by  a  single-stage  centrifugal,  a  multiple  combustion  chamber  system,  a 
two-row  Curtis  turbine,  and  then  the  multiple  exhaust  nozzles.  Of  course,  in  practice 
we  only  use  one. 

This  drawing  was  deliberately  made  to  look  a  little  bit,  1  would  say,  "childish," 
because  I  thought  the  thing  ought  to  be  on  the  secret  list,  so  I  thought  the  best 
thing  was  a  drawing  that  was  good  enough  for  patent  purposes,  but  which  wouldn't 
give  away  to  much  to  anyone  who  got  hold  of  the  patent. 

Well,  nothing  happened  for  five  years,  and  no  one  would  believe  it.  You  see,  the 
great  problem  was  this:  in  those  days  the  efficiency  of  compressors  was  so  poor  that, 
say  if  you  needed  1000  horsepower,  ideally,  in  fact,  you  had  to  use  1500  horsepower  to 
drive  it.  And  turbine  efficiencies  were  low  too,  so  where  ideally  a  turbine  might  pro- 
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Figure  1.  The  Turbojet  Engine  in  Whittle's  1930  Patent 


duce,  say,  2000  horsepower,  it  only  produced  1500.  So  when  you  deducted  the  work  required 
to  drive  the  compressor  from  the  work  produced  by  the  turbine,  you  didn't  have  much  left. 
In  fact,  you  were  very  lucky  if  the  thing  would  even  turn  itself  over. 

And  incidentally,  what  we  call  the  positive/negative  ratio  and  the  work  depended 
very  much  on  the  maximum  temperature  you  could  use  in  the  turbine,  and  that  in  turn  was 
governed  by  the  materials  available,  and  that  was  a  very  severely  limiting  factor.  But 
I  was  pretty  confident  that  we  would  get  materials  which  would  stand  the  temperature, 
and  I  was  even  more  confident  that  I  could  improve  compressors  and  turbines  so  that  this 
big  objection  to  the  gas  turbine  would  disappear.  And  I  further  realized  that  as  we 
got  bigger,  the  cold  air  of  high  altitude  would  make  a  very  big  difference  to  the  pros¬ 
pects  where  efficiency  was  concerned. 

Well,  as  I  said  just  now,  I  was  very  confident  about  being  able  to  get  an  efficient 
compressor.  .  .too  confident,  I  will  admit.  I  was  equally  confident  to  do  much  better 
on  the  efficiency  of  the  turbine,  but  I  was  far  from  sure  that  we  would  get  satisfactory 
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combustion,  and  that  did  in  fact  prove  to  be  our  biggest  difficulty.  And  this  (figure] 
(Figure  2)  shows  you  the  kind  of  conditions  we  had  to  work  under  in  trying  to  develop 
the  combustion  system. 


Figure  2.  Early  Laboratory  Used  to  Develop  the  Combustion  System 

This  was  outside  the  BTH  turbine  factory,  in  the  open  air.  Above  it  there  was 
an  office  where  the  planning  engineers  and  their  staff  worked,  and  we  were  most  unpopular 
with  them  because  that  thing  here  (see  Figure  2)  is  the  combustion  chamber,  that's  the 
fuel  container,  and  here  is  a  motor  driven  fuel  pump,  and  the  burner  was  in  here. 

Well,  that  thing  used  to  behave  more  like  a  machine  gun  than  a  combustion  chamber. 

In  fact,  on  one  occasion,  the  people  in  the  office  above  had  gone  away  to  lunch,  and 
when  they  came  back  and  found  that  everything  on  the  desk  -  the  inkpots,  their  pens, 
papers,  and  everything  else  -  was  on  the  floor,  they  thought  some  practical  joker  had 
been  around  and  done  all  this,  until  we  started  it  up,  and  of  course  the  concussions 
coming  up  from  below  were  causing  everything  to  "walk"  in  the  office.  On  other  occasions 
we  smoked  them  out. 

Well,  the  whole  job  was  hung  up  for  five  years  because  nobody  would  take  any  interest 
in  it,  but  it  did  happen  that  someone  who  had  been  a  cadet  with  me  at  Cranwell  -  this 
is  very  much  an  Air  Frrce  job,  you  know  -  Williams  (he  is  now  Sir  Dudley  Williams),  he 
was  in  my  term  at  Crowell,  and  we  were  in  the  same  hut  together,  and  he  knew  about  the 
thing.  In  the  Royal  Air  Force,  later  on,  it  was  rather  jokingly  referred  to  as  "Whittle's 
flaming  touchhole."  No  one  took  it  very  seriously,  but  after  five  years  Williams  and 
another  officer  who  had  been  invalided  out  after  a  bad  crash,  came  along  and  said  they 
thought  they  could  raise  some  money.  They  managed  to  interest  a  firm  of  city  bankers. 
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and  a  little  company  called  Power  Jets  Limited  was  formed.  At  the  time,  I  was  at  Cambridge 
University  being  trained  as  an  engineer  officer  following  my  flying  career.  We  managed 
to  get  enough  money  to  start  the  company,  and  we  got  enough  of  it  to  put  an  order  with 
the  British  Thompson  Houston  Company  to  build  an  engine  to  my  design  while  I  was  still 
at  Cambridge. 

The  Air  Ministry  still  didn't  believe  in  it,  but  fortunately,  the  Air  Member  for 
Personnel  agreed  to  let  me  stay  at  Cambridge  for  a  post-graduate  year,  knowing  that  in 
reality  I  was  working  on  the  job,  though  on  paper  1  was  doing  aerodynamics  research  for 
Professor  Jones  -  Sir  Melville  Jones,  as  he  was  later.  1  was,  in  fact,  doing  research 
work  with  Sir  Melville  Jones,  but  I  was  spending  most  of  my  time  on  the  jet  en¬ 
gine. 

[Figure  3  shows]  the  first  engine.  For  some  reason  or  other,  when  the  photograph 
was  taken,  they  blotted  out  the  stand  on  which  it  was  mounted,  but  you'll  see  that  in 
a  later  [figure].  This  one  had  a  single  combustion  chamber,  because  I  was  very  nervous 
about  using  a  multiple  combustion  chamber  system  for  reasons  I  will  explain.  [This 
figure  also  shows]  auxiliary  things  -  a  water  pump,  because  we  had  a  water  cooling  jacket 
each  side  of  the  turbine,  and  a  fuel  pump,  and  so  forth  -  and  this  was  a  piece  of  plywood 
with  the  instruments  mounted  on  it,  and  here  was  the  control  valve.  All  this  was  in 
the  open  factory  in  the  BTH  initially,  and  it  was  later  thrown  out  by  the  chief  engineer 
because  he  thought  it  was  too  dangerous. 

When  we  first  tried  to  start  up  the  engine,  1  gave  a  signal  to  light  up,  and  some¬ 
one  looking  through  a  quartz  window  in  the  combustion  chamber  told  me  that  it  had  lit 
up.  Then  we  speeded  the  thing  up  with  an  electric  motor  to  2000  RPM,  then  I  opened 
the  throttle,  and  the  thing  accelerated,  and  it  accelerated  out  of  control.  And  so  did 
everyone  standing  around  it.  I  didn't,  because  I  was  paralyzed  with  fright.  It's  a 


Figure  3.  First  Engine 
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funny  thing  -  if  you  are  actually  in  control,  that’s  the  kind  of  reaction. 

Well,  this  happened  several  times,  and  I  began  to  think  that  we  had  something  which 
was  fundamentally  uncontrollable.  The  next  night,  which  was  a  Tuesday  evening,  not  only 
did  the  thing  accelerate  out  of  control  on  this  second  occasion,  because  the  heat  had 
weakened  the  spring  in  the  fuel  injector,  (but]  there  were  leaking  joints  at  the  top 
of  the  thing  and  flames  in  midair.  It  was  pretty  terrifying.  The  thing  went  off  making 
a  noise  like  an  air  raid  siren.  I  badly  needed  half  a  bottle  of  wine  that  night. 

The  next  day  I  was  walking  down  toward  the  turbine  factory  with  the  deputy  chief 
engineer  of  the  turbine  factory.  He  was  a  north  countryman  with  a  broad  Yorkshire  accent 
which  is  probably  not  very  familiar  to  you  -  and  he  hadn’t  seen  the  engine  run.  1  told 
him  how  frightened  1  had  been  the  night  before.  He  said,  "Ee,  lad,  you  don't  know  what 
it  is  to  be  frightened.  You  should  have  been  standing  next  to  the  Blackdown  job  when 
it  jumped  out  of  its  bearings."  That  was  a  vertical  turbine. 

Anyway,  he  was  present  when  we  made  our  first  experiment  with  a  vaporizer  in  the 
exhaust  pipe,  using  vaporized  kerosene  instead  of  diesel  fuel. 

The  thing  went  fine,  at  first.  It  went  up  to  4000  RPM,  beautifully  under  control. 
Then  the  vaporizer  primed,  so  that  instead  of  getting  kerosene  vapor  in  the  combustion 
chamber,  we  got  liquid  kerosene.  And  off  she  went  again.  And  off  went  Collingham  at 
high  speed,  and  he  got  jammed  in  a  door  leading  out  of  the  factory,  competing  with  some¬ 
one  else  trying  to  get  out  at  the  same  time.  And  the  chief  engineer  of  the  turbine  fac¬ 
tory  was  fighting  for  the  possession  of  a  turbine  exhaust  casing  with  one  of  the  test 
hands.  It  was  pretty  exciting! 

On  another  occasion  only  a  few  days  later,  the  compressor  fouled  the  casing,  and 
the  whole  thing  came  to  a  screeching  halt  from  12,000  RPM  in  one  and  a  half  seconds. 

One  of  my  colleagues  turned  to  run  when  this  terrible  screech  happened.  Before  he  had 
even  taken  a  pace  and  a  half,  the  thing  was  stopped.  The  shop  superintendant  put  a  hand 
on  his  shoulder  and  said,  "It's  no  use  running,  sir.  It'll  soon  catch  you  if  it  wants." 
What  made  me  laugh  was  this  man,  Walter  Smith,  had  himself  down  the  factory  at  high  speed 
four  times,  at  least. 

Anyway,  after  that  the  chief  engineer  of  the  whole  BTH  turned  us  out  -  he  thought 
it  was  too  dangerous  having  this  thing  running  in  the  open  factory.  But  we  did  manage 
to  rent  an  old  iron  foundry  3even  miles  away,  and  there  we  carried  on  with  our  test  work. 

Well,  that  first  engine  wasn't  very  satisfactory.  The  first  run  was  April  1937. 

It  was,  in  fact.  .  .the  first  jet  engine  ever  to  run,  though  the  Germans  were  only  a 

few  weeks  later.  Their  first  engine  ran  in  September  of  1937. 

We  reconstructed  it  twice  in  1937  and  1938.  The  first  reconstruction  didn't  work 
at  all  well.  We  had  a  turbine  failure.  So  we  had  to  reconstruct  it  a  second  time,  and 

this  time  I  took  the  bold  step  of  having  a  multiple  combustion  chamber  system  (see  Figure 

4),  You  wouldn't  think  it  was  the  same  engine,  but  it  is,  except  for  the  sheet  metal 
work.  All  the  sheet  metal  work  was  new.  Of  course,  we  had  to  have  a  new  turbine  because 
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Figure  4.  Reconstructed  Engine  with  Multiple 
Combustion  Chambers 

we'd  just  wrecked  one,  but  otherwise  all  this  front  end  is  the  same,  the  compressor  cas¬ 
ing  is  the  same,  the  rotor  shaft  and  all  the  rest  of  it  were  the  same  as  in  the  first 
engine,  because  we  couldn't  afford  to  do  anything  else.  We  lost  a  lot  of  time  because 
we  couldn't  afford  to  have  new  parts  as  we  should  have  done. 

We  were  now  starting  up  with  a  10-horsepower  BSA  car  engine  which  was  connected 
by  a  pulley.  The  petrol  tank  (from  a  motorcycle)  [was]  for  the  car  engine,  [and]  the 
radiator  [was]  for  the  starter  motor. 

Well,  we  found  about  a  year  later  that  we  could  start  this  thing  up  with  the  starter 
motor  of  this  starting  engine.  The  horsepower  [of  the]  electric  motor  was  enough  to 
start  up  the  experimental  engine. 

We  couldn't  get  much  support  from  the  Air  Ministry  at  first  -  they  looked  upon  this 
as  long-term  research.  They  did  put  me  on  the  special  duty  list  after  the  end  of  my 
post-graduate  year  at  Cambridge,  but  at  first  they  wouldn't  do  anything  else  to  help 
us.  Then  they  agreed  to  pay  about  5000-fc  for  so  many  hours  of  running  and  a  report 
and  the  reconstruction,  and  then  a  short  time  later  they  agreed  to  pay  for  those  test 
houses.  The  test  houses  were  built  for  lO.OOOt  -  two  test  houses  (see  Figure  5). 

We  had  a  large  Venturi  tube  inside  so  we  could  measure  the  airflow  into  the  engine. 
But  the  real  point  of  showing  you  (Figure  5]  is  to  show  you  the  conditions  we  had  to 
work  under.  .  .all  this  rubble  all  over  the  place.  It  was  like  working  in  a  junk  yard, 
and  the  engine  wasn't  much  more  than  junk  either. 

[Figure  6  does  not  show]  the  actual  rotor  of  the  experimental  engine,  but  it's  suf¬ 
ficiently  similar  to  explain  it.  As  you  can  see,  we  had  a  single-stage,  double-sided 
centrifugal  compressor,  driven  by  a  single-stage  turbine.  That  was  really  the  main 
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rotating  part.  Since  no  lady  is  present,  I  think  I  can  make  a  comment  I  have  made  at 
various  times:  that  reciprocating  motion  is  fine  in  biology,  but  I  have  no  use  for  it 
in  engineering. 

As  you  can  see,  it's  pretty  simple.  It  was  simple  then.  They  are  not  so  simple 
these  days.  In  fact,  when  Hives  of  Rolls-Royce  came  to  see  the  engine,  when  first  of 
all  we  told  him  the  thrust  was  about  1000  pounds,  he  said,  "Poo!  That's  not  enough  to 
blow  the  skin  off  a  rice  pudding."  But  one  of  his  engineers  -  now  Sir  Stanley  Hooker  - 
pointed  out  to  him  that  that  made  it  more  powerful  than  the  Merlin,  which  was  the  main 
engine  of  the  Royal  Air  Force  in  those  days.  And  that  made  Hives  change  his  mind  about 
jet  engines.  And  then,  when  I  was  explaining  to  him  the  simplicity  of  the  thing,  he 
said,  "We'll  soon  design  the  bloody  simplicity  out  of  it."  And  they  succeeded,  because 
the  more  complicated  it  is,  the  more  you  can  charge  for  it,  you  see. 

[Figure  7  shows]  what  happened  to  the  experimental  engine  after  a  very  useful  four 


Figure  7.  Failed  Turbine  with  de  Lavale 
Blade  Fittings 

years  of  development  running.  The  first  turbine  had  the  de  Lavale-type  [fittings]  for 
the  turbine  blades,  but  they  just  weren't  good  enough.  Before  this  failed,  we  had  de¬ 
vised  what  we  called  the  "fir  tree"  [fitting],  which  is  now  universally  used.  This  was 
the  end  of  about  a  100-hour  test  -  the  blades  came  out.  It  was  getting  pretty  bad  before 
it  actually  failed;  so  bad,  in  fact,  that  every  time  the  engine  came  to  rest  you  could 
hear  a  tinkling  noise  as  the  blades  flopped  over  in  their  slots. 

At  the  end  of  June  1939,  the  Director  of  Scientific  Research,  who  up  to  then  had 
called  this  a  long-term  research  thing,  came  up  to  see  it  run,  and  he  was  greatly  impressed. 
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I  had  the  curious  experience,  driving  him  to  catch  his  train,  when  all  the  way  to  the 
station  he  was  telling  me  all  the  advantages  of  it:  how  smooth  it  was,  that  it  could 
run  on  almost  any  fuel,  et  cetera,  et  cetera.  1  was  only  a  squadron  leader  then,  so 
I  wasn't  going  to  argue.  He  decided  there  and  then  that  we  must  have  a  flight  engine 
and  an  airplane  to  put  it  in,  and  they  later  ruled  that  the  flight  engine  must  be  as 
similar  to  the  experimental  engine  as  possible,  and  this  is  it.  [Figure  8  shows]  the 
W— 1 ,  the  first  flight  engine.  It's  very  similar  to  the  experimental  engine,  except  that 
is  was  a  lot  lighter.  The  experimental  engine  weighed  about  1000  pounds,  whereas  this 
weighed  about  620  pounds.  It  was  designed  for  1240  pounds  of  thrust,  but  it  was  derated 
to  860  pounds  for  the  first  flight.  It  flew  once  at  1000  pounds  during  flight  trials. 

It  did  eventual !v  reach  1240  pounds  before  we  eventually  retired  it. 


Figure  8.  The  W-l  Engine 

[Figure  9  shows]  the  airplane  that  was  built  at  the  same  time,  the  Gloster-Whitt le 
E28.  It  was  a  highly  successful  airplane.  It  flew  with  several  engines  one  after  the 
other,  at  first  with  the  Wl-X,  then  with  the  W-l,  then  with  the  Wl-A,  then  with  the  W2-B. 
then  the  W2-500,  and  eventually  the  W2-700,  which  was  too  powerful  for  it  in  fact,  and 
couldn't  be  opened  up  fully.  With  the  W-l  it  reached  a  speed  of  about  370,  which  was 
faster  than  the  Spitfire.  With  the  W2-700  it  did  well  over  500,  and  then  started  to 
run  into  compressibility  troubles,  and  is  intact  in  the  Science  Museum  in  London  today. 

The  other  one,  the  thin-winged  one  -  the  second  one  that  was  built  -  was  lost  because 
the  ailerons  seized  up  through  contraction  at  the  height  of  32,000  feet,  and  Squadron 
Leader  Davey,  the  pilot,  had  to  bail  out. 
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Figure  9.  The  Gloster-Whittle  E28 
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Figure  10.  Fuselage  Arrangement  of  the 
Experimental  Airplane  E28/39 
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[Figure  10]  shows  you  the  general  arrangement  of  that  airplane.  As  you  can  see, 
the  air  was  taken  in  right  at  the  front  and  passed  around  the  pilot  through  two  ducts 
and  into  the  engine  chamber,  which  was  behind  the  fuel  tank.  I  don't  know  how  the  pilot 
felt  about  having  a  fuel  tank  right  behind  him,  but  there  it  was.  It  never  did  him 
any  harm. 

Now  everyone  sees  engines  in  the  tail  these  days,  but  1  always  wanted  them  to  be 
at  the  back  of  the  airplane  (see  Figure  11).  That  seemed  to  me  to  be  the  most  appropriate 
place.  But  few  people  realize  that  we  did,  in  fact,  fly  the  first  airplane  with  an  engine 
in  its  tail.  That  was  the  Wellington,  used  as  a  flying  test  bed.  We  were  lent  a  Welling¬ 
ton. 

By  this  time,  we  had  our  own  flight  testing  facilities,  including  the  Wellington, 
and  we  replaced  the  tail  gun  turret  with  the  test  cell.  The  Wellington,  as  you  probably 


Figure  11.  Jet  Engine  Installed  in  the  Tail  of  a  Wellington  Bomber 

know,  was  a  twin-engined  bomber,  powered  by  two  Bristol  radial  engines  (with  propellers, 
of  course).  Well,  we  took  out  the  tail  gunner's  turret  and  put  in  the  W2-B  engine,  which 
was  the  production  version  of  the  engine,  and  used  this  as  a  flying  test  bed. 

Well,  the  story  goes  -  1  don't  know  how  true  it  is  -  a  B-17,  which,  as  you  know, 

was  the  standard  day  bomber  of  the  U.S.  Air  Force  (the  U.S.  Army  Air  Force,  as  it  was 

in  those  days),  a  four-engine  bomber,  was  flying  along  when  along  came  a  Wellington  with 

its  two-piston  engines  going  and  the  W2-B  in  the  back.  The  story  goes  that  the  pilot 

for  a  joke  cut  his  two-piston  engines  and  feathered  the  propellers,  opened  up  the  W2-B 
and  overtook  the  B-17.  So  the  crew  of  the  B-17  saw  a  Wellington  with  both  engines  failed, 
as  they  thought,  go  shooting  past  them.  They  probably  decided  to  take  more  water  with  it. 
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At  the  same  time  as  we  were  doing  all  this,  the  Germans  were  at  it  too.  I've  often 
said  that  inventions  are  something  which  happen  in  several  places  at  much  the  same  time. 
Ideas  ripen  rather  like  fruit;  you  get  several  people  getting  the  same  idea  at  the  same 
time,  and  no  doubt  many  other  people  thought  of  jet  engines  at  the  same  time  as  I  did. 
The  Germans  were  certainly  doing  it  too.  A  lot  of  people  have  written  to  me  saying  that 
they  had  thought  of  the  jet  engine,  and  I  don't  doubt  they  did,  but  thinking  about  it 
is  one  thing,  and  making  it  work  is  quite  another. 

Well,  von  Ohain  was  at  work  in  Germany,  and  [Figure  12  shows]  his  engine,  which 
was  quite  similar  to  ours  in  many  ways,  except  that  he  used  a  radial  in-flow  turbine  in¬ 
stead  of  an  axial~flow  turbine.  And  though  his  engine  ran  a  little  later  than  ours, 
they  did,  in  fact,  get  an  airplane  into  the  air  before  we  did.  They  did  about  three  or 


Figure  12.  Schematic  of  von  Ohain's  Jet  Engine 

four  flights.  I  think  I  am  safe  to  say  it  wasn't  anything  like  as  successful  as  the 
E28,  because  the  type  of  engine  used  in  the  E28  went  into  production,  but  von  Ohain's 
engine  never  did.  The  Germans  then  changed  over  to  the  axial-flow  type  of  engine,  the 
Junkers  004,  which  powered  the  Messerschmid t  Me  262. 

[Figure  13  shows]  the  German  airplane  which  flew  in  August  1939,  the  Heinkel  He  178. 
An  this  is  what  von  Ohain  wrote  on  this  thing  when  he  gave  it  to  me.  Very  nice! 

[Shown  in  Figure  14  is].  .  .a  later  version  of  von  Ohain's  engine.  I  can't  really 
tell  you  much  about  it.  Whether  it  actually  flew  or  not,  I  don't  know.  1  don't  think 
it  did,  but  I  may  be  wrong. 

[Figure  15  shows]  the  Me  262,  which  was  really  the  first  jet  airplane  to  go  into 
production.  I  wouldn't  call  it  a  very  safe  airplane.  It  did  a  lot  of  damage  to  American 
day  bombers.  It  didn't  do  any  damage  to  the  RAF  night  bombers  because  they  never  sent 
them  up  at  night.  But  it  was  a  real  menace  to  the  American  day  bombers,  so  much  so  that 
our  Meteors  were  withdrawn  from  combat  duty  for  tactical  exercises  to  try  and  find  out 
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Figure  13.  German  He  178  Aircraft  with  Inscription  by  Hans  von  Ohain 
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Figure  14.  Later  Version  of  von  Ohain’ s  Engine 

a  way  of  dealing  with  jet  fighters,  and  they  didn't  succeed  in  finding  out  any  way. 

But  nevertheless,  I  think  it's  true  that  these  things  killed  more  of  their  own  pilots 
than  they  did  anyone  else's  because  they  had  a  very  high  wing  loading  -  it  took  a  hell 
of  a  long  runway  to  get  off  the  ground  -  and  the  engine  had  a  most  unfortunate  character¬ 
istic:  if  you  tried  to  open  up  quickly,  the  compressor  surged  and  you  lost  power  instead 
of  gaining  it.  And  believe  me,  when  you  are  coming  in  and  are  under  shooting,  and  you 
want  to  open  up  your  engine,  and  instead  of  gaining  power  you  lost  it,  there  just  is 
no  future  in  it.  An  that's  what  happened  several  times,  I  understand,  to  these  airplanes. 


USAFA-TR-81-4 


MESSERSCHMITT  Me  262A 


Figure  15.  Messerschmitt  Me  262A 

After  the  war,  we  had  a  captured  one  at  Farnborough.  It  took  off,  taking  the  whole  of 
the  Farnborough  runway  and  just  cleared  the  trees  at  the  far  end.  A  wing  commander, 
a  friend  of  mine,  turned  to  me  with  a  shudder  and  said,  "Now  we  know  why  Hitler  wanted 
to  extend  the  Third  Reich." 

In  1941,  before  America  entered  the  war,  the  British  and  American  governments  got 
together  to  arrange  for  the  production  of  jet  aircraft  in  America.  The  main  reason, 
from  our  point  of  view,  was  that  we  were  afraid  of  being  bombed  out,  so  we  wanted  to 
back  up  what  we  were  doing  with  production  in  the  States. 

On  October  1st,  1941,  we  flew  over  the  Wl-X  engine  with  an  engineer  and  two  techni¬ 
cians  to  the  General  Electric  Company.  We  selected  them  because  they  had  already  gotten 
a  lot  of  experience  with  turbo-superchargers,  which  were  rather  similar,  and  they  got 
on  with  the  job  very  quickly.  By  April  of  1942,  before  I  went  over  in  June,  they  had 
gotten  the  engine  on  the  bench  -  the  1-14,  which  was  their  version  of  our  W2-B.  Bell 
Aircraft  were  getting  on  so  fast  with  the  airplane  that  they  had  it  flying  on  October 
2nd  of  1942.  In  fact,  the  Bell  airplane  flew  before  our  twin-engined  Meteor,  but  it 
never  went  into  combat. 

Well,  [Figure  16]  is  it.  As  you  can  see,  it  had  its  two  engines  tucked  in,  and 
for  the  first  time  we  used  these  plates  to  prevent  the  boundary  layer  from  fouling  up 
the  flow  into  the  intake.  I  was  responsible  for  that  particular  feature  on  this  thing. 

[Figure  17]  is  one  of  our  later  engines,  the  W2-700,  which,  as  I  have  already  told 
you,  was  rather  too  powerful  for  the  E28.  It  developed  i  thrust  of  2000  pounds  initially 
and  reached  2500  by  the  time  we'd  developed  it. 

[Figure  18]  shows  the  production  airplane  which  was  powered  by  Rolls-Royce  Wellands, 
the  Rolls-Royce  version  of  the  W2-B. 
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BELL  XP-59A 


Figure  16.  Bell  XP-59A 


Figure  17.  The  W2-700  Engine 

Throughout  1942  we  had  developed  terrible  problems  all  the  time,  until  Rolls-Royce 
took  over  from  the  Rover  Company,  who  had  been  selected  for  production.  When  Rolls- 
Royce  took  over.  Sir  Stanley  Hooker  became  the  boss,  as  it  were,  of  the  technical  side, 
and  things  really  began  to  move,  and  we  soon  had  the  Meteor  in  production. 

They  went  into  operation  in  July  of  1944.  Their  first  use  was  a  quite  unexpected 
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Figure  18.  First  Production  Airplane,  The  Meteor, 

Powered  by  Two  W2-B  Engines 

one  -  they  were  used  against  the  flying  bombs,  because  they  were  the  only  airplanes  which 
could  catch  them.  Then,  after  this,  as  I've  already  told  you,  they  were  withdrawn  for 
operational  practice  and  exercises  to  see  whether  jet  airplanes  could  be  dealt  with. 

I  flew  the  Meteor  myself,  but  not  until  1945,  and  then  I  had  to  do  it  by  stealth 
because  the  Ministry  wouldn't  give  me  permission  to  fly  it.  But  we  had  our  own  airfield 
then,  and  all  the  people  there  were  subordinate  to  me,  so  I  told  them  I  was  going  to 
get  into  the  airplane  and  taxi  it.  There  was  no  one  there  who  was  able  to  say  "no," 
and  you  know,  it  just  took  off. 

Well,  Rolls-Royce  really  got  on  with  the  job.  They  built  the  Nene  engine,  which 
had  a  thrust  of  5000  pounds.  One  evening  Hooker  and  another  Rolls-Royce  engineer  and 
myself  were  having  drinks  in  a  pub  in  Lancashire,  and  one  of  us  -  I  don't  know  who  - 

said,  "Why  don't  we  scale  down  the  Nene  until  it  will  fit  the  Meteor?" 

So  they  did  that,  and  six  months  after  that  conversation  they  built  the  Derwent  V 
(see  Figure  19),  with  which  we  got  the  airspeed  record  in  November  1945. 

I  flew  the  Meteor  again  -  the  Meteor  3  on  that  occasion;  not  the  one  with  the  Derwent 

V,  but  one  with  the  Derwent  1.  In  order  to  keep  the  press  quiet  while  a  delay  was  going 

on  in  the  attempt  at  the  speed  record,  they  asked  me  to  fly  the  Meteor  up  and  down  the 

high  speed  course.  I  had  the  great  pleasure  of  flying  at  50  feet,  flat  out,  up  and  down 
the  high  speed  course.  I  had  the  wind  taken  out  of  my  sails,  however,  when  my  younger 
son  remarked,  "What?  Only  450?  It  will  do  600,  you  know,"  to  me.  The  Derwent  V,  inci¬ 
dentally,  was  made  in  thousands. 

Now  [Figure  20]  is  just  to  show  you  that  we  were  also  working  on  more  advanced  con¬ 
cepts.  This  was  the  basic  patent  for  the  turbofan  which  1  took  out  in  1936.  It  doesn't 
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look  much  like  the  turbofans  you  see  these  days  because  in  order  to  avoid  the  complica¬ 
tions  of  shafts  going  through  shafts,  1  had  [the]  gas-generating  unit  cross-wise  -  cross¬ 
wise  in  the  fuselage  -  and  two  nacelles  with  the  fans  driven  by  their  own  turbines. 

Not  many  peope  realize  that  Power  Jets  were  responsible  for  turbofans  and  aft-fans 
as  well  as  the  simple  turbojet.  The  big  disadvantage  of  the  simple  turbojet  is  that 
its  propulsive  efficiency  is  very  low.  So  I  always  wanted  to  be  able  to  gear  down  the 
jet;  in  other  words,  to  turn  a  low-mass,  high-velocity  jet  into  a  high-mass,  low-velocity 
jet,  and  that’s  the  whole  object  of  the  turbofan. 

To  make  the  turbofan  engine  was  a  major  proposition.  We  couldn't  spare  the  manpower, 
and  we  didn't  have  the  money,  so  we  rigged  up  a  scheme  that  we  could  tack  onto  an  existing 
engine.  We  called  that  the  af t-fan,  [ shown  in  Figure  21].  It's  simply  something  you 
add  to  the  back  end  of  the  engine  and  use  the  exhaust  of  the  main  engine  to  drive  turbines 
independent  of  the  main  engine,  which  drive  a  low-pressure  compressor,  which  gears  down 
your  jet  more  or  less,  so  you  get  a  cold  jetstream  there  and  a  hot  jetstream  here  and 
a  much  better  propulsive  efficiency.  This  we  call  the  No.  2  Thrust  Augmenter.  We  never 
actually  built  it,  but  we  sent  the  drawings  to  General  Electric,  and  they,  some  time 
later,  did  in  fact  build  it.  The  engines  in  the  Convair  990  have  them,  and  the  Phantom 


Figure  21.  Schematic  of  Aft  Fan,  the  No.  2  Thrust  Augmenter 

Fanjet  has  it,  and  other  General  Electric  engines  have  it. 

[Figure  22  shows  a  different  turbofan  design.]  It's  not  a  good  thing  to  have  the 
turbine  running  slower  than  the  fan,  which  is  inevitable  in  the  arrangement  you've  just 
seen.  It's  much  better  to  have  it  the  other  way  around.  So  we  wanted  to  do  a  tip  tur¬ 
bine  job,  and  have  a  tip  turbine  driving  the  fan  instead  of  the  other  way  around.  The 
big  snag  there  was  ducting  the  hot  exhaust  from  the  main  engine  to  the  tip  turbine  with¬ 
out  fouling  up  the  flow  into  the  fan. 
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Figure  22.  Alternate  Turbofan  Design 

[Figure  23]  shows  you  what  we  were  trying  to  do.  We  built  this  one.  It  did,  in 
fact,  give  a  40-percent  improvement,  but  since  we  were  expecting  far  more  than  that, 
we  didn't  go  on  with  it  because  [of]  the  ducts  leading  to  the  blades  of  the  tip  turbine. 
It  was  a  single-stage  job.  This  is  the  tip  turbine,  with  its  blading  on  the  fan  blades, 
but  the  ducting  did,  in  fact,  foul  up  the  flow  into  the  fan.  We  could  have  done  some¬ 
thing  about  it  if  we  had  the  time  and  the  manpower,  but  we  had  the  much  more  urgent  job 
of  getting  on  with  the  simple  jet. 

The  turbofan  was  obviously  desirable  for  a  long  range  airplane.  [Figure  24]  shows 
the  civil  version  of  an  airplane  which  was  based  on  what  we  then  regarded  as  what  was 
needed  for  bomber  operations  in  the  Pacific.  I  went  to  the  Staff  College  in  1943,  and 
I  took  the  opportunity  of  having  a  meeting  of  the  students  and  the  directing  staff.  We 
had  a  conference  to  try  and  de  ;ine  what  was  needed  for  a  Pacific  bomber,  and  we  came 
up  with  4000  miles  still  air  range,  a  speed  of  450  miles  an  hour,  and  so  on,  a  bomb  load 
of  10,000  pounds,  which,  of  course,  also  made  very  good  civil  airliner  specification. 

And  that  was  the  airplane  on  which  we  based  the  design  of  the  LR1  engine  (Figure  25). 

[The  LR1  ducted-fan  engine]  was  the  last  engine  I  designed.  This  is  just  a  diagram 
of  it,  and  you  might  say  it  was  the  first  turbofan  to  be  designed.  Unfortunately,  Power 
Jets  by  then  had  been  nationalized,  and  one  consequence  of  nationalization  was  that  the 
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Figure  24.  Transatlantic  Civil  Airliner 
with  Four  LR1  Ducted-Fan  Engines 
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Figure  25.  LR1  Ducted-Fan  Engine 


contract  for  this  thing  was  cancelled.  In  my  opinion,  this  cost  20  years  in  the  develop¬ 
ment  of  the  turbofan. 

As  you  can  see  (see  Figure  25),  this  engine  has  a  nine-stage  axial-flow  compressor 
with  a  single  centrifugal  compressor  stage.  That  assembly  was  driven  by  a  turbo-row 
turbine,  which  also  drove  a  front  fan  -  a  two-stage  front  fan  -  via  a  gearbox.  This 
diagram  is,  of  course,  very  diagrammatic.  It's  a  great  pity  that  we  were  never  able 
to  complete  that.  It  was  about  half  done  when  the  contract  was  cancelled.  We  had  the 


Figure  25.  A  Simplified  Turbine-Fan  Arrangement 
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wretched  experience  of  being  told  that  we  were  to  concentrate  on  research  and  development. 
We  were  no  longer  to  be  allowed  to  design  and  build  engines. 

[Figure  26  shows]  another  arrangement  by  which  we  avoided  the  complications  of  one 
shaft  going  through  another.  Here  we  had  four  simple  jet  engines  arranged  round  the 
shaft  of  a  fan  assembly,  [with]  the  turbine  driving  that  fan. 

So,  gentlemen.  I'll  tell  you  one  or  two  more  amusing  stories  before  we  go  on  to 
questions. 

When  that  experimental  airplane  first  flew,  of  course  it  caused  an  awful  lot  of 
surprise  because,  though  you  are  used  to  seeing  airplanes  without  propellers,  no  one 
had  ever  seen  that  sort  of  thing  in  those  days.  One  story  that  came  back  to  me  was  when 
two  officers  were  watching  at  Cranwell,  where  we  did  our  first  flight  trials.  One  said 
to  the  other,  "How  the  hell  does  that  thing  work,  John?"  And  the  other  one  replied,  "Oh, 
it's  easy,  old  boy.  It  just  sucks  itself  along  like  a  Hoover." 

And  at  a  neighboring  Air  Force  station,  Waddington,  which  only  about  six  miles  from 
Cranwell,  one  officer  was  sitting  in  the  mess  and  he  said,  "You  know,  chaps,  I've  seen 
a  very  strange  airplane.  It  made  the  most  peculiar  noise  and  there  was  something  very 
wrong  with  it,  and  I  can't  think  what  it  was."  And  then  he  said,  "Oh,  my  God,  I  must 
be  going  round  the  bend.  It  didn't  have  a  propeller." 

One  more  story  I'll  tell  you.  You've  probably  heard  it. 

One  of  the  early  American  jet  airplanes  was  the  P-80,  which  looked  rather  like  a 
duck,  the  Lockheed  P-80,  and  the  story  is  that  two  ducks  were  flying  along  when  the  P- 
80  went  whistling  past  at  500  miles  an  hour.  One  duck  said  to  the  other,  "Say,  that 
guy  is  sure  in  a  hurry  to  get  someplace."  And  the  other  duck  said,  "I  guess  you'd  be 
in  a  hurry  to  get  any  place  if  you  had  a  ten-foot  blowtorch  where  he's  got  it." 

Well  now,  gentlemen,  if  you'd  like  to  ask  a  few  questions.  I'll  do  my  best  to  answer 

them. 

Question:  Why  do  you  figure  it  was  that  although  the  American  version  of  the  Bell 
P59A  was  built  early  during  the  war,  it  was  never  used  in  combat? 

Sir  Frank:  Well,  you'd  better  ask  your  people  about  that.  I  just  don't  know  the 
answer  to  that  one. 

Question:  What  was  it  like  when  you  took  off  in  a  jet?  How  did  it  handle? 

Sir  Frank:  What  was  it  like?  Oh,  delightful.  It  was  far  simpler  t*-«n  today.  On 
the  way  here  I  was  in  the  cockpit  of  the  DC-10,  which  is  a  very  nice  airplane  from  the 
passenger's  point  of  view,  but,  my  God,  is  the  cockpit  complicated!  And  didn't  I  offer 
to  land  it  for  them.  But.  .  .in  the  days  of  the  Meteor,  all  you  had  was  the  throttle 
and  just  one  or  two  simple  instruments  -  the  air  speed  indicator,  the  altimeter,  the 
attitude  indicator,  or  whatever  you  call  it  these  days.  The.  .  .what  is  it  called  these 
days?  [Flight  director?]  Ah,  yes.  Well,  it  was  a  very  simple  airplane  to  fly.  Even 
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in  those  days,  with  the  piston  engine,  you  had  to  go  through  the  routine  TMP  FFG,  you 
see,  when  you  were  taking  off.  "X"  meant  trim,  "M",  mixture  control  -  make  sure  you 
were  in  rich  -  "P",  propeller  pitch,  "F",  flaps  -  make  sure  your  flaps  are  in  the  right 
position  -  and.  .  .I've  forgotten  now,  but  I  had  to  remember  it  in  those  days.  But  you 

didn't  have  to  worry  about  things  other  than  how  to  set  the  undercarriage  and  flaps  on 

the  Meteor  -  no  radiators,  no  pitch  control,  no  mixture  control. 

Question:  Sir,  what  did  the  turbodrill  have  to  do  with  the  turbojet? 

Sir  Frank:  What  did  it  have  to  do  with  the  turbojet?  Well,  its  mechanical  engi¬ 
neering  has  a  lot  in  common  with  the  turbojet  -  you  know,  things  like  bearings  and  gears  - 
and  it  has,  of  course,  a  turbine,  and  a  pretty  efficient  turbine  too,  but  it's  a  hydraulic 
turbine.  But  beyond  that,  of  course,  it  is  a  very  different  thing,  naturally,  and  has 
a  very  different  purpose.  But  it  does  drill  about  four  or  five  times  as  fast  as  the 
conventional  method.  Unfortunately,  it's  now  "on  the  shelf"  since  Rolls-Royce  went  bank¬ 
rupt  -  through  a  jet  engine,  incidentally.  That  brought  that  project  to  an  end. 

Question:  What  kind  of  engine  life  did  the  early  engines  in  the  Meteor  have?  How 

often  did  you  have  to  change  them? 

Sir  Frank:  The  official  life  was  125  hours.  The  Germans  were  lucky  if  their’s 
lasted  25  hours. 

The  program  concludes. 
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